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CHAPTER 1 



Introduction 



1.1 Importance of particle physics in thermal background 

The study of elementary particles in vacuum is based on quantum field theory. How- 
ever, there are situations in which it is necessary to consider these particles in a non-trivial 
background consisting of a large number of particles at finite temperature and/or density. 
Such temperatures and densities of significance occur in stellar interiors and were also 
present in the early universe. Another relevant situation occurs in the laboratory during 
heavy-ion collisions. The formalism needed for the study of these situations is provided 
by thermal field theory, also known as finite temperature field theory H, ^. 

From a theoretical point of view (as we will see in the course of the present work), two 
major features are induced by the background medium. These are the occurrence of an 
additional four-vector in the form of the four-velocity of the medium, and the breaking 
of discrete symmetries like C, P and T even when the underlying Lagrangian respects 
these symmetries. As a result, self-energies, vertex functions and other Green's functions 
involve new tensor structures and associated form factors, together with a relaxation of 
some of the constraints restricting the form factors. A third important feature induced 
by the medium (also to be encountered in the present work) is the breaking of flavour 
symmetry. All these features give rise to rich physics which was absent in the vacuum. 
A modified self-energy, for example, gives rise to modification of the dispersion relation 
(energy-momentum relation). This, in turn, among other things, can open up the phase 
space, and allow processes which were kinematically forbidden in the vacuum to take 
place in the medium. 

Let us now turn to some specific illustrations of changes in the properties of elementary 
particles due to a thermal background. Gauge symmetry spontaneously broken in vacuum 
can be restored by finite temperature corrections. The masses of the fermions and gauge 
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bosons of the Standard Model vanish when this happens [Q. In the phase of restored 
chiral invariance, one finds the generation of thermal effective masses of the particles, 
different effective masses of left-handed and right-handed leptons, and different dispersion 
relations of transverse and longitudinal modes of gauge bosons 0, Modification of 
photon dispersion relations leads to Cerenkov radiation (emission of a photon) by even 
a single neutrino in a medium The dispersion relation of neutrinos become flavour 
dependent in a medium containing electrons but not the other charged leptons [El M. 



In such a flavor asymmetric background, the electromagnetic |jTO| and gravitational [|TT|] 
interactions of the neutrinos in response to external fields are also flavor dependent, while 
the radiative decay of a massive neutrino is strongly enhanced compared to the decay 
in vacuum (as the leptonic GIM mechanism becomes inoperative) [|12|. When neutrinos 
travel in a background of electrons with a static external magnetic field also present. 



the latter causes an anisotropy of the neutrino dispersion relation [|T^. An entirely new 
property of the Standard Model neutrino acquired in a background of electrons is an 
effective electric charge A CPT-asymmetric background medium causes rotation of 



the plane of polarization of light passing through it |T^, a phenomenon known as optical 
activity. 

There are a number of extremely important applications of background-induced parti- 
cle physics, some of which we now mention. The temperature dependence of the effective 
potential for the Higgs field gives rise to the inflationary scenario in the early universe 
|T^. The flavour dependence of the neutrino dispersion relation in matter leads to the 



MSW effect, with consequences for the solar neutrino problem ||I6[. The decay of a pho- 



ton propagating through a medium, known as plasmon, into neutrino-antineutrino pair is 



the major cooling mechanism in white dwarfs and red giants [T^, [T3]. The anisotropy of 
neutrino dispersion relation in matter with background magnetic field has been suggested 
as a possible explanation for the large birth velocities of pulsars . 



1.2 The problems considered in this thesis and their motivations 

The major problem considered in this thesis is the gravitational interaction of charged 
leptons in a medium. Gravitational interactions (in the vacuum) are universal in the sense 
that the ratio of the inertial and the gravitational masses of any particle is a constant. This 
fact, expressed in the form of the equivalence principle, is one of the basic axioms of the 
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general theory of relativity. Although this is a feature of the theory at the classical level, it 
has been shown by Donoghue, Holstein and Robinett (DHR) |^0[ that the corresponding 
linearized quantum theory of gravity also respects this ratio, at least to 0{a). 



However, it was also shown by DHR [20, 21] that this property is lost when the particles 
are in the presence of a thermal background rather than the vacuum. To arrive at this idea, 
the inertial and the gravitational masses must be defined in the context of quantum field 
theory. We consider in Chapters ^ and |^ their precise definitions in terms of the particle 
propagator and the gravitational vertex, which we will need in the subsequent work. For 
the moment, let us denote the mass in vacuum by m, the inertial mass by M and the 
gravitational mass by M', and summarize the results of DHR. These authors considered 
the case of the electron, and calculated Mg and M'^ to 0{a), in a background with a 
temperature T ^ rrie and zero chemical potential. Thus, the background contained only 
photons, but not electrons or any other matter particles. Although in those calculations 
only the case of the electron was considered explicitly, the results are equally applicable 
to other charged leptons, namely, the muon and the tau lepton. Thus for any charged 
lepton i, the DHR results imply 

«7rT^ , , 

Me = me + - , (1.1 

6m£ 

M = me-- , (1.2) 

Sme 

Mi - I ^"^^^ 



^-^^ (1-3) 



to 0{a). Therefore, not only the inertial and gravitational masses of a given charged 
lepton cease to be equal when the background effects are taken into account, but, in 
addition, the ratio of these two quantities is no longer the same for all the particles; i.e., 
universality is lost as well. 

Now, the DHR results, as mentioned before, were obtained by considering only pho- 
tons in the thermal background. But in a matter background with a non-zero chemical 
potential (such as the Sun), one can no longer neglect the matter particles in medium, 
even when T -C me- It is then possible that the matter contribution to the charged lepton 
masses will dominate over the photon background contribution (the latter being, as the 
above results indicate, quite small for T <^ me). Moreover, since a medium with non-zero 
chemical potential is asymmetric with respect to charge conjugation, the matter-induced 
corrections will not be the same for the corresponding antileptons. 



Motivated by these considerations, in this work we calculate the leading matter- 
induced QED corrections to the inertial and gravitational masses of charged fermions 
in a medium that consists of a photon background and a matter background of electrons 
and nucleons ||2^. These represent the dominant corrections for charged leptons and 
antileptons (weak interaction corrections being negligible in comparison). For strongly 
interacting particles such as the quarks, gluon exchange corrections are expected to be 
even stronger and our results will not apply. After finding the general expressions for the 
matter-induced corrections in a generic matter background, we give explicit formulas for 
the corrections in terms of the macroscopic parameters of the background medium for a 
few special cases of the background gases. 

While these calculations are performed for low temperatures and densities, in the last 
part of the thesis, we turn to a background at high temperature. We consider high- 
temperature QED with exact chiral invariance and at zero chemical potential, and inves- 
tigate gauge independence of th one-loop electron dispersion relation in this situation. 

It is well-known that the one-loop electron dispersion relation is gauge independent to 
leading order in the temperature T. This was demonstrated by Weldon by showing 
that the leading part of the one-loop electron self-energy (the self-energy being the field- 
theoretic input to the dispersion relation) goes like and is gauge independent, while 
the gauge dependent part goes like T. Now, a somewhat similar situation occurs for 
the one-loop self-energy of neutrinos in material medium, which is gauge independent 
at 0(1/M^) but gauge dependent when the 0(1/M^) terms are included. However, it 
has been shown |§ that the one-loop dispersion relation following from the latter is still 
independent of the gauge parameter. This raises the possibility that the gauge dependence 
of the part of the one-loop electron self-energy subleading in T may also not show up in 
the one-loop dispersion relation. This possibility is explored in the present work in two 



different limiting cases One of theses cases is the limit of momenta much larger than 
eT. The other is the zero momentum limit, when we investigate the gauge independence 
of the effective mass. 

It may be mentioned in this context that the demonstration of gauge independence of 
the properties of elementary particles in material medium has turned out to be of extreme 
importance in thermal field theory. In particular, the gauge dependence of the gluon 
damping rate at one loop was a long-standing problem. In Ref. [^, the Ward identities 
determining the gauge dependence of the gluon dispersion relations were deduced and 



4 



used to prove gauge independence in a self-consistent perturbative expansion. In Refs. 



25[| and it was shown that such a consistent expansion requires resummation of 



the hard thermal loops. The use of this idea finally led to a gauge-independent gluon 



damping rate to leading order in the coupling constant |27|. It will be seen that the gauge 



dependence identities of Ref. [Q, although dealing with gluons, are still relevant for the 
case studied by us. 

1.3 Outline of the thesis 

In Chapter we introduce the real-time formulation of thermal field theory. After 
working out the thermal averages of bilinears of ladder operators, we deduce the time- 
ordered thermal propagators for various fields. Additional vertex and propagators in the 
real-time formulation are then briefly mentioned. 

In Chapter we discuss the linearized theory of gravity. The linearization of the 
metric and how it gives rise to the graviton are explained. Photons and spinors in curved 
space are discussed, leading to the QED Lagrangian of curved space. This Lagrangian 
is then linearized, and the Feynman rules for the lowest order interaction of the graviton 
with photons and fermions are deduced. 

In Chapter ^, we consider charged leptons in thermal QED. General expressions for 
the dispersion relation, the inertial mass, the spinor and the wavefunction renormaliza- 
tion factor are obtained at O(e^) for a fermion as well as an antifermion. Relations 
between the fermion and antifermion properties are then deduced for both C-symmetric 
and C-asymmetric backgrounds. Explicit O(e^) calculation of the inertial masses and 
wavefunction renormalization factors of charged leptons and antileptons are performed. 

In Chapter |^, we turn to gravitational couplings in a medium. We show in detail 
how the gravitational mass is determined from the gravitational vertex function, and 
obtain an operational formula for the gravitational mass at O(e^). We then calculate the 
gravitational vertex of the leptons, and use it to find the gravitational masses of charged 
leptons and antileptons in a medium. Since the terms involving the fermion distribution 
function cannot be evaluated exactly, we evaluate the corrections in two different limits, 
viz., the classical and the degenerate limits for the electron gas. 

In Chapter ^, we consider some aspects of high-temperature QED. The expression for 
the one-loop self-energy of the electron in a general linear covariant gauge is deduced. 
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The gauge dependence relation for the one-loop electron dispersion relation is obtained, 
and gauge independence is proved for momenta much larger than eT. After obtaining a 
formula for the effective mass of the electron, gauge independence of the one-loop effective 
mass is investigated. 

In Chapter ^ we present our conclusions. 

1.4 Some notations and conventions 

We now summarize some notations and conventions which will be frequently used in 
this thesis. Although these will be mentioned again in the appropriate places, it was felt 
that a separate summary of them may be useful as a ready reference in going through the 
present work. 

Throughout this work, we will denote four-momenta by lower case letters, and the 
corresponding three-momenta by upper case letters. For example, we write = {p^, P). 
The magnitude of P is denoted by P. Also, we use the notation Ep = V-P^ + 'rn'^, where 
m is the mass of the particle concerned. 

The mass in the vacuum, the inertial mass in a medium and the gravitational mass in 
a medium are denoted by m, M and M' respectively. They come with various subscripts: 
/ denotes a generic fermion, p a proton, i a charged lepton and £ a charged antilepton. 
Further, i can be e or /i or r. 

The notations for the thermal corrections to m, as defined in Eqs. ( [4. 1041 ) and ([ 



contain the additional subscripts 1 and 2. Here 1 stands for the contribution from the 
photon background, and 2 stands for the contribution from the electron background. The 
symbols 1 and 2 are used in this way for various other medium effects also. Thus, it 
is used to denote the contributions from a diagram: the Bl contribution, for example, 
denotes the photon background contribution to the diagram B, and the A2 contribution 
denotes the electron background contribution to the diagram A. 

Finally, in our work on high-temperature QED, the effective mass of the electron is 
denoted by M, while Mq denotes the effective mass to leading order in temperature. 
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CHAPTER 2 



Thermal field theory 



2.1 Introduction 

Thermal field theory is the theoretical framework for studying a large number of 
quantum mechanical and relativistic particles in thermodynamic equilibrium. The two 
well-known versions of this subject are the imaginary-time formalism ^ ^ |2^, ^ and 



the real-time formalism 0, ^, |2^, ^ , of which we are going to discuss only the latter. 
The real-time formalism has the following advantages. The calculations are manifestly 
covariant, the effects of the medium are separated out at the very beginning from the 
part already present in the vacuum, and it can be easily understood that the medium 
does not introduce any new ultraviolet divergence. The particular approach to the real- 
time formalism to be discussed by us is the canonical approach This approach 
involves a generalization of the canonical quantization formalism of ordinary quantum 
field theory, by defining the Green's functions as thermal averages instead of vacuum 
expectation values. 

2.2 Thermodynamic equilibrium 

Thermal field theory makes use of quantum field theory as well as statistical mechanics. 
As mentioned in Sec. |2.1| , we are going to use the canonical quantization formalism of 
quantum field theory. For doing statistical mechanics, we have to choose a particular 
ensemble according to our convenience. In the present case, each of the systems making up 
the ensemble is in some state of the Fock space, and the calculations are the simplest (see 
below) if the ensemble contains all the states of the Fock space, without any restriction on 
the total number of particles in a system. Consequently, we consider the grand canonical 
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ensemble, which is characterized by the temperature 

T^i (2.1) 

and the chemical potential /i. (We have chosen our units such that the Boltzmann constant 
is equal to unity.) The thermal expectation value (ensemble average) of any operator O 
is 

{O) = TiipcO) , (2.2) 
with the density operator pc in the grand canonical ensemble given by 

P-/3J 

where 

J = H~pQ. (2.4) 

Here H is the normal-ordered hamiltonian, and Q the normal-ordered charge operator, 
i.e., the total number operator for the particle minus that for the antiparticle. 

We shall now consider the case of free fields. Let us begin with a real scalar field, for 
which (and, in general, for any self-conjugate field) = 0. For the sake of convenience, we 
shall consider quantization with discrete momenta in this section. (From Sec. P]B| onward. 



we shall consider quantization with continuous momenta; Appendix ^ relates the two 
approaches.) We therefore have 



H = Y.EKa^.afi. (2.5) 

K 

Here 

Ek = Vi^2 + ^2 ^ (2.6) 

m being the mass of the quantum of the scalar field, and a^a^ is the number operator 
for quanta of momentum K. 

We are now going to find out the thermal average of the bilinears of the ladder oper- 
ators. First consider 

Tr[atap exp(-/? Exa^^ag-)] 
alap) = L (2.7) 

K 



where Eqs. (2.2)- (|2!5| ) (with /i = 0) were used. The trace is in the Fock space, which is 



spanned by the direct product of the eigenstates of the number operators a^^a^ for the 



various modes K. Therefore for any operator O 



(2.8) 



Kl K2 



Here {nj}^\{nj^J ■ ■ ■ denotes the direct product of the states (?^^J, {^j^J, ■ ■ ■ ■ Also, since 
we are considering a bosonic field, the sum over each eigenvalue n^^ runs from zero to oo. 
Use of Eq. ( |2.8D in both the numerator and denominator of Eq. ( ^.7[ ) leads to 



K 



^-f3Epnp 



d 



(2.9) 
(2.10) 



The sum over np is an infinite geometric series, and equals 1/(1 — e ^^p). So 



a}pap 



ol3Ep 



f2.111 



which is the Bose distribution function. Using 



ap,ap 



1, this leads to 



apa'^p 



1 + 



Also using Eqs. ( p.2|) -( ^3D and Eq. (|2.^ ), it follows at once that 



f2.12) 



alap,) = (apal) = for P ^ P' 



(apap,) = (^alal) = 



(2.13) 
(2.14) 



For a complex scalar field with the number operators o^Oj^ for the particle and b^^bj^ 
for the antiparticle of momentum K, Eq. ( p.4|) gives 



J = J2[iEK - f^)aUg + {Ek + fi)blb 



(2.15) 



K 



The Fock space states are now the direct product of the eigenstates of a^a^ and 

for various K. The only change in Eqs. ( p.ll| ) - ( p.l4| ) is then found to be Ep — > Ep — /i. 



For the bilinears of the antiparticle ladder operators, relations similar to ( 2.11 ) - ( 2.14|) 
hold with Ep Ep + fi. Also, averages of mixed operators, like {apbp,), are all zero. 

We next turn to the case of the Majorana fermion field, which is analogous to the 
real scalar field, and has /i = necessarily. The important difference from the case of 
the real scalar field relevant for our purpose is that the number operator ^ for any 

given mode can only have the eigenvalues and 1. Consequently, the geometric series 
represented by the sum in Eq. ( [2.1(J| ) now contains only two terms and equals 1 + e~^^^. 
So 

4/p.) = ;^l^' (2-16) 



which is the Fermi distribution function. Using 



1, this leads to 



The other thermal averages are zero, as for the real scalar field. 
Finally, for a fermion field which is not self-conjugate, 

J = j:[iEK - /x)ct ^c^,. + {Ek + f^)d^g/K,s] (2-18) 

K,s 

(analogous to Eq. ( p.l5| )), so that the replacements Ep — ^ Ep =f /i result, just as in the 
case of the complex scalar field. 

2.3 Time-ordered thermal propagators 
2.3.1 Scalar field 

Time-ordered thermal propagators are defined by replacing the vacuum expectation 
value in the definitions of the vacuum propagators by a thermal average. Beginning with 
scalar fields, we shall first work out the more complicated case of the complex scalar field 
in detail, and then state the result for the case of the real scalar field. We shall use the 
Fourier expansion 



(a(P)e-^P-" + 6VP)e^P-"), (2.19) 
•/(27r)32Ep^ ^ 



for the complex scalar field, where p'^ = {p^, P) and = Ep. Since we are now considering 
continuous momenta, the thermal averages read 

[a\P)a{P')) = 6'^'\P-P')fB{Ep,fi), (2.20) 
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aiP)a\P')) = 5^'\P-P')[l + fB{Ep,i^)], (2.21) 
b\P)b{P')) = 6^'\P-P')fB{Ep,-f,), (2.22) 
b{P)b\P')) = S^'\P-P')[l + fB{Ep,-fi)], (2.23) 



with the other thermal averages equalhng zero. Here the Bose distribution function for 
the particle is given by 

fB{Ep,f^) = (2.24) 

Eqs. (g) and (^^, consistent with Eq. (|X1| ), correspond to Eqs. (|2Tl|)-(|27[3|). 
The time-ordered thermal propagator for the complex scalar field is defined by 

tD{x-y) ^ {T[<Pix)<p\y)]) (2.25) 
= ^(xo - yo){m4>\y)) + 0{yo - xo){(P\y)(P{x)) ■ (2-26) 

Using Eqs. ( p.l9| ) - (|2.23|) , we can obtain {(j){x)(t)\y)) and (0^(?/)0(x)). After substituting 
these in Eq. ( p.26| ), let us spht up the resultant expression as 

iD{x -y)= iDf{x -y) + Dt{x - y) , (2.27) 

where iDp{x — y) is the part independent of /b, and Dt{x — y) is the part arising out 
of the medium. Now, if were set equal to zero in Eqs. ( |2.20| ) - (|2.23|) , the thermal 
averages would reduce to the vacuum expectation values, and the thermal propagator to 
the propagator in the vacuum. So, we can immediately write down 



Df{p) = ^y— -. (2.28) 



On the other hand, 



r d^P 

Dt{x -y) = j j^-^^[f^(Ep, /x)e-^^ (--^) + fsiEp, -^^)e^^<^-y^] . (2.29) 

(One can also get Eq. (p.29|) by using Eq. ( |A.9D and the discretized versions of Eqs. 
( p. 201 ) - (|2.23| ), and finally using Eq. ( |A.4D in reverse.) Now, Eq. ( p.29|) can be rewritten 

as 

Drix -y) = I (2^')32Ep " /i)e-^P«(^°-^°)+^^-(^^-^^) 

+5{po + Ep)fB{-po, -;i)e-*Po(^'«-2^o)-*^ (^-27)] . (2.30) 
11 



Changing P to —P in the second term, and using 

1 



we obtain 



2Ef 



6{poTEp) 



12.31) 



Dt{p) = 2'k6{p^ - m^)r]B{p) , 
Vb{p) = 0{po)fB{po, /i) + 0{-po)fB{-po, -/i) • 



(2.32) 



Thus, tlie tliermal propagator is given by Eqs. ( |2.27| ), ( p.28| ) and ( |2.32| ). However, 
this form is not manifestly covariant since the distribution function ( ^.24[ ) was written 
in the rest frame of the medium. Therefore we introduce the four-velocity f ^ of the 
centre-of-mass of the medium. In the rest frame of the medium, 



(2.33) 



Then carrying out the replacement 



Po-^p-v 



(2.34) 



in Eq. ( p.32| ), we arrive at the manifestly covariant form of the thermal propagator. 
For the real scalar field, we have to use the Fourier expansion 

d^P 



(f){x) 



:(^a{P)e-'P'' + a\P)e'P'' 



(2.35) 



'{27t)^2Ep 

and the thermal averages (|2.2CI| ) and ( p.21|) with n = 0. This leads to the thermal propa- 



gator 

(T[0(a;)0(2/)]) 



d^p 
(2^ 



-ip-(x-y) 



p2 _ ^2 _|_ 



+ 2n6{p'^ - 



ol3\p-v\ _ I 



, (2.36) 



which is seen to be the thermal propagator of the complex scalar field at /i = 0. 



2.3.2 Fermion field 

For the fermion field, we shall use the Fourier expansion 
r d^ P 

H^) = J ^^E(c.(^)«.(i^)e-^^-^ + ciI(P)t;.(P)e^^-^ 



(2.37) 



12 



and normalize the spinors as 



(2.38) 



(This normahzation will be later extended to the spinors in a medium; see Eqs. ( ^.27| ) 
and (|4.4(j ).) The thermal averages then read 



ci{p)cAP': 
cs{P)cUp': 

di{p)dAP': 
ds{p)dUp': 



Sss'5^'\P-P')ff{Ep), 
S,,,6^^\P-P')[l-ff{Ep)] 
6ss'6^'\P-P')ffiEp), 
6ss'S^'\P-P')[l-ff{Ep)] 



(2.39) 
(2.40) 
(2.41) 
(2.42) 



with the other thermal averages equalling zero. Here the distribution functions for a 
fermion and an antifermion are given by 

ffj(E-) = ,,iEX) + 1 ' (2.43) 

respectively. (We have denoted the fermion chemical potential by /ij, since we will later 
deal with several species of the fermion /.) Eqs. ( |2.39| ) and p.4CI| ), consistent with Eq. 
( |CT| ), correspond to Eqs. (|23^ ) and ( pT7|) . 

The thermal propagator for the fermion is defined by 



iSfabix) = {T[ilJa{x)'ilJb{y)]) 

= 0{.xq - yo){^aix)^IJb{y)) - Oivo - xo)(V'fe(i/)V'a(a;)) 



(2.44) 
(2.45) 



Using Eq. ( p.37|) , Eqs. ( p.39|) - (|2.42|) , and the spin-sum relations 

S 

J2v,iP)vsiP) 



2Ep 

Ep-yo — P ■ ^ — nif 



2Ei 



(2.46) 



following from Eq. ( p. 381) , we can find the R.H.S. of Eq. (|2.45|) . The expression is then 
split up into the vacuum part and the thermal part: 



Sf{x-y) 



(2vr)4 

iSf{p) = iSpfip) + Srfip) ■ 



(2.47) 
(2.48) 
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The vacuum part is 

— mj + le 

while, following the derivation in the case of the complex scalar field, the thermal part is 
found to be 

STf{p) = -27r{^ + mf)6{p^ -m})7]f{p), (2.50) 
Vfip) = ^ / — + J , — . 2.51 

2.3.3 Photon field 

In this case, the thermal propagator in the Feynman gauge is the same as that for 
the real scalar field, except for the facts that the photon is massless, and that there is an 
extra factor of —rji^u (as in the vacuum). Thus, we have the propagator 

iD^,{k) = -T]^, [iAp{k) + Arik)] , (2.52) 

where 

A,(*) = (2.53) 

Arik) = 27rS{k^)7]^{k) , (2.54) 

The factor r]^{k) arises out of the photon distribution function 

/7(^) = > (2-56) 

where K is the magnitude of the three-momentum of the photon. 

We end this section with the following observation. At any finite temperature, the 
distribution function contained in the medium-dependent part of the thermal propagator 
becomes an exponential damping factor for momenta much larger than T. Consequently, 
the corresponding loop integration will have an effective cut-off of 0{T). If T = 0, we can 
still have a Fermi gas; in that case, the Fermi momentum provides a cut-off. Thus, the 
medium does not introduce any new ultraviolet divergence, and rcnormalization in the 
vacuum suffices to make the theory finite in the presence of thermal background. 
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2.4 Additional vertex and propagators 



The full structure of the real-time formulation of thermal field theory is more compli- 
cated than the replacement of vacuum propagators with time-ordered thermal propaga- 
tors. However, as will be seen later, the additional features of the full structure do not 
play a role in the present work. So, in this section, we shall just outline the main results 
involving those additional features |]32|| . 



Ultimately, one has to calculate an n-point Green's function like (T[0(xi)0^(x2) 
■ ■ ■ (j){xn)]), where (plxi) are the Heisenberg picture fields. (We shall consider only complex 
scalar fields in this section. Similar results hold for other fields.) In a medium, as in the 
vacuum, one can prove Wick's theorem and reduce the n-point function to a sum of fully 
contracted terms. But the fact that the Green's function in the present case involves a 
thermal average over all the Fock space states, rather than a vacuum expectation value, 
leads to four types of contractions or "propagators." They are defined as 

tD^,ix-y) = (T[0(x)0t(^)]), (2.57) 

tD,2ix-y) = (T[0(x)0t(^)]), (2.58) 

iDu{x-y) = (0t(i/)0(x)), (2.59) 

tD2i{x-y) = (0(a;)0t(y)). (2.60) 



Here T stands for anti-time ordering. In Eqs. ( p. 57] ) - (|2.60|) , the fields are the interaction 



picture fields, which can be expressed in terms of the ladder operators. Thus, in addition 
to the time-ordered propagator of Sec. we also have anti-time ordered propagator, as 
well as propagators without any ordering in time. As given by Eqs. (|2.27| ), ( p.28| ) and 

O) 



tDuip) = + 27r5(p2 _ m')r]Bip) . (2.6i; 



p2 _ jyi^Z _|_ 



The others are found to be 



iD22{p) = 5 - + 27r(5(p^-m^)r/B(p), (2.62) 

iD^2{p) = 27i6{p^ -m^)[r]B{p) + 9{-p-v)], (2.63) 
iD2i{p) = 27t6{p^ -m^)[7]B{p) + e{p-v)]. (2.64) 



Here it is to be understood that the replacement (|2.34|) has been carried out in ?7_b(p). 
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The Feynman diagrams now contain two types of vertices: the 1-type vertex is the 
usual one, while the 2- type vertex has a vertex factor which differs from that for the 
1-type vertex by a sign. Green's functions can now be evaluated with these two types of 
vertices and four types of propagators. However, the external legs representing physical 
particles can only be connected to 1-type vertices. 

The importance of this doubling of degrees of freedom manifests itself in diagrammatic 
calculations in the following way. If one uses only the 1-type vertex and "11" propagator, 
one sometimes ends up with singular expressions involving products of delta functions 
with the same argument. However the use of the additional vertex and propagators leads 



to the cancellation of such singularities |33 
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CHAPTER 3 



Linearized theory of gravity 



3.1 Introduction 

Gravity is governed by the general theory of relativity which tells us how the geometry 
of space-time, as given by the metric tensor g^^i^x), is affected by the presence of matter. 
The Einstein field equations for g^j^y in the absence of external source, derived from the 
action for pure gravity, are 

- \R9^.u = 0. (3.1) 

Here -R^^^ is the Ricci tensor and R is the curvature scalar. One important point is 
that these equations are non-linear in gfj_^. This is in contrast with the linear free- field 
equations for the scalar, fermion and photon fields, where we use the principle of linear 
superposition to write down the general solution before quantizing the field. The reason 
for this non-linearity in the case of gravity is that the gravitational field possesses energy 
and hence effective mass, and therefore produces gravitational field. The gravitational field 
thus contributes to its own source. To ignore this feedback effect, we have to consider 
a weak gravitational field. One thus considers only small deviations from the flat space 
metric, and writes 

Qf^uix) = 7]^^ + 2Kh^^{x) (3.2) 

with \Kh^i^\ <^ 1. Here k is related to the Newton constant G. The approximation stated 
in Eq. ( p. 21 ) is the starting point of the linearized theory of gravity [Q, Now r]^^, and 



hence from Eq. (|3.2|) , /i^j^, are symmetric flat space tensors, but not coordinate tensors (in 
contrast with g^,y, which is a symmetric coordinate tensor, and SO, cLS cL special case, also 
a tensor in fiat space). Consequently, a// tensors in an expression should be regarded only 
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as flat space tensors after we use the approximation of Eq. (3^) to perform an expansion 
of the expression in powers of k, and we will adopt the convention of raising and lowering 
the indices of such tensors with the fiat space metric. 

The precise value of k, to be deduced in the next section by identifying h^i, with the 
graviton field, is given by 



K 



VmG. (3.3) 



Since v G = 1 / mpianck in natural units, it is clear that the weak field approximation 
given by Eq. ( p.2|) is valid when all mass-scales in the problem are much smaller than the 
Planck mass. These mass-scales include not only particle masses, but also temperature 
and chemical potential, as we shall later perform calculations in the presence of a thermal 
background. Since mpianck is 0(10^^) GeV, the weak field approximation holds in all 
physical situations except in the very early universe. 

At this point, it is interesting to consider the case of the Yang-Mills theory for com- 
parison, because the equation of motion for the non-Abelian gauge field derived from the 
pure Yang-Mills Lagrangian is also non-linear. However, there the equation of motion 
is linear when one considers the free gauge fields, setting the coupling constant equal to 
zero, and so the need to linearize does not arise. In contrast, the fiat space metric tensor 
rj^i,, being a constant, does not have any dynamics (and the LHS of Eq. ( p.l|) reduces to 
zero in this case). So, to do quantum field theory with gravity it is essential to impose 
the condition (|3.2|) , which considers linear deviations from the fiat space metric case. 



3.2 Graviton 

The field equations ( p.l| ) follow from the Einstein-Hilbert action given by 

cf^x if, (3.4) 



The curvature scalar, the Ricci tensor and the affine connection are respectively given by 

R = g^^R^,, (3.6) 

rj, = ^g^^id^g.^ + d^g^^-d^g,,), (3.8) 
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while 



g = detg^ 



flU- 



(3.9) 



In anticipation of the transition to quantum field theory, here and elsewhere we always 
write the action in the form ( |3.4| ) (although =Sf is not a coordinate scalar). 
From the weak field condition (13. 21), we obtain 



(3.10) 



(since g^'^Qux = S'^), and 



-g = 1 + ^/1 + 0(^2) 



(3.11) 



where 



(3.12) 



Note that we are raising and lowering the indices of the graviton field with the fiat space 
metric [^]; this is in accord with the convention stated in Sec. [O . 
Use of Eqs. (|3) and {^J^) in Eq. (gj) leads to 



K 



(3.13) 



Putting this in Eq. p.7|), one can determine i?^,^ to 0{k'^). Let us write 



(3.14) 



with R^u^^^ and R^u^'^^ denoting the parts linear and quadratic in k respectively. 

Following again the convention stated in Sec. |3.1| , we shall now raise and lower the 
indices of Rfj,u^^\ Rfj.u^'^^ and d/dx^ with the fiat space metric Then 



R%^'^ + KhR%^'^ - 2Kh^'^R^^^^ + + 0{k^). (3.15) 



The first term on the RHS is the only term of 0{k). In fact, 



(3.16) 
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and this, being a surface term, does not contribute to the Einstein-Hilbert action. So we 
now turn to the three 0{k'^) terms on the RHS of Eq. ( |3.15| ). First consider R'^^^'^K After 
leaving out surface terms, we obtain 

+ {d,h){d^h-2dxh^^)]. (3.17) 

Carrying out of several partial integrations with discarding of surface terms [for example, 
{dxh^^){d^h^^^) ^ -h^^{^x^^h^^^) {d'^K^){dxh^^^)] leads to 

= -n\d^h,,x){d^^h^^) - 2{dxh^^^){d^h^^) + 2{d^h){dxh^^) - {d^h){d^^h)]. (3.18) 

Similarly one can write KhW^^^^^ and 2/t/;,^'^i?^j,*-^'' in terms of the combinations contained 
on the R.H.S. of Eq. (|3.18|) . Putting everything back into Eq. ( p.l5| ), we get 



WttG^ = f^^[{d,hr,xWh^^) - 2{dxh^^Wh^^) + 2{d,h){dxh^^) - {d,h){d^h)]. (3.19) 

Here and henceforth it is to be understood that the expression for 167rG=Sf has been 
written neglecting the 0{hC') terms. 

With the aim of arriving at the propagator, we carry out further partial integrations 
to obtain 

IQtxG^ = K^Kp^'^^^'^h^^, (3.20) 

+7]"^'d^d'') +4r]'"^d^d'']. (3.21) 

We have explicitly symmetrized ^"/^/^'^ with respect to a ^ f3 and with respect to fi ^ u 
(such symmetrization helps one to guess the form of the inverse; see later). But this 
operator does not possess an inverse! This is easily seen from the fact that 

^"^^^a„a^a^9^A(a;) = (3.22) 

where A(x) is any function of x, i.e., the operator has a zero eigenvalue. 

The problem is similar to the one encountered in quantizing the photon field, and the 
solution also is the same, viz., gauge- fixing. To this end, one chooses the harmonic or de 
Bonder gauge 

F^" = d^h"'' - ^d^h = (3.23) 
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and adds to WtcG^ the gauge- fixing term p7 



IGttG^gf = 2k^F^F^. (3.24) 

Before proceeding furtlier, let us dwell a little on gauge invariance. 

While the Einstein-Hilbert action is invariant under general coordinate transforma- 
tions, the allowed coordinate transformations in the present case are only those which 
sustain the weak field approximation. Thus, we consider 

_^ ^'i^ = x'' + 2/te^(x) (3.25) 



and demand that 



Using these relations in 



g'''" =7]^"" -2Kh'>"' + 0{n^). (3.26) 



we arrive at the following gauge transformation of /i^j^ which keeps the action invariant: 

= V - (^M^i^ + ^^^^m)' (3-28) 

where = rj^^e'^. It can easily be checked that if h^u does not satisfy the gauge condition 
stated in Eq. ( ^.231 ), /i^^ does, provided 

□e'^ = F^". (3.29) 

So it is always possible to implement the harmonic gauge condition. Another point to be 
confirmed is that the gauge-fixing term given in Eq. ( p.24| ) does break the gauge invariance 
of the action. This readily follows from 



9' 



/^^r^^ = -4kF^ + 0{k^). (3.30) 



Since the affine connection is not a coordinate tensor, the square of the L.H.S. is not 
a coordinate scalar. (Note that ^/qf involves d^x which also is not a scalar. But the 
transformation of d'^x cannot compensate the breaking of gauge invariance, as IGttGJ^gf 
is already of 0(fi;^).) 
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We now explore the consequences of fixing the gauge. Addition of the gauge non- 
invariant term given in Eq. ( |3.24| ) immediately yields the correctly normalized kinetic 
terms for the massless spin-2 field, to be identified as the graviton field: 

^ + ^GF = \{d^Kx){d^h-^) - -^{d^h){d^'h), (3.31) 

provided we set k equal to the value stated in Eq. ( ^.31 ). 

We conclude this section by obtaining the graviton propagator. Partial integration of 
the last expression leads to 

s^ + s^GF = J d'x hi^p^^^^'h^,, (3.32) 

^<^0^^- = _l(ry"/^/- + ry"-//^ - r/"^r/'^'^)n. (3.33) 

So the graviton propagator iD^p^^u is i/^^ times the inverse of 

M^^^^"^ = -(r^"'^/^ + 7]"^//'^'^ - r/^^r]^'^). (3.34) 
The inverse must be determined from 

= 1(5^5; + 5^^5\). (3.35) 

Guided by the fact that must be symmetric under a <-> /3 and (separately) under 
A ^ p, one can write down 

{M~^)c,i3.xp = A{r]axrit3p + VapV^x) + BrjaisVxp- (3.36) 
We then obtain A = 1/2 = —B. Therefore the graviton propagator is given by 



li^aPM - i ^ • (3.37) 



One should now compare the photon propagator in the Feynman gauge, given by iDp^ = 
—iViJ.^ 1^'^ ■ Clearly Dqo is negative, while -Doo,oo is positive. This has the consequence that 
the fermion-fermion interaction in the non-relativistic limit is repulsive when mediated 
by the photon and attractive when mediated by the graviton |3^ . 



3.3 Photons in curved space 

For the electromagnetic field, the Lagrangian in curved space is given by 

= -\^F,,F^pg^''g''^ (3.38) 
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where 



F^^u = d^A, - d,A^. (3.39) 

This has been written down following the standard prescription of passing from flat space 
to curved space, i.e., by replacing all flat space tensors with coordinate tensors, all or- 
dinary derivatives with covariant derivatives (by which we shall mean gravity-covariant 
derivatives) and the flat space metric tensor with g^j_y. The additional a/— g in Eq. (|3.38| ) 



comes because df^x in the action gets replaced with the coordinate scalar d^x^—g. It may 
be noted that the use of covariant derivative and ordinary derivative in Eq. ( p.39|) are 



equivalent due to the symmetry of the Christoffel symbol in the two lower indices. 

3.4 Spinors in curved space 

The procedure of replacing flat space tensors with coordinate tensors cannot be ex- 
tended to the the case of spinors. The procedure works with tensors because the ten- 
sor representations of GL(4,R), the group of 4 x 4 real matrices, behave like tensors 
under the subgroup S0(3,l). Thus, considering the vector representation as an exam- 
ple, V'^{x') = {dx' ^ / dx^)V''{x) becomes V''^{x') = A^^^V^^x) if we restrict ourselves to 
x'^ = K^yX'^ . But there are no representations of GL(4,R) which behave like spinors 
under S0(3,l) (there is no function of x and x' which reduces to -D(A) for x'^ = A'^^x"). 

Consequently, spinors have to be considered as scalars under general coordinate trans- 
formations. But the problem of incorporating the flat space spinorial property into curved 
space remains. This is tackled by making use of the tangent space formalism, to be de- 



scribed now [^, |4T 



At any given point of the curved space, we can choose a locally inertial coordinate 
system. Let the locally inertial coordinates at the point x'^ be be ^"(x) with a = 0, 1, 2, 3. 
Since the locally inertial coordinate system at each point may be redefined by an arbitrary 
Lorentz transformation (LT), the action must be invariant under the local LT 

C{x)^^'^{x)=AMx)e{x), (3.40) 

with A°'b{x) a real matrix satisfying 

riabA\{x)A'dix) = r],d. (3.41) 
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The action, of course, has to be also invariant under the general coordinate transformation 
x'^ — > x''^, which is independent of the local LT. The requirement of invariance under local 
LT does not affect the actions discussed earlier, since coordinate tensors like ^f^i, and 
are taken to be local Lorentz scalars (i.e. scalars under local LT). 
Let us now consider the vierbein or tetrad, defined by 



dxi^ 



(3.42) 



x=X 

These constitute a set of four coordinate vectors, forming a basis for the (flat) tangent 
space to the curved space at the point x = X. Under x^ — x'^, they transform as 

e'Mx') = ^,e\{x), (3.43) 

since ^'"(x') = ^"(x). Under local LT, the vierbein transforms just like ^"(x): 

e'%(x)=A^(x)e^(x), (3.44) 

since x'^ does not transform. From the definition given by Eq. ( p.42|) , it follows that 
raising or lowering of the local Lorentz index a (in general, a latin index) and the general 
coordinate index fi (in general, a greek index) of the vierbein must be done with rjab and 
g^u respectively. 

The metric g^i, is given by 

g^,Ax) = Vabe%ix)e\{x). (3.45) 

Therefore 5^ = Cb^e^v, i.e., C},^ is the inverse of e'^ ^. We can rewrite this as l^i, = 
(E~^)^feEbi,, where the indices now denote matrix indices. It then also follows that 1^6 = 
Ea^(E~-^)^;,. This means that 5^ = e^^Cb'^, or, 

V''' = g'"'{x)e\{x)eK{x). (3.46) 

Eq. ( p.45| ), or equivalently, Eq. ( p.46| ), expresses the orthonormality of the vierbeins. 

The vierbeins turn out to be necessary for constructing the curved space Lagrangian 
for the Dirac fermion. This will involve modifying the fiat space Dirac Lagrangian 

■i 



(/) 



-ip'-f^df^ip + h..c. —mfifjilj (3.47) 



such that the action is both a coordinate scalar and a local Lorentz scalar. With this aim 
in mind, we turn to the transformation properties of the spinor field. 
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Spinor fields are coordinate scalars which transform under local LT as 



(3.48) 



Here DaplA) is the spinor representation of the Lorentz group. Note that we use a and 
P to denote the components of ip, D{A) and other quantities associated with the spinor 
representation; they must not be confused with general coordinate indices. Eq. ( |3.48|) is 
the relation which incorporates the fiat space spinorial property into general relativity. 

Since A is a function of x'^, d^ipa does not transform like ipa under local LT. Hence 
one defines the covariant derivative 



(3.49) 



with the connection matrix Qn transforming as 



% = D{A)n,D-\A) - {d,D{A))D-\A), 



(3.50) 



so that 



(3.51) 



The connection matrix can be written as 

[n^Uix) = ^[SabUto^'"'ix) (3.52) 

where Sab are the generators of the Lorentz group in the spinor representation: 

1 i 

Sab = -(^ab = -[la,lb] (3.53) 

(with 7a denoting the ordinary gamma matrices), and cj^"* is a coordinate vector field 
antisymmetric in a and b, called the spin connection. 

From Eq. (|3.51|) , it follows that is a local Lorentz vector (just like %l)'^°"ip). It 

is also a coordinate vector (as ip is a, coordinate scalar). So when we contract it with Ca^, 
we get a quantity which is a local Lorentz scalar and a coordinate scalar. We can now 



modify Eq. (|3.47|) to write down the curved space Lagrangian 



(/) 



I - 



-^lj-i''ea^^^,ip + h.c.j - m/V^V 



^ ■ 4 
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(3.54) 
(3.55) 



The latin indices of the spin connection are by definition raised and lowered with the flat 
space metric. 

We now explain why we started with the explicitly hermitian form of the flat 

space Dirac Lagrangian. As long as we stay in flat space, the form ^'q^ ' = ip{i'')^d^—mf)il} 
can also be used, since it differs from ^[^"^ by ^d^{il)'-)^il)) which contributes a surface term 
to the action. But this difference is modifled to ^\/—gea^d^{ip'j""ip) in curved space, which 
is no longer a surface term. Thus we would have missed some contribution if we naively 
modifled ^[/^ ' to make the transition to curved space. 

Our next aim is to express the spin connections in terms of the vierbeins. The proce- 
dure to be adopted by us will involve the covariant derivative of the vierbein. So we flrst 
discuss the covariant derivative of a local Lorentz vector. 

We have already discussed the covariant derivative of a local Lorentz spinor, and 
similar considerations apply for other representations of the Lorentz group. The analogues 
of Eqs. ( p.48|) -( p.52|) for a local Lorentz vector fleld V"" are 



V^ix) V'ix) = A\ix)V\x), (3.56) 



^^V" = d^V^ + [nl]\V\ (3.57) 



K']\ = ^^Kr,i^-')\ - id,A\){A-')\ (3.58) 
%V^{x) ^ A\{x)%V\x), (3.59) 



K]\i^) = l[U\<'i^)- (3.60) 

Here is the connection matrix and Jed the Lorentz group generators in the vector 
representation: 

[^cJ"'' = ^(«-5^5^). (3.61) 

The spin connection is representation- independent and the same as that in Eq. ( pj.52| ). 
Eqs. (^) and ( |33lD yield 

= (3-62) 
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and so Eq. (|3.57|) reduces to 



(3.63) 



For completeness we now give the transformation of uj^j_ab (which are the analogues of 
the gauge fields of Yang-Mills theory) under an infinitesimal local LT, parametrized by 



K\{x) = 5l + e\{x). 



Eq. ( |3.41[ ) then gives the condition 



Oha — —dab- 



(3.64) 



(3.65) 



where Oab = Vac0%- Using Eq. (TO) , Eq. (|33^) and (A"^)";, = K- in Eq. (|38 
obtain 



, we 



One can check that the same transformation is obtained on using Eq. ( p.52|) 



D(A(x)) = 1 - -s,,r^(x), 

and the fact that Sab satisfy the Lorentz algebra 

[Sab-, Scd\ = i{f]bcSad — VacSbd — VbdSac + VadSl 



be 



(3.66) 



(3.67) 



(3.68) 



in Eq. ( 3.50 ). (It is necessary to make ui'^^^ manifestly antisymmetric in a and b for having 
a complete agreement with Eq. ( |3.66|) .) Note that Eq. (|3.67|) is consistent with Eq. ( |3.64|) , 
since the analogue of the former, viz., A°-b{x) = 5^ — ^[Jcdl^'b^'^'^i^) leads to the latter on 
use of Eq. (^M )- 

Let us now return to Eq. ( p.63| ). From this equation, it is straightforward to write 
down the covariant derivative of the vierbein: 



tl. b^ " 



(3.69) 



Here we have added the usual affine connection term for the covariant derivative of a 
coordinate vector. Since the vierbein is like the square-root of the curved space metric, 
it is natural to demand that its covariant derivative be zero. Put in another way. 



(3.70) 



ensures that the covariant derivative of the metric is zero, since Eq. ( p. 45 ) holds and the 
covariant derivative obeys the Leibniz rule. The spin connections can now be determined 
in terms of the vierbeins by making use of the above condition. 

Actually, instead of using the condition ( p.70|) , it is more convenient to use the weaker 
condition 



afj, 



0. 



(3.71) 



The antisymmetrization on the L.H.S. of Eq. (|3.71|) ensures that the Christoffel symbol 
does not appear when we use Eq. ( |3.69|) in it: 



0. 



(3.72) 



Now we multiply with ec^e/e°'\ to get 



(3.73) 



Here we used the vierbein orthonormality relation Eq. ( p.46| ). Next we use Eq. ( p.72| ) to 
replace w^dae^A and ujuca^^'x in Eq- (p.73| ). After using Eq. ( p. 461 ) again, this gives 

^^f,ab = ^ ea^id^Cbf, - df^Cbu) - eb'idyCa^, - d^Cau) + eaei,^e''^{d^ecx - dxe^u) ■ (3.74) 

Finally we note that the Lagrangian given in Eq. ( p.55| ) accounts only for the gravita- 
tional interaction of the Dirac fermion. Electromagnetic interaction is introduced in the 
usual way by further demanding invariance under ?7(1)em gauge transformation: 



(/) 

G,EM 



-g 



{^i^l"ea^[d^ + ieQfA^ + -a^'u^k^ip + h.c. 



(3.75) 



where Q/ is the charge of the fermion. Since is a local Lorentz scalar and a coordinate 
vector, the new term maintains the symmetries we already had. 



3.5 Interaction vertices 

QED in curved space is governed by the Lagrangians and =5?'g em given by Eqs. 
( p.38| ) and ( |3.75| ) respectively. In this section, we expand these Lagrangians to 0{k) in 
the approximation of linearized gravity, and obtain the Feynman rules for the lowest order 
interaction of the graviton with photons and fermions. 
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— inC^j^\p{k, k ) 



Figure 3.1: Feynman rule for the coupling of the photon with the graviton. Cf^^Xpik, k') 
is given in Eq. ( ^.791 ). 



3.5.1 Photon-photon-graviton vertex 

On use of Eqs. (|3Tl| ) and ( ^lOD in Eq. ([CTI) , the 0(k) term in ^[j^ turns out to be 



with the energy-momentum tensor operator for the photon given by 



T 



(3.76) 



(3.77) 



Writing in terms of the photon field, 



f{A) 
J- Xp 



-yVXp 



{d^Ap){d"A^) - (9«A^)(a^A")l + {dxA^){d'^Ap) + {dpA^){d-A, 



-{do^Ax){d''Ap) - {dxA^){dpA'') 



(3.78) 



Let us now consider the vertex given in Fig. Each term of Jt^AAh gives two 

contributions to this vertex, depending on which of the two photon field operators in 
the term destroys the photon and which creates it. Let us take, as an example, the 
{daAf^){d'^A^) part in Eq. (|3.78|) . When A/3 destroys the photon of momentum k and A^ 
creates the photon of momentum k', the contribution is {—ika)f]i3fi{ik"^)6^; while for the 
reverse case, it is {ik'^ripy{—ik°')5^. Adding up, we find that the total contribution to the 
vertex from the — Kh^P]^rjxp{daAp){d" A^^) term of .^aah is —it^VxpVfMuk ■ k'. Carrying out 
a similar exercise for the other terms of ^aah, we finally obtain ^ 



CpyXpik, k ) 



VxpiVi^uk ■ k' - k'pk^) - r]^y{kxk'p + k'^kp) + k^{T]xpk'p + r]ppk[) 

+k' {rjxvkp + rjpykx) - k ■ k'{r]xpr]piy + VxuVpf^)- (3-79) 



29 



3.5.2 Fermion-fermion-graviton vertex 

This vertex emerges from the Lagrangian ^q^-^u given by Eq. ( |3.75| ), which involves 
the vierbein. In the weak field approximation, the vierbein can be written as 

e%{x)=5';, + K(l>\ix). (3.80) 

Here and below, we drop terms of higher order than what is necessary to determine the 
vierbein to 0(k). Note that 5^, and hence 0°^, are neither local Lorentz vectors nor 
coordinate vectors, but only fiat space tensors of rank two. Consequently, all tensors in 
an expression should be regarded only as fiat space tensors after we use the approximation 
of Eq. (|3.80| ) to perform an expansion of the expression in powers of n. The distinction 
between greek and latin indices would then be dropped, and they would be raised and 
lowered with the fiat space metric. 

The vierbein is determined from the orthonormality relation given by Eq. ( ^.45| ) or 
Eq. Thus, use of Eqs. {^^E^ and in Eq. leads to 



+ 0.^ = 2 V. (3.81) 



Eq. ( 3.81 ) does not determine the antisymmetric part of (f)^^^. This freedom corresponds 
to the fact that the curved space action for the Dirac field is invariant under local LT: a 
new vierbein, given by Eq. (|3.44|) , is as good as the old one. The ambiguity is removed 
by imposing the condition 



(3.82) 

This is called the Lorentz symmetric gauge [Q. Let us now dwell a little on this condition. 

The allowed local LT's in the weak field approximation are only those which sustain 
the closeness of the vierbein to the unit matrix. So we consider the local LT 

A\{x) = 6''h + Ke\{x). (3.83) 

Using this and Eq. (|3.8CI|) in Eq. (|3.44|) , we arrive at the transformation 

^'a, = <Pa, + Oa,. (3.84) 



Because 6a^ is antisymmetric (see Eq. (|3.65| )), we can always choose it to cancel out the 



antisymmetric part of and implement the gauge condition given by Eq. (|3.82|) . An- 
other point to note is that the gauge condition does not allow any further transformation 
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like ( |3.84|) , and so it does break the local LT invariance. Finally we remark that the 
harmonic gauge condition, which was implemented earlier using general coordinate trans- 
formation invariance, is not disturbed by the implementation of the Lorentz symmetric 
gauge condition, since Q/j^^, is a local Lorentz scalar. 

Let us now return to Eq. ( |3.80| ). Since Eqs. ( |3.81| ) and ( p.82| ) give (p^u = h^^, we have 



e\{x) = 5l + Kh\{x). (3.85) 

This implies that = rja^ + i^ha^. And use of Eq. ( p.lOj ) in = g^^Cav gives 

e/ = - Khf. (3.86) 

Now we are in a position to expand the Lagrangian -^'g em 0{k). First consider 
the terms in Eq. ( |3.75| ) involving the spin connection: 

1 



^e/w^fec(^7V''"^/^ + h.c.) 



On using the weak field expressions for the vierbeins in Eq. (|3.74D, we obtain 



Also, 



(3.87) 



(3.88) 



(3.89) 



usmg 



(3.90) 



The product of the expressions given by Eq. ( p.88|) and (|3.89|) vanish because the graviton 
field is symmetric. Thus the expression (|3.87| ) vanishes at 0{k,). 
The other terms in Eq. ( |3.75| ) not involving are 



-g 



(3.91) 



Using Eqs. ( |3.11| ) and ( p.86| ), the 0(k) term in the above expression turns out to be 



(/), 
Xp ' 



X] 



(3.92) 
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f{p) f{p') 



Figure 3.2: Feynman rule for the tree-level gravitational vertex of a fermion. Vxp{p,p') is 
given in Eq. ( |3.95| ). 



with the energy-momentum tensor operator for the fermion given by 

^ip (^) = {li^i^)bxdp + ipdxMx) + h.c.) - vxp^i^\x). (3.93) 

(see Eq. ( p.47|) for ^l/^). Here we have explicitly symmetrized the term within square 
brackets in Eq. ( 3.93|) , making use of the fact that it comes multiplied with /i^^. Writing 
out in full, 

Txp = ^ {iplxdptp + tpipdxtp - {dpij)-fx'ip - (dx^hptp) 

-Vxpilii^l'^dp^ - id,^P)ri^) - mfijij). (3.94) 
Let us now consider the vertex given in Fig. ^.2|. Each term of ^ffh gives a contribution 



to this vertex, with ip destroying the fermion and ip creating it. We thus obtain HI, ITI 

Vxpip, p') = ^[7a(p + p)p + lp{p + p')x] - \r]xp\[i> - m) + (/ - m)]. (3.95) 

3.5.3 Fermion- fermion-photon-graviton vertex 

The term involving Ap in Eq. ( |3.75| ) is 

-eQ/v^e/V'TV^M- (3-96) 

Using Eqs. (|3.11| ) and ( p.86|) , the Oi^n) term in this expression turns out to be ^ 

S£SSAK = -CKQfap^xph^yYM'', (3.97) 

with 



ai^uXp = VpuVxp - ^{VfixVup + VppVux)- 



(3.98) 



Here we have explicitly symmetrized a^^xp in A and p. The Lagrangian of Eq. (p.97|) gives 
rise to the fermion-fermion-photon-graviton vertex with the vertex factor —leK.QfapvXpl'^ ■ 
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Figure 3.3: Diagrams of 0{eK) with two fermions, one photon and one graviton in the 
external hnes. These are the tree-level diagrams for the photoproduction of gravitons off 
electrons. 



3.5.4 Important features of the vertex factors 

All the three vertex factors deduced by us are seen to be symmetric under the exchange 
of A and p. This is because they represent the coupling to hxp, which is a symmetric tensor. 
We now proceed to discuss the other important features of these vertex factors. 

The photon-photon-graviton vertex factor —inC^^xpikjk'), with C^^xpik, k') given by 
Eq. ( p.79|) , is also symmetric under the simultaneous exchange of /i and z/, and, k and k'. 
This reflects the Bose symmetry under the interchange of the two photons. The vertex 
factor also satisfies electromagnetic gauge invariance: 

k^C^,Xp{k, k') = = k'^C,,xp{k, k') . (3.99) 

It also satisfies gravitational gauge invariance: 

[k - k')\^'{k)e''{k')Cp,xp{k, k') = = {k- k'ye^{k)e''{k')Cp,xp{k, k') , (3.100) 

k — k' being the four- momentum of the graviton. In writing Eq. ( p.lOO|) , we have taken 
the polarization vectors of the photons to be real. To prove Eq. (|3.100| ), one has to use 
the on-shell conditions fc^ = = k'"^ and k ■ e{k) = = k' ■ e{k') for the photons. 

It may be noted that to prove a given kind of gauge invariance, the particles carrying 
the corresponding charge must be put on-shell. Since neither the photon nor the graviton 
carries electric charge, it was not necessary to put them on-shell in Eq. (|3.99|) . On the 



other hand, the photon has gravitational coupling, i.e. carries "gravitational charge" (this 
is in fact true for any particle, as all particles possess energy). So the photons were put 
on-shell in Eq. (|3100| ). 
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We now turn to the fermion-fermion-graviton vertex factor —iK,Vxp{p,p'), with 
V\p{p,p') given by Eq. (|3.95|) . Here again one can check gravitational gauge invariance: 



{p-p'ru{p')VxMp'HP) = = {p-p'yi^{p')Vxp{p,p')u{p) 



(3.101) 



with the electrons put on-shell for reasons stated just now. 

Finally we consider the fermion-fermion-photon-graviton vertex factor iena^yXp'-)'^ with 
o-iivXp given by Eq. (|3.98|) . The checking of gauge invariance is now less trivial. This is 
because there are three other diagrams of 0{en) with the same external lines. These four 



diagrams, given in Fig. |3.3| , represent the photoproduction of gravitons off electrons at 
tree-level (note that the third diagram from the left is just the vertex under consideration). 
When all the four diagrams are taken together, it is possible to prove electromagnetic and 
gravitational gauge invariance . 
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CHAPTER 4 



Charged leptons in thermal QED 



4.1 Dispersion relation and inertial mass 

We saw in Chapter ^ that in the real-time formulation of thermal field theory, there 
are four free-field propagators. Similarly there are four full or exact propagators with 
the Heisenberg picture fields replacing the interaction picture fields in the definitions of 
the free-field propagators. There are also four self-energies J2ab with —iT,ab (a, b = 1,2) 
denoting the IPI self-energy diagram in which the momentum enters at a 6-type vertex 
and leaves at an a- type vertex. (We are now dealing with the specific case of fermions.) 
Thus the free-field propagator, the full or exact propagator and the self-energy are all 



2x2 matrices [^, 0. The poles of the full propagator for a fermion of momentum 



occur at the poles of the function 

S'f{p) = y>-mf-^fip)]-' (4.1) 
where Sj(p) is related to the components of the self-energy matrix. Let us define 

ReS^ = i(S; + 7°4r) (4.2) 



and 

ImS^^l(S^-7°47°), (4.3) 
and likewise for Hat- Then the aforementioned relations are 3D 



ReSj(p) = ReSn(p), (4.4) 
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with rif{p) given by Eq. ( 2.51 ). 

Considering the analogy of Eq. ( |4. 1| ) to the vacuum case, we shall call the 



self-energy of the fermion. It has the general form 

T,f{j)) = ai) + bi) + c (4.6) 

at one loop in an isotropic medium with four-velocity f^. Here a, h and c are Lorentz- 
invariant functions of the two Lorentz scalars 

ui = p ■ V, (4.7) 

P= [(p-^;)2-p2]i/2_ (4 8) 
uj and P are the Lorentz-invariant energy and three-momentum respectively, which satisfy 



Eq. ([4.6|) can contain an additional term proportional to (j'^'^p^Vy in the more general case. 
However, such a term does not appear at the level of the one-loop calculations that 
we are considering in this work, and therefore we omit it. 

We shall be interested only in the dispersive part of the self-energy. This means that 
the corresponding term ipT^fip of the Lagrangian in momentum space is hermitian, which 
translates to the condition 

S^(p)=7°4(p)7°. (4.10) 

Thus, the dispersive part of the self-energy is given by ReSj. On the other hand, when 
one is also interested in the absorption of the fermion in the medium, the Lagrangian 
must no longer be hermitian, and so one has to consider the entire Sj. For the relation 
between Im S j and the damping rate of the fermion, see, for example, Refs. and ||2^ . 
We shall deal only with Re and byJ^f, we shall henceforth refer only to Re . Clearly, 
the Lorentz scalars a, b and c of Eq. ( [4.6|) are real in our case. 



Now S'j:{p) of Eq. (|]T) can be written in the form 

S'fip) = N{p)[D{p)]-' (4.11) 

with 



N{p) = (1 - a):^ -b^ + {mf + c), 
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(4.12) 



D(p) = [(1 - a)p - bv]"^ - {ruf + cf. (4.13) 

The poles of Sj{p) are given by D{p) = 0. At O(e^) (i.e., on keeping only the terms 
independent of or linear in a, b and c), this leads to 

a;2 = P"^ + mj + 2[a{uj'^ - P^) + bu + nifc] (4.14) 

= P^ + m} + ^TT[{^ + mf)J:f]. (4.15) 

Taking the square root, we obtain the pair of solutions 

1 



uj± = ±JP^ + m}±^ Tr[(|^ + my)Ej]. (4.16) 



4,/P2 + 



m 



f 



These are implicit equations for uj as function of P, since the R.H.S. contains a, b and c 
which are again functions of uj. We now write 

u^iP) = ±EfjiP) (4.17) 



Eq. ( 4.16 ) can then be used to solve for Ef{P) and Ej{P). These are the dispersion 
relations corresponding to the fermion and the antifermion respectively. In perturbation 
theory, the justification for Eq. ([4.17|) is that it reduces to the usual dispersion relations 
at the tree-level. (See also Sees. ^]3| and ^!4| on charge conjugation. For justification at 
a non-perturbative level, see the last paragraph of this section.) The inertial masses are 
defined as 

Mfj = EfjiO) = ±uj±{0). (4.18) 

Eq. ( [4.16| ) then gives 



Mfj = + ^1 Tt[{^ + mf)Ef] | ^^^^^ , (4.19) 

P=0 

using the tree-level values of uj in the O(e^) part. 

We shall perform all calculations in the rest frame of the medium, in which has the 
components given by Eq. (|2.33| ), and in that frame, we define the components of p'^ by 
writing 

p^' = {J)\p). (4.20) 

(Eq. ( [4.8|) shows that the magnitude of the rest-frame three-momentum P is indeed equal 
to the Lorentz-invariant three-momentum P. So our notations are consistent.) Then 
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a, b, c are functions of the variables and P, which we will indicate by writing them as 
a{p^, P), and similarly for the other ones, when we need to show it explicitly. In the same 
way, N{p) and D{p) of Eqs. ( ^121) and ( ^A3| ) will be written as A^(/,P) and D{p^,P). 



Since the condition uj = ±mf now reduces to p° = ±m/, we have 



Defining 



+ UTr[(l±/)S/]| (4.21) 



Sjj{p',P) = lTr[(7°±l)S,] (4.22) 
= ap° + 6±c, (4.23) 



the inertial masses to O(e^) are given by 



Mf = rrif - Sf{-mf,0). (4.24) 



Eq. ( [4.22| ) is a useful formula that allows us to extract the matter-induced corrections to 
the inertial mass directly from the one- loop expression for Sj. As we will see later, the 
wavefunction renormalization factor is determined in terms of the same quantities Sf and 

It may be noted CPT symmetry, if present, ensures that Ef{P) = Ej{P), and con- 
sequently, Mf = Mj. However, normal matter is CPT-asymmetric, and, as we shall see 
from explicit calculations, Mf no longer equals Mj. 

The antifermion dispersion relation of Eq. ( |4.17| ) and the hole dispersion relation of 
Ref. 1^ (see also Ref. |^) are based on the same concept, namely, the removal of a 



particle from the thermal background. However, Ref. |^ considers a different situation: 
high-temperature massless QED (or QCD) at zero chemical potential. In this situation, 
while one can decompose the propagator as D{p) = (tu — co'_|_)(co' — ct;_) such that lj± = ±P 
at the tree-level, there are some important differences from our case. Firstly, one gets 
two dispersion relation from uj^ itself. One relation, of course, is obtained by writing 
07+ = Ef{P). The other is obtained by writing uoj^ = —Eh{P)- The form of uoj^ is 
such that there is an Eh{P) which is positive for all P. This is the hole state, also called 
plasmino . Secondly, U- does not give dispersion relations distinct from those obtained 
from since the medium is CPT symmetric. 

38 



4.2 Spinor and wavefunction renormalization factor 



Consider the momentum-space Dirac equation 

(pV -P-l-nif- ^f{p°, P))^{P) = (4.25) 

in a medium. Requiring Eq. (|4.25|) to have non-trivial solutions yields the condition 
D{p^, P) = (see Eq. (|4.13| )). This equation, as we have seen, possesses the two solutions 
Pq = Ef{P) and Pq = —Ej{P). Now i{P) for pq = Ef{P) is the fermion spinor Us{P)- 
On the other hand, C,{P) for po = —Ej{P) is the antifermion spinor at three-momentum 
-P. 

Thus, the fermion spinor Us{P) in a medium satisfies the equation 

{^-mf-J:fip))Us{P) = (4.26) 
with pf^ = {Ef{P), P). We choose the normalization 

UliP)Us{P) = l. (4.27) 
The hermitian conjugate of Eq. ( [4.2(j| ) is 

Us{P)U-mj-J:jip))=0 



(4.28) 

where the condition ( 4.10 ) was used. We now use the general form of Eq. (^4.6|) . Multi- 
plying Eq. ( [4.26| ) from the left with Us{P)'-ffj, and Eq. ( [4.28|) from the right with 'y^Us{P), 
and adding the two resulting equations, one then arrives at the identity 

'1 - «)Pm - ^^A* 



UsiPh.UsiP) 



Us{P)Us{P). 



(4.29) 



nif + c 

The time-component of Eq. ( 4.29 ), together with the normalization adopted in Eq. ( 4.27 ), 
give 

rrif + c 



which, when fed back into Eq. ( |4.29| ), yields 



(4.30) 



il-a)Ef-b 

For future reference, we note that in particular, in the frame specified by Eq. ( p.33| ) 



(4.31) 



Us{P)Us{P) 
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p=o 



1 



(4.32) 



use Eq. (|]23]) and the first of Eqs. (|]23|), or use Eq. (|]3|) below), and 



(4.33) 



From Eq. ( |4.26| ), the exphcit form of Us{P) can easily be worked out in the rest frame 
of the medium, using any particular representation of the gamma matrices. The condition 
for Eq. ( [4. 261 ) to possess non-trivial solution is given by 



[(1 - a)Ef - by = (1 - ayP' + {ruf + c 



(4.34) 



which, of course, is nothing but the equation D{po, P) = with pq = Ej. The [/-spinors 
in the Dirac-Pauli representation are then given by 

(4.35) 



U±iP) 



(1 — a)Ef — b + ruf + c 



\ 2[il -a)Ej-b] 



X± 

'1 -a)a-P 



V (1 - a)Ef -b + mf + c^^ j 
with the prefactor chosen in accord with the normalization condition of Eq. ([4.27|) . Here 

(4.36) 



X+ = I ^ ) , X- 
The solutions in Eq. ( |4.35| ) can be used to find that 

N{Ef,P) 



Y.Us{P)Us{P) 



(4.37) 



2 [(1 - a)Ef - b] 

This relation is representation- independent. 

Eq. ( [4.37|) is useful in finding out the fermion wavefunction renormalization factor 
Zf{P). In a medium Zj is defined through the behaviour of the propagator at the one- 
particle pole po = Ef{P) |]46| : 



S'fip) 



Zf{P)Y.Us{P)U,{P) 



(4.38) 



1- particle Po — Ef 

(See App. |B|for a proof of this relation in vacuum.) Note that we have explicitly indicated 
the fact that Zf depends on P, but only on the magnitude of it, in an isotropic medium. 
This will be checked below. Near the pole considered in Eq. ( |4.38| ), Eq. ( [4.11| ) reduces to 



S'fip) 



N{Ej,P) 



(4.39) 
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The requirement that the residues of these two expressions coincide, then yields 
Zf{P) = {2\{l-a)Ef-b 



my 



(4.40) 



pO=Ef 



where we have used Eq. ( [4.37| ). Note that since a, 6, c and (hence) D depend on P only 
through P, so also does Zj. Now use 



D{po,P) = {l-2a)ipl-P^)-2bpo-mj-2mjc 



(4.41) 



to O(e^). Then, to this order. 



ZjiP) = l + 



( P'^\ da db nit dc 

V PO / OPQ OPO PO OPO 



(4.42) 



- po=Ef 



For the particular case P = in which we shall be interested, Eq. ( ^.42] ) reduces to 



(4.43) 



where 



C/ 



dSf 



(4.44) 



with Sf given by Eq. (|4.22| ) or Eq. ( [4.23|) . From now on whenever we omit the dependence 
of Zf on P in an equation, it is to be understood as the quantity evaluated at P = 0. 

In this context, we note that in vacuum one can use the result Zf~'^ = 1 — (dS j / d-p) 
just as well as Eq. ( |4.4CI| ) in order to determine the wavefunction renormalization factor. 
But since in a medium a, b and c depend on p-f , one cannot define dY^fjdj) unambiguously. 



and so Eq. ([4.40|) is the result to be used in the present case. 

We now consider the case of the antiparticles. We have already defined the antifermion 
spinor in the beginning of this section (see after Eq. ( [4.25| )). Thus, the antifermion spinor 
K(P) satisfies the equation 

(4.45) 



(|^ + m; + S;(-p))K(P)=0, 
where = (Ej(P\ P). The normalization is chosen to be 

Vl{P)V,[P) = \. 



(4.46) 
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The analogy of Eq. ( |4.29| ) in the present case is 

;i - a{-p))p^ + b{-p)v^ 



V,{Ph,V,{P) 



V,{P)V,{P). 



(4.47) 



m/ + c{—p) 

From Eqs. ( [4. 47] ) and ( [4.46D , one can obtain the expressions for V s{P)Vs{P) and 
V s{P)'^^Vs{P) just as in the case of the [/-spinors. It may be noted that since a{p) 
depends on P only through P, the quantity a{—p) essentially translates to a{p) evaluated 
at = —Ej. Similar remarks hold for h and c. 

The condition for Eq. ( [4.45|) to possess non-trivial solution is given by 



[(1 - a{-p))Ej + h{-p)Y = (1 - a{-p)YP' + {mf + c{-p)) 



(4.48) 



which is the equation D{pq,P) = with pq = —Ej. The normalized K-spinors in the 
Dirac-Pauli representation are then given by 



V±{P) 



(1 


- <-p))Ef 


+ h{-p) + mj + c{-p) 


i 


2[(l-a( 


-p))Ej + b{-p)] 


j 


(1- 


-a{-p))a-P 



(4.49) 



Here 



X- 



X- 



(4.50) 



(The choice of is motivated by consideration of charge conjugation; this will be ex- 
plained at the end of Sec. ^75| .) Eq. ( [4.49] ) leads to 

N{-Ej,-P) 



j:Vs{P)VsiP) 



Zf{P) is defined through the relation 



:i-a{-p))Ef- + bi-p) 



(4.51) 



S'fip) 



Zf{P)EVsi-P)Vsi-P) 



(4.52) 



1— antiparticle Po ~\~ E j 

(See App. ^ for a proof of this relation in vacuum.) Near the pole considered in Eq. 
(|3g), Eq. ( plD reduces to 

N{-Ej,P) 



S'fip) 



{Po + Ej) 



(4.53) 
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Eqs. (|3l|), and (|33|) give 



ZjiP) = -l2 



[l~a{p))Ej + b{p) 



'dD' 



-1 



(4.54) 



pO=-Ef 



Hence we obtain, to O(e^) 
Zf-iP) = 1 + 
At P = 0, this reduces to 

where 



/ P'^\ da db mr dc 
a+Po h^ + ^ + 



Po / dpo dpo po dpo 



(4.55) 



J po=-Ef 



1 + C/- 



C/ 



dS 



(4.56) 



(4.57) 



The tree-level spinors Ug and Vs are the limiting cases of the spinors Us and Vs when the 
effects of the medium are neglected. They satisfy the free Dirac equation in the vacuum, 
as well as the relations 



rrif 
Ef 



(4.58) 
(4.59) 



(where po = \j P"^ + ^/); with similar relations for Vg but with the substitution p^ —>■ — 
in the above equations. 

We conclude this section by pointing out why a medium is expected to impart new 
features to the wavefunction renormalization factor which are absent in vacuum. Firstly, 
had the wavefunction renormalization factor been a function of momentum in vacuum, it 
would have led to new kinds of divergence, spoiling the renormalizability of the theory. 
But there is no such problem in a medium, since the thermal corrections do not involve 
ultraviolet divergences. Secondly, CPT invariance in vacuum ensures that Zj = Zj. But 
since a medium can be CPT-asymmetric, this equality may not generally hold. 



4.3 Effect of charge conjugation invariance 



In this section, we consider the consequences of charge conjugation invariance for the 
self-energy, the dispersion relations, the wavefunction renormalization factors and the 
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antifermion spinor. C invariance is present when the background contains only photons, 
and also when it contains, in addition, charged leptons at zero chemical potential. So 
these consequences can be checked in such cases. 

We begin by writing down the effect of the charge conjugation operator C on a spinor 
field: 

%l)c{x) = C^{x)C = C-fo^ip*{x) . (4.60) 

(Note that = 1, so that = C.) C invariance of the free Dirac Lagrangian is 
ensured by the relation 

C^S'^C = -7;.^ (4.61) 

where the matrix C is unitary: 

C'^ = d. (4.62) 

Using (V'c)c = '4^1 oiie can then show that the matrix C is also antisymmetric: 

C^ = -C. (4.63) 

It will be useful to first go through the demonstration of the C invariance of the 
tree-level Dirac action 

s^f=f d^xtlj{x){iYdf, - mf)^{x) . (4.64) 



Let us define the Fourier transform of the spinor field by 

W) 



74 

^{x) = I 77^e-'P-"^(p) . (4.65) 



Putting this in Eq. ( 4.64 ), we obtain 



^/ = / J^Mp)(^ - ^fMp) ■ (4-66) 



Now from Eqs. j ^Ml and (j^Wj ) 



and so. 



Mx)= / (^^''■^C7oV(p), (4.67) 



Mx) = ^Pc\xho = -j ^e-^^-V(p)C-^ , (4.68) 
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where we used Eqs. (^4.61| ) and ( [4.62|) . Substitution in Eq. (^4.64|) with ip ipc gives 



^fC 



-^T(p)C"i(-|^-m^)C7oV(p). 



(4.69) 



Note that we have — / — in the place oi — mj of Eq. ( ^.661 ) , due to the presence of 
^ip-x ( [4.67] ) in the place of e"*^'^ in Eq. ( (4.651 ). A little rearrangement in Eq. ( |4.69| ) 

leads to 



^fc 



^{p)[C-\^ + mf)Cf^{p) 



(4.70) 



which, on using Eq. (|4.61| ), is found to equal as it should. 
Now consider the full Dirac action 

Following the derivation of Eq. ( [4.69|) , we can write down 



(4.71) 



J (27r) 

Rearranging as in the case Ajc and demanding that 



(4.72) 



/ 



we arrive at the condition for C invariance: 



fC ' 



(4.73) 



S/(P) = [C-'^s{-pn 



(4.74) 



Using the general form of Sj(p) given in Eq. ( [4.6|) , and also using Eq. ( [4.61|) , it follows 
that the Lorentz scalars a, h and c then satisfy 



a{-p) = a{p) , h{-p) = -h{p) , c{-p) = c{p) . 



(4.75) 



Substitution of these conditions in Eq. ( [4.48| ) and comparison of the resulting equation 
with Eq. ( |4.34| ) show that Ef and Ej satisfy identical equations under C invariance, i.e., 
the fermion and antifermion dispersion relations are then the same. 

Turning to the wavefunction renormalization factor, we rewrite Eq. ( [4.55| ) as 



ZAP) = 1 + 



a{-p) + (po - —] ^^^-P^ _ ^^(-P) ^ dc{-p) 
V ° Po J dpo dpo po dpo 



(4.76) 



Po=Ef 
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Substitution of the conditions of Eq. (|4.75| ) and comparison with Eq. (^4.42| ) then shows 



that Zj{P) = Zf{P) (on using Ej = Ej). So C invariance also leads to the equality of 
the fermion and antifermion wavefunction renormalization factors. 

Such symmetries between the properties of the particle and the antiparticle are, in 
general, the consequence of CPT invariance. However, the form of the self-energy as 
given by Eq. ( ^.6| ) already implies invariance of the action under parity and time reversal. 
So the additional imposition of C invariance through Eq. ( |4.75| ) is, in fact, equivalent to 
the requirement of CPT invariance. 

Finally we consider the antifermion spinor. At the tree-level, the fermion spinor is 
related to the antifermion spinor by 

t;,(P) = C7o^<(P). (4.77) 

To verify this, start with the equation (j^ — mf)us{P) = 0. Then Vs{P), as defined above, 
satisfies 

(p'^C7o%*(C7o^)-' - mf)vs{P) = . (4.78) 

This, on using the hermitian conjugation property of the gamma matrices and then Eq. 
( [4.61| ), reduces to the familiar equation + mf)vs{P) = 0. 



Let us now try to generalize this exercise in the presence of radiative corrections. That 
is, let us start with Eq. ( |4.2(j| ), and write 

K(P) = C7oX*(P). (4.79) 



Then one can find an equation for Vs{P) just as we found Eq. ([4.78|) . Demanding that this 



equation for Vs{P) be the same as Eq. ( [4.45| ) leads, on using Ef = Ej, to the condition 

S/(p) = C7o^S/(-p)7o^C-i . (4.80) 

Use Eq. (^l0|) and then Eq. (WM) in Eq. (WM) - This reduces it to Eq. (ITT^) . We 
therefore conclude that the tree-level definition ( [4.77| ) of the antifermion spinor can be 
extended to ( [4.79|) in the presence of C invariance. In the Dirac-Pauli representation. 



one can choose C = 17^7° (thereby fixing the overall phase in C which is left arbitrary by 
Eqs. ( 4.61|) , ( 4.62 ) and (^4.63|) ), while the explicit solutions for the spinors have already 



been given in Eqs. ( [4.35| ) and ( [4.49|) . One can then verify Eq. (|4.79| ), making use of the 
conditions of Eq. ([4.75|) together with Ef = Ej. 
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4.4 Relations when background is not C invariant 



Let us now consider the case when tlie background contains charged leptons at non- 
zero chemical potential. This breaks the invariance under charge conjugation, even though 
the underlying Lagrangian, being that of QED, is C invariant. 



Recall from Sec. |2.2| that the effect of the medium is introduced through the operator 
J given by Eq. ( p.4|) , which occurs in the density operator. The operator J is made up of 
H and Q, which transform under charge conjugation as 

CHC = H, CQC = -Q . (4.81) 

The relations ( [4. 811 ) can be explicitly checked at the tree-level, as follows. Considering the 
contribution from the charged leptons only, the expressions for H and Q are those given 
in Eq. ( |A.11|) . Let us now consider Eq. ( [4. 601) and substitute the expansion (|2.37|) into it. 
Using Eq. ( [4.771 ), and Us{P) = C'yo^v*{P) (following from Eq. ([4.771) and the properties 
of the matrix C), we obtain 

Ccs{P)C = ds{P), Cd,{P)C = Cs{P) , (4.82) 

so that 

CcliP)cs{P)C = dl{P)d,{P) , CdliP)dsiP)C = cliP)cs{P) . (4.83) 

(For these steps, it is convenient to use = C, which follows from = 1 and the fact 
that C is a unitary operator.) Use of Eqs. ([4.83[) in Eqs. ( |A.11| ) at once leads to Eqs. 

(EM)- 



From Eqs. ( [2. 41 ) and ( [4.81[ ), it is clear that CJC does not equal J when yU 7^ 0. Note, 
however, that J stays the same if charge conjugation is followed by a reversal of the sign 
of fi. Consequently, while the relation in Eq. ( [1.73D breaks down, the relation 

holds. Eq. ( [4. 75] ) is then modified to 

H-P)] M = ' H-P)] M = -Kp) , [c{-p)] 
Retracing the subsequent steps of Sec. |4]^, one then arrives at 

[EfiP)]„ = Ef{P), 



c(p) . 



(4.84) 



(4.85) 



(4.86) 
(4.87) 
(4.88) 
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£(p) i{p) 
(B) 



Figure 4.1: One-loop diagrams for the self-energy of a charged lepton £ in a medium. 



4.5 Calculations at O(e^) 
4.5.1 Self-energy 

The self-energy diagrams are shown in Fig. In view of Eq. ( [4.4| ), we need to 

concern ourselves only with the 1-type vertex. (Actually, the upper vertex of Fig. |4.1| B 
can still be of 2-type, but we shall see that there is no contribution in this case.) The 
propagators for the internal lines are therefore the "11" propagators of Chap. For a 
fermion, the propagator is thus given by Eqs. ( [2.4^ ), ( [2. 491 ), ( [2.5(J| ) and ( [2. 51] ), while for 
the photon, the propagator in the Feynman gauge is given by Eqs. (|2.52| ), ( p.53|) , (|2.54|) 
and (|235|) . 



The contribution for Fig. [4.1[ B vanishes because the photon tadpole is zero. There are 
several reasons for the vanishing of the photon tadpole in vacuum, but, as we shall see, 
only one of these arguments continue to hold in normal matter. First consider Lorentz 
invariance. The photon tadpole must equal a four-vector, but there being no four-vector 
available in vacuum (the photon line of the tadpole carries zero momentum), the tadpole 
has to vanish. However in a medium, we do have a four- vector, namely, the four- velocity v'^ 
of the medium, and the tadpole can be proportional to it. Next consider invariance under 
the action of the charge conjugation operator C. Using C\0) = and CA^C = —A^, it is 
clear that (0|A^|0) vanishes. But in a medium, the expectation value is with respect to the 
thermal state which need not be C-invariant (normal matter, e.g., contains electrons but 
no positrons) and in that case this argument also is invalidated. Finally one must realize 
that the photon tadpole, since it has no other external lines, represents the coupling of 
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the photon to the net current in vacuum or in the medium, as the case may be. This net 
current is zero not only in vacuum, but also in an electrically neutral medium, thereby 
ensuring the vanishing of the photon tadpole. 

As an illustration of the formal reasoning presented above, let us now explicitly eval- 
uate the photon tadpole contained in Fig. [4.ip . The tadpole is given by 

= E / ^ Tr [teQjY^Sfik)] . (4.89) 

Here the sum is over all species of fermions in the medium, the charge of each species 
being denoted hj Qf with the convention that Qe = —1. The vacuum part of Eq. ( [4.89| ), 
coming from iSpfip) of Eq. ( |2.48| ), vanishes separately for each fermion in the loop. On 
the other hand, the thermal part, coming from Srfip) of Eq. ( |2.48| ), equals 

T^ = -^AeY.Qf J -^,6{k' - m))r^^{k)k^ . (4.90) 

Let us write T'^ = v^T where T is a scalar, so that T = v^T^^. It is easy to evaluate T by 
going over to the rest frame given by Eq. (|2.33|) , and performing the /cq integration. We 
thus obtain 

_ r r rl^K 1 

(4.91) 



f 



\ff{EK)-fj{EK)] 



(27r)3' 

The distribution functions are given by Eq. ( |2.43| ). Here, as elsewhere, we write the 



components of a four-vector as k^ = {k^,K), and use the notation Ex = ^jK'^ + m'j. 
Since the number densities are given by 

^fJ = ^Jj^s ff^ ' (4-92) 

Eq. ( ^.91| ) shows that is proportional to the total charge density of the medium, and 
indeed vanishes for an electrically neutral medium such as normal matter. 

At this point, we consider the case when the upper vertex in Fig. ^.1| B is a 2-type 
vertex. Because the thermal part of the "22" propagator is the same as that of the "11" 
propagator, and the 2-type vertex differs from the 1-type vertex by a sign, the photon 
tadpole in this case is just the negative of the value given in Eq. ([4.91|) , and therefore 
again vanishes in a neutral medium. 

However, in anticipation of the calculation of the gravitational interaction of the 
fermion in Chap. ^, we mention here that in the presence of the gravitational poten- 



tial, the diagram of Fig. [4.ip is not zero by itself. This is because the condition of 
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the vanishing of the photon tadpole (which is equivalent to require that the medium be 



electrically neutral) involves other diagrams. This will be discussed in detail in Sec. |5.4.3 . 

Let us now return to our present calculation of the self-energy (in the absence of 
gravitational potential), for which we need to consider only Fig. |4.1| A. Thus the charged 
lepton self-energy is given by 

- 2S,(p) = {lef I ^ Y^Seip + k)YtD,.{k) . (4.93) 

When Eqs. ( p.48| ) and (|2.52|) are substituted into Eq. ( [4.93|) , four terms are produced. 



Since we are interested in the background induced contributions only, we disregard the 
term involving both Spe and Ap- Among the other three, the one involving both Ste and 
A J' contributes only to the absorptive part of the self-energy (see Eq. (^.5|)), and does 
not concern us here. The contributions to the dispersive part of the self-energy (see Eq. 
]4|)) arises from the remaining two terms, which can be written in the form 

Kip) = ^^ip) + K2ip) , (4.94) 



where 



Kiip) = / mv,ik) '"^^ , , (4.95) 
J [iTiy p'^ + 2k ■ p — mj 

K^ip) = -2e' f^,5{k'~mj)veik) , ^ ~ (4.96) 
J (ivr)'^ p-^ — 2k ■ p + mf 



We shall now write 



{p) = a, (p)^ + b, {p)f + ci (p) (4.97) 



and check the conditions of C invariance stated in Eq. (^4.75|) . Eq. ( [4.97|) can be inverted 



to express ai(p), hi{p) and Ci(p) in terms of Tr[|^S^]^(p)], Tr[|^S^]^(p)] and Tr[S^]^(p)]. Then 
using the expression of Eq. ( |4.95| ), we obtain 

2e2 r d^k ,,Ap-v){k-v + p-v)-{k-p + p^) 

ai[p) = r / -— TT S{k )vAk)- 



(p ■ vy — p"^ J (27r)^ p"^ + 2k ■ p — mj 

2e^ f d^k ^,j2\ Ap ■ v){k ■ p) — p^k ■ V 



^'^^^ = {p.vf-p^lj2nf^^^ p- + 2k.p-m} ' 

cM = -^e'mJ ^^6{e)vAk)^—r^ 2" (4-98) 

J (ztt)'^ p^ + 2k-p — mf 

Changing p to — p, and, simultaneously, k to —A;, in each of these expressions, one can 
readily check the properties given by Eq. ([4.75| ). 
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Let us next write 



= «2 W + Hp)^ + C2{p) . (4.99) 

We then have 

r d^k ^ ^ ■ p - {p ■ v){k ■ v) 



{p ■ vY — p"^ J {2txY p'^ — 2k-p + mj 

2e^ f d^k 2 2\ .p^k ■ V — {p ■ v){k ■ p) 



CM = 4eV/g3^(*.--mj),,W ^,_,,^^^„g . (4.100) 

To check the conditions of Eq. (|4.85| ), we first have to change /i^ to — /i^. Eq. (|2.51|) shows 
that this is equivalent to changing rii{k) to rii{—k). After making this change, one only 
has to change p to — p, and, simultaneously, k to — /c, and the conditions check out. 

4.5.2 Calculation of inertial mass 

Using Eq. (|4.22| ), and according to the decomposition given in Eq. ( [4.94|) , we write 



8t = £n + £t2. (4.101) 

where 

We can make a similar decomposition of The quantities and £^^2 ^'^^ obtained from 
£11 and £i2 by replacing by — m^. 



The inertial mass is determined by applying Eq. (|4.24|) , and, according to the decom- 
position given in Eq. ( [4.101| ), we write it as 



Mt = mi + ma+mg2-, (4.104) 



and similarly for the anti-leptons. Substituting p^ = (m^, 0) in Eq. ( |4.102| ), and using the 



fact that the terms in the integrand that are odd in k yield zero, we obtain 

— * 6^ c d^k 

ma = £ei{me,0) = — j7rT^^{k'^)ili{k) 
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This is the contribution to the inertial mass from the photons in the background, in 
agreement with the result quoted in Eq. ( |1 . 1| ) , and it is non-zero for any the charged 
lepton propagating through the medium. In a similar fashion we find 

mil = 0) = ^ , (4.106) 

and therefore the photon contribution for the anti-particle is the same as for the corre- 
sponding particle. 

In order to justify the usage of the term "inertial mass," we now deduce the O(e^) 
dispersion relation when only photons are present in the background. For this purpose, 
Eq. (|4.95|) has to be used in Eq. (f4.15|) . Since the tree- level dispersion relation = 
^2 _ p2 _ ^2 ^g^^ used in O(e^) terms, we obtain 

TT[i^ + me)Ki] = ^ , (4.107) 

after a term odd in the integration variable drops out, leading to the same integral as in 
Eq. ( |4.105|) . Consequently, Eq. ([4.15|) reads 

uj^ = P^ + m\ + —— , (4.108) 
o 

and the dispersion relations are 



E(,{P) = Ei{P) = ^P^ + (4.109) 

with 

mtT = mi + — — . (4.110) 

Thus the dispersion relations in photon background are of the same form as in vacuum, 
and reduce to Ei = Eg = itlit + -P^/Sm^ in the non-relativistic limit. From this, it 
is clear that itlit is the inertial mass. Now, from Eq. ( [4.109| ), itlit also equals -^^(O) or 
£^^(0), thereby justifying our definition of inertial mass given in Eq. ( CT ). Note that Eq. 
( [4.1101 ) is in agreement with the results given in Eqs. ([4.105 ) and ( 4.106| ). 



When the medium also contains fermions, one finds that the dispersion relations no 
longer possess the form given in Eq. ([4.109[) . This is because, unlike in Eq. ([4.107[) , 



Tr[(j^ + m£)S^2] is a function of momentum. We have chosen to call -E'^(O) and Ei{{]) the 
inertial masses even in this case, by a generalization of the usual application of the term. 
Sometimes these quantities are called the phase-space masses pO| . 
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In passing, we remark that the first equahty in Eq. ( 4.109| ) is as expected from C 
invariance, which is present in the case of a photon background. 

We now turn to the fermion background term given in Eq. ( [4.103P , and find its effect 
on the inertial mass. In a background that contains electrons but not the other charged 
leptons, the distribution functions for the muon and the tau vanish. As a result, 



(4.111) 



For the electron, we obtain 

nip J ZTT r 



kn — 2mp 



ko — rrif 

Performing the integrations over ko and the angular variables, we obtain 



(4.112) 



me2 



27r'^me Jo 



oo /<'2 r 

dK 



2E 



K 



^"-'^^Ve(i^.)+r^^)/e-(i^K; 



(4.113) 



Ek -rrie J \ Ek + 

Here Ek is given by Eq. ( [^.61 ) (with m replaced by rrie) and the distribution functions by 
Eq. ( CT ). Similarly, 



me2 



£e2i-'me,0) 

2 T{2 

dK- 



2E 



K 



J^K + me / 



Ek 



rrip 



(4.114) 



The integration over K can be performed only when the distribution functions are 
specified, and we will consider some examples in Sec. |5]^. Here we only note the general 
result of Eq. ( |4.8(j| ), obtained by charge conjugation, is borne out at P = by Eqs. ( |4.113| ) 
and ( |4.114| ): the inertial mass corrections for the particle and the antiparticle are the same 
only after we implement fie — /^e, or equivalently. 



(4.115) 



in one of them. 



4.5.3 Calculation of wavefunction renormalization factor 



We decompose 



£2 



(4.116) 
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with a similar decomposition for the anti-leptons, where, from Eq. ( 4.44 ) 





QpO 



for i = 1, 2, 



and, from Eq. ( |4.57] ), 



From Eq. ( [4.102| ), we obtain 



pf=(— ,0) 



£1 



QpO 

which imphes 



n 



QpO 



for i=l,2. 



(27r)3 K 



+ 



12m? 27r2 



°° dK 
IT 



UK) 



(4.117) 



(4.118) 



(4.119) 



(4.120) 



the photon distribution function being given by Eq. ( p. 561 ). Notice that the integral in 
Eq. ( [4.120| ) is infrared divergent. This divergence will turn out to be of importance in 
Chap. ^, since there it will cancel a similarly divergent term in the gravitational vertex 
contribution to the gravitational mass (see Eq. ( |5.59| )). 

Since the electron background terms do not contribute to the self-energy of the muon 
or the tau, it follows that 



Cm2 — Cr2 — Cm2 — Cr2 — . 

For the electron, Eq. ( 4.1U^ ) implies 



(4.121) 



dS, 



e2 



QpO 
d£e2 



pi^={me,0) 



QpO 



which yield 



p'' = (— me,0) 

Ce2 = 
Ce-2 = 



me 
nip 



me2 



(4.122) 



Trie 

ms2 



Trip 



(4.123) 



with and mg2 given in Eqs. ( [4.1 13| ) and ( |4.114| ) respectively. The result ( [4.120| ) shows 
that Zf. = Zg in a photon background, while the result ( |4.123| ) shows that transforms 
to Zg under /g ^ /g (so that becomes 171^2) as expected from considerations of charge 
conjugation. 
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CHAPTER 5 



Gravitational couplings of charged 
leptons in a medium 



5.1 Introduction 

The underlying Lagrangian for the calculations of this chapter is that of the weak 
gravitational field version of QED in curved space, which was discussed in Chapter |^. 
Thus we shall use the graviton couplings of Sec. KB, in addition to, of course, the flat 



space QED coupling. Since the coupling constant n is of mass dimension —1, the theory 
we are considering is non-renormalizable. But this will not be a problem, since our 
calculations will be confined to the lowest order in k, and the radiative corrections to be 
considered will involve an expansion in the electromagnetic coupling constant e. All our 
calculations will be valid at energies low compared to k,~^ ~ ""^pianck, as stated in Sec. 
|3.1| . (This is comparable, to some extent, to tree-level calculations using the four-Fermi 
interaction, which are valid at energies much less than m^y.) 

5.2 Gravitational mass 

The gravitational mass is a measure of the strength of the coupling of the fermion 
to the graviton. It can be determined in terms of the fermion's vertex function for the 
gravitational interaction. For this purpose, we are going to consider the S'-matrix element 
of quantum field theory for the case of scattering of a slowly moving fermion, by a weak 
and static external gravitational field (these conditions define the Newtonian limit of 
gravitational interaction). Since our aim is to calculate gravitational mass in a medium, 
we must consider this scattering to take place in the presence of a thermal background. 

So, let us first consider the relevant part 

^eff = ~Kh^'{x)fxp{x) (5.1) 
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of the interaction Lagrangian. Here Txp{x) is the total energy- momentum tensor operator, 
involving quantized fermion and photon fields, such that represents the possible 
interactions of these fields with the classical field h\p{x) to all orders in e. The S'-matrix, 
to lowest order in is then 



S=l-iKj d^x h^P{x)Txp{x). (5.2) 

Eq. ( p. 21 ) follows from the well-known S'-matrix expansion 

S = T[eM^ j d^^int)] (5.3) 

where T denotes time-ordering. The standard derivation of this formula assumes that the 
interaction does not contain derivatives of fields, so that the commutation relations of 
the interacting theory are the same as those in the free-field theory, and also J'^mt = —^int- 



In the present case, does contain derivatives through T\p (see, for example, Eqs. (|3.78 



and ( p.93| )). However the formula ( ^.3|) can be shown to be generally valid, irrespective 



of S^int containing derivatives or not, provided one uses the commutation relations of the 



free theory in calculating the S'-matrix elements |49|. This justifies our use of Eq. ([5.3| 



and also the derivation of the interaction vertices in Sees. |3.5.1j and p. 5. 2 . 



Returning to Eq. ( [5.2| ), we now consider the S-matrix element between incoming and 
outgoing fermion states: 

% ^ (/(P',.)|S|/(P,.)) (5.4) 
= -m(/(P',s)|f,,(0)|/(P,s)) U^xe-'^-^h^^ix), (5.5) 



where we have used the relation 

%{x) = e*^■^•fAp(0)e-*^•^ (5.6) 

Pp being the momentum operator, and also defined 

q = p — p' . (5.7) 

Here = Ef — E'j, each of the energies satisfying the dispersion relation of charged 
leptons in thermal QED. 

At zeroth order in e, the matrix element of Eq. ( ^.5|) equals Us{p')Vxp{p, p')us{p) with 
Vxp{p,p') given by Eq. (|3.95|) . When higher orders in e are taken into account, we write 



{f{P',s)\Txpmf{P,s)) = ^Zf{P)Zf{P') Us{P')TxMp)Us{P) (5. 
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following the LSZ reduction formula of vacuum field theory. Here rxp{p,p') is the IPI 
vertex function: it is the sum of all connected, amputated diagrams contributing to the 
fermion-fermion-graviton vertex. Though the vertex function can be defined in general 
for off-shell momenta (see Eq. ( Cll )), the arguments of Txp in Eq. ( ^.8|) are on-shell 



momenta, and satisfy the dispersion relation of the charged lepton. It may be noted that 



the replacement of the tree-level spinor Us{p) with Zf{P)Us{P) in going beyond the 
lowest order, as done in Eq. ( |5.^ ), also gives the residue of the one-particle pole of the 
dressed propagator in Eq. (|4.38| ) |pO| . 

Using Eq. ( ^.8|) in Eq. ( |5.5|) , and also taking h^'^ to be independent of time, we obtain 

Sff = -in{27i)5{Ef - E'f)Zf{P)UsiP')Txp{p,p')Us{P)h^'{Q). (5.9) 

Here we have denoted the three-vector part of g'^ by upper case letter, and also defined 
the Fourier transform by 



/(^) = / i^jiQy-"^-' (5-10) 



(27r)3' 

for any function f{x). Note that Ef = E'j implies that P = P', so that the wavefunction 
renormalization factors for the initial and the final states are the same. 

Now we relate /i'^'' to an external gravitational potential (p^^^. First consider Poisson's 
equation 

VV*(:?) =47rGp""*(f), (5.11) 

the static mass density p'^^^{x) being the source of the gravitational potential. In momen- 
tum space, this becomes 

-2gV*(Q) = /€V^*(g). (5.12) 

Next, the energy- momentum tensor corresponding to p*^^* is 

Txp = vxvpp'''' (5.13) 

where v\ is the four-velocity of the mass density. This T\p determines the desired h\p 
through the linearized Einstein field equations, which, in the harmonic gauge, read 

Uhxp = -k{Txp - \vxpT\). (5.14) 
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(To obtain Eq. ( 5.14|) , one starts with the field equations in the form 



Rxp = -SnGiTx, - ^gxpT%) , (5. 15) 

and evaluates each side to the lowest order in k. Thus, one puts Rxp correct to 0{k) 
(from Eqs. (|3.7|) and (|3.13|) ) on the L.H.S. Since the R.H.S. is already of 0(/€^), one can 
put gxp = Tjxp on this side. Then the gauge condition ( p.23|) is used.) After noting that □ 



becomes —V for the static case, we write Eq. ( 5.14 ) in momentum space, and use first 



Eq. ( ^.13|) and then Eq. ( |5.12| ) in it. In this way, we finally arrive at 



Kh^p{Q) = (2t;V - v^'WiQ) (5.16) 

It may be noted that in the frame in which the mass distribution is at rest, i.e. for 
= (1,0), Eq. ( ^.16|) immediately leads to goo = 1 + 20'^''*, which can also be deduced 
by demanding that the geodesic equation reduces to Newton's equation of motion in the 
appropriate limit Another important point is that we shall take equal to the 

four-velocity of the medium, which occurs in the thermal propagators. This means that 
we shall take the medium to be at rest relative to the mass distribution producing the 
external gravitational field (a special case being that the medium itself is the source of 
the gravitational field). Eqs. ( |5.9| ) and ( |5.16| ) constitute our end-result for the S'-matrix 



element. 

We now define the gravitational mass of the fermion. The mass density operator for 
the fermion, pf{t,x), is determined by writing another effective Lagrangian 

^'eS = -Pf{t,x)^'''\x) (5.17) 



such that it reproduces Sff as given by Eqs. ( |5.9| ) and (|5.16| ). The gravitational mass Mj- 



is then defined through the matrix element of the mass operator (obtained by integrating 
pf at any particular instant tp over all space) in the zero-momentum limit as follows: 

lim(/(P',s)| / d'xpf{tp,x)\f{P,s)) = {27r)-Y'\P-P')]p^oM'f (5-18) 

A number of remarks on this definition are in order. Firstly, the normalization of one- 
particle states adopted by us, as given by Eq. (|B.8|) , is consistent with the definition given 
in Eq. ( ^.18|) . Secondly, the mass operator is independent of x (as it involves integration 
over all space) and so carries zero momentum. Consequently, on the L.H.S. of Eq. (|5.18|) 
we have P = P'. The delta function on the R.H.S. embodies this fact (see also the 
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derivation of Eq. ( 5.21|) below.) Thirdly, we shall see that Mj- is independent of the choice 
of tp. 

The remaining task is to find an expression for the gravitational mass. Eqs. ( [5.3| ), 
and (|5.17|) give, to lowest order in interaction, 



where we used 



pf{t,x) 



-iH(t-tp)+iP-x 



and the definition of Fourier transform given in Eq. ( p.l(J| ). Again Eq. 
using Eq. ( |CT| ) with t = tp) 



M'f = lim 



{fiP',s)\pfitp,6)\fiP,s)) 



P'=p 



(5.19) 

(5.20) 

gives (on 

(5.21) 



Now P' = P implies that E'j = Ej. So the matrix element occurring on the R.H.S. of 
Eq. ( ^.21|) can be determined by using Eq. ( ^.19| ) on one hand, and Eqs. and (|5.16|) 
on the other hand. We are thus led to 

M; = {2v\^-7^^P)\imZf{P)[us{P')V^p{p,p')Us{P)\^, ^ . (5.22) 



p'^p 



Note that in the last step we have relaxed the equality P' = P to a limiting procedure, 
while retaining the weaker condition E'jr = Ef. The significance of this exercise will be 



explained in detail at the end of Sec. 5.3 



We end this section with the comment that our aim was only to get the expression for 



the gravitational mass in a medium Calculation of scattering cross-section, which 
also involves modification of phase space factors at finite temperature, does not concern 
us here. 

5.3 Operational definition at O(e^) 

We write the complete IPI vertex function at one-loop in the form 

Tap = Vxp + T'^p . 



(5.23) 



Since Vxp denotes the tree-level vertex function, F^^ is the O(e^) contribution. Then, 
using Eqs. (|5.23| ) and ( [4.43|) , the formula given by Eq. (|5.22| ) can be rewritten at O(e^) 
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in the form 



(2^V - 7]^^) X lim Us{P)Vxp{p,p)Us{P) + CfUs{P)Vxp{p,p)us{P) 



+ 



Us{P')r'{p,p')us{P) 



(5.24) 



Since C/ and T'^^ are O(e^), in any term that contains either of these factors we have 
substituted the tree-level expressions for the other quantities. Now, the terms involving 
Vxp can be evaluated immediately with the help of the identities given in Eqs. ( [4.32|) and 
( [4.33| ). Remembering that Ef{0) = Mf, we finally obtain the operational definition to 



M'j = 3Mf - 2m f + CfTUf + (2 



V V' 



^f) lim 



p'-^p 



(5.25) 



m?j in the last term. 



where we can set Ef = ^P^ 

We now turn to the derivation of an operational definition at O(e^) for the gravitational 
mass of the antilepton. We first aim to write down an expression for {f{P, s)\S\f{P', s)) 
from Eq. ( |5.9| ) by crossing. To this end, let us note the action of the field operators 
ip and ip, contained, through J^'es, in S. In the S'-matrix element of Eq. (|5.9|), ip kills 



/(P, s) and gives J Zf{P)Us{P)e-'P'' , while in (/(P, s)\S\f{P', s)), ifj creates /(P, s) and 



gives y Z j{P)Vs{P)e'^^'^ (see the tree-level Fourier expansion for ip given in Eq. ( |2.37] )). 
This does not only identify the external leg factors, but also shows that the incoming 
momentum p in the amputated vertex corresponding to the first case gets replaced in the 
second case with an outgoing momentum p, equivalent to an incoming momentum —p. 
Similarly, one can go through the action of ip. Note also that the graviton comes out with 
three-momentum P' — P in {f{P,s)\S\f{P',s)), while the fermion exchange rule gives 
an extra minus sign. We therefore have 

{f{P,s)\S\f{P',s)) = {-i){-iK)2n5{Ej-E'j) 

xZjiP) Vs{P')Txpi-P, -p')Vs{P)h^''{P' - P) . (5.26) 



Note that p and p\ of which T\p is a function of in Eq. ( p.26| ), are on-shell, and satisfy 
the antifermion dispersion relation. 
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Next, we simply interchange P and P' in Eq. ( ^.26 ) to obtain the antiparticle equation 
corresponding to Eq. (|5T 



S 



ff 



{f{P',s)\S\f{P,s)) 

{-l){-iK)27r6{Ef- - E'j)Zj{P) Vs{P)Txp{-v\ -p)VsiP')h^'{P - P')(5.27) 



This leads to an equation that is analogous to Eq. ( p.24| ), but with an extra minus sign 
in front and some obvious changes in the corresponding symbols, which in turn lead to 
the O(e^) formula 



3Mf -2mf + Cf-mf 



Vs{P)T'^,{-p',-p)vs{P') 



E'-=Ef 



(5.2J 



Using the usual relation between the free particle and antiparticle spinors given in Eq. 
( [4.77] ), and the properties of the charge conjugation matrix C given by Eqs. ( |4.61| ), ( |4.62| ) 
and ( |4.63| ), the spinor matrix element that appears in Eq. ( p.28[ ) can be rewritten in the 
form 



VsiP)T'{-p',-p)vsiP') = -u,iP')CT'l{~p',-p)C~'usiP) 



(5.29) 



Finally, using Eq. ( |C.2CI|) (which actually holds at all orders and even for off-shell mo- 
menta) in Eq. (|5.29| ), and substituting the relation in Eq. (p.28|) , we arrive at 



3Mf-2mf + Cj-ruf 
+ {2v\'' - ri^P) lim 



Us{P')T'{p,p')Us{P) 













p'^p ) 



(5.30) 



From Eqs. ( |5.25|) and (|5.3CI| ), it follows that M'j and Mj- are related by 

[M'j] „ 



(5.31) 



Since the gravitational mass of the antiparticle should equal that of the particle for C 
symmetric Lagrangian and medium, Eq. ( ^.31] ) is expected from the discussion of Sec. 
for C symmetric Lagrangian and C violating background. The consequence of Eq. ( |5.31| ) 
is that we can just obtain Mj from Mj by implementing the prescription of Eq. ( |4.115| ) 
in the O(e^) correction. 
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We take the opportunity to emphasize the following point. In the calculations that 
follow, we will find expressions for the various contributions to r^^(j9, p'), which are given 
as integrals over the propagators and thermal distribution functions. In general, such 



expressions do not have a unique limiting value as we let p' — p in an arbitrary way |53 
More importantly, some of the integrals involve ill-defined expressions if the limit is not 
taken properly. In our case, the precise order in which the various limits must be taken 
has been dictated by the physical issue at hand. Thus, since we are interested in the 
interaction of the particle with a static gravitational potential, the quantity that enters 
is r^^(p, p') evaluated for E'j = Ej. Next we set P' = P since we want the gravitational 
field to vary slowly over a macroscopic region. Finally we set P ^ to obtain the 
coupling at zero momentum, which determines the gravitational mass. This justifies the 
somewhat cumbersome notation regarding the limits in Eq. ( |5.25| ), but it is meant to 



indicate precisely what we have just explained. We will see that this prescription allowed 
us to evaluate all the integrals involved (including those that superficially seem to be 
ill-defined) in a unique and well-defined way, without having to introduce any special 
regularization technique. Our result for the case of the photon background agrees with 
the result of DHR [^, who used a special regularization technique, while our result for 
the case of the electron background agrees with the answer obtained by the use of a 
regularization technique similar to that of DHR (see Appendix [E.4D. 



5.4 Gravitational vertex 
5.4.1 Irreducible diagrams 



The irreducible one-loop diagrams for the vertex function are given in Figs. |5.1| and 



57^ . The four- vector q^, which has been defined in Eq. ( |5.7|) , gives the momentum of the 
outgoing graviton. 

We calculate only the terms that contribute to the dispersive part of the vertex func- 
tion. To determine the condition for this, let us return to Eq. ( ^.1|) , and demand that 
^eff be hermitian. This means T^^(x) = Txp{x). Next consider Eq. ( ^.8|) . Exchange of p 
and p' in this equation, followed by complex conjugation, gives 



(/(P,s)|f,,(0)|/(P',,))* = ^Zf{P)Zf{P') Us{P'hoTl{p',phoUs{P). (5.32) 
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Figure 5.1: Bubble diagrams for the one-loop gravitational vertex of charged leptons. These 
diagrams contribute in a background of photons and electrons. The braided line represents 
the graviton. 




Figure 5.2: Tadpole diagrams for the one-loop gravitational vertex of charged leptons. These 
diagrams contribute to the vertex of any charged lepton in a background of electrons and 
nucleons. 
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But since Txp is a hermitian operator, the left hand sides of Eqs. ( ^.8|) and ( |5.32| ) are equal. 
So equating the right hand sides of these equations, we are led to the desired condition 

^xp{p,p') =lo'^Xp{p',Pho (5.33) 

for the dispersive part of the vertex function. The absorptive part, on the other hand, 
corresponds to antihermitian ^eff and satisfies Txp{p,p') = —1o^\p{p' ^p)lo- It contributes 
to fermion damping, with which we are not concerned in the present work. 

When the formulas given in Eqs. (|2.48| ) and ( |2.52| ) for the propagators are substituted 
in the expressions corresponding to the diagrams, we obtain terms of different kind. One 
of them is independent of the background medium, in which we are not interested. Those 
involving two factors of the thermal part of the propagators contribute to the absorptive 
part of the vertex, while those involving three factors of the thermal part vanish because 
of the various 5-functions appearing in it. Thus, the background induced contribution to 
the dispersive part of the vertex contains the thermal part of only one of the propagators, 
and they are the only kind of term that we retain. 

We have omitted the one-particle reducible diagrams in which the graviton line comes 
out from one of the external fermion legs, because they do not contribute to F^p. The 
proper way to take them into account in the calculation of the amplitude for any given 
process, is to choose the external spinor to be the solution of the effective Dirac equation 
for the propagating fermion mode in the medium (instead of the spinor representing the 
free-particle solution of the equation in vacuum), and multiply it with the square- root of 
the wavefunction renormalization factor, as discussed in Sec. 157 



5.4.2 Bubble diagrams given in Fig. U7l 
Diagram |5.1[ A 

The amplitude of the diagram in Fig. p.l| A can be written as 
where 

k' = k-q. (5.35) 

As already explained, to determine the contribution to the dispersive part of the vertex 
function we need to retain the terms that contain the thermal part of only one of the 
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propagators. Any of them contains some combination of the form 



After some straightforward algebra, this can be written as 

Axp{ki,k2) = {ki + k2)p{hl\h + ^h>) + {^1 + ^2)\{hlph + "^hp) 
+Vxp {kl - ml){l/i2 - 2mi) + {kl - m^)(|^i - 2m^) 
+2mi{ki + k2)\{ki + k2)p . 

For the sake of convenience, we divide the total contribution into two parts 



(5.36) 



(5.37) 



V{A1) 



,{A2). 



(5.38) 



where r';^^^^ contains the distribution function of the photon and therefore contributes 
to the gravitational vertex for all charged leptons, and T'^'^'^'' contains the distribution 
function of the electrons and contributes only to the vertex for the electrons. Changing 
the integration variable from to A; + p, we obtain 



and similarly. 



AA2) 
Xp 



{p,p') 



—e 



d^k 



6{k^)v,{k) 



(27r)3 [{k + p'y - mj][{k + py - mj 



-Axpik+p,k + p')i5.39) 



d^k 



5{e - n^Mk) 



(27r)3 

Axp{k,k-q) 
[{k — qy — mj]{k — py ' [{k + — mj]{k — p') 



+ 



Axp{k + q, k) 



. (5.40) 



Diagram |5.1| B 

For this diagram 



iKVf^ip.p') = I -0^^ie^-iS,{p-k) ie^^{-iK)Cp,xp{Kk')iD^"{k)iD^^^{k'),{bAl) 



and we decompose it in analogy with Eq. ( ^.381) . The part that contains the photon 
distribution function is 



d^k 
f2^ 



^"Sp.ip - k)YCp,xp{k, k') [Ap(fc)AT(fc') + A^(fc')Ar(fc)l(5.42) 
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Making a change of the integration variable in one of the terms, this can be written as 



iP,p') 



e2 / 5{k^)7^,{k) 



{2n) 



r 



[{p' - ky - mf\{k + g)^ 



Y{i)-]/, + m,)YC^,Xp{k,k-q) 
[{P ~ " ^Tlik — qY 



(5.43) 



while 



V(B2) 
\p 



{P,P') 



d^k 

d^k 



{2n) 



YSTi{k)rC,,xp{p -k,p'- k) Af{p - k)AF{p' - k) 
6{k^-mjMk) 

(5.44) 



(p — ky{p' — ky 



gives the lepton background part. 

Diagrams |5.1| C and |5.ip 

For these two diagrams the manipulations are similar and, omitting the details, the 
results are 

d^k 



e-a^^xp / 777-7:7 S{k^)r]^{k) 



(2vr)= 



{k + p'y — mj 



{k + pY — m| 



(5.45) 



and 



p/(C2+D2) 
^ Ap 



(p,p') 



e a 



puXp 



d^k 

(2^ 



6{e - mj)r],{k) 



{k-p'Y {k-pf 



(5.46) 



5.4.3 Tadpole diagrams given in Fig. |5.2 
The question of the photon tadpole 



The calculation of the diagrams of Fig. |5.2| leads to an immediate problem. While the 



contributions from the diagrams ^.2P , |5.2| C and |5.2| D appear to vanish due to the charge 
neutrality of the medium (see Sec. [4.5|) , the contribution from the diagram ^.2| A appears 
to diverge for g = (due to the divergence of the photon propagator in this limit). The 
following discussion is devoted to the setting up of a framework of calculation which leads 
to the disappearance of the aforementioned divergence. 
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Figure 5.3: The one-loop diagrams that contribute to order k, to the photon tadpole in a medium, in 
the presence of a static and homogeneous {q = 0) gravitational potential. The fermion loop involves 
a sum over all the species of fermions present in the medium. 



Let us recall from the discussion of Sec. [4.5| that the electrical neutrality of the medium 
requires the vanishing of the photon tadpole. The one-loop diagrams that contribute to 
the photon tadpole in the presence of a static and homogeneous gravitational potential 
are shown in Fig. |5.3| , where the graviton line represents represents the g = background 
field. In the absence of the background field, only the diagram OA contributes to the 
photon tadpole. In that case, the requirement that the tadpole vanishes yields the familiar 
condition (see Eq. ( [4.9 1|) ) 



f 



■[fj{EK)-ff{EK)] 



0. 



(5.47) 



and the quantity is identified with the total charge of the medium. However, in 

the presence of the background field, and to the order that we are calculating, we have to 
take into account the contributions of the diagrams ^.3P and |5.3| C to the photon tadpole 
or, equivalently, to the total charge of the system. If we denote them by and 
respectively, the condition for the vanishing of the photon tadpole is 



Q 



(5.3A) 



+ Q») + Q 



(5.3C) 







(5.48) 



instead of Eq. ( |5.47] ). Physically, this means that the number density of the particles are 
not determined by their free distribution functions. The particle distributions rearrange 
themselves in a way that depends on the background gravitational field. 
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Consider now the problem of finding out the dispersion relation of the fermion in 
an external gravitational field. The dispersion relation is determined by the poles of 
the function given in Eq. ( [4.1| ), with now denoting the fermion self-energy in an 

external gravitational field. Thus, the one-loop diagrams contributing to at 0(/t) 



are the diagrams of Fig. |4.1| , Fig. f).l\ and Fig. [5.2| . Let us confine our attention to 
the tadpole diagrams, i.e., the diagram of Fig. |4.1| B and the diagrams of Fig. |5.2| , and 
consider the implications of the charge neutrality condition of Eq. ( 5.48| ) for the sum of 
the contributions arising from these diagrams. 



Firstly, the unadorned tadpole of Fig. |5.3| A is now itself of order k because of the 
condition (|5.48|) . Since the diagrams |5.2| C and |5.2| D contain an explicit factor of k apart 
from the unadorned tadpole, their contribution is actually of order and therefore we 
can neglect them. Secondly, the diagram [4.ip cancels the g-independent contributions 
from the diagrams ^.2[ A and p.2p . Since the loop in diagram p.2p in independent of q, 
this diagram is totally cancelled. 

In summary, the only contribution from the various tadpole diagrams arises from the 
g-dependent part of the tadpole subdiagram of Fig. |5.2| A (by subdiagram, we mean that 
the lower part, beginning with the photon propagator, is amputated). It seems reasonable 
to assume that this holds true also when we consider the gravitational vertex function and 
its contribution to the gravitational mass. It then turns out that the g-dependent part of 
the tadpole subdiagram of Fig. |5.2| A, when multiplied with 2v^v'' — rj^^, and evaluated at 
go = and Q 0, is proportional to v^Q^ (see Eqs. (|^ ), (^^) (|KiT|) and (^^)). 
When multiplied with the photon propagator at go = 0, it gives zero for the 5(Q^) part 
of the propagator, and cancels the l/Q"^ in the other part. This latter contribution will 
be labeled by the letter 'X' in order not to confuse it with the contributions of Fig. ^.1| A. 

The non-vanishing contribution 

We denote the vertex contribution coming from Fig. ^.2|A by 



ri?(p,p') 



(5.49) 



where X\pa{q) is the photon-graviton mixing diagram with external momentum g 



Tr 



(5.50) 
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Then, taking the above discussion above into account, the quantity which will appear in 
the expression for the gravitational mass is given by 



{p,p') 



(5.51) 



(As mentioned in Sec. O, the sum in Eq. ( 5.5CI| ) is over all species of fermions in the 
medium, the charge of each species being denoted by Q/ with the convention that Qe = 
— 1.) The medium-dependent contribution to Xxpa{(l) can be written as 



Axpa{k, k-q) ^ Axpa{k + q, k) 



q^ — 2k ■ q 



q'^ + 2k ■ q 



(5.52) 



where, for arbitrary 4- momenta ki and ^2, we have defined 



Axpa{ki,k2) 



Tr 



Vxp{ki, k2){h + "^/)7a(l^2 + ruf) 



{2kixkip + kixk2p + k2xkip)k2a + {m) - ki ■ k2){r]xakip + r]pak 
-2r]xp{kf - m))k2a + ki ^ k2 



(5.53) 



Putting k'^ = m'j, we obtain 



Axpa{k, k-q) = [Skxkp - A{kxqp + kpqx) + 2qxqp]ka - [4:kxkp - {kxqp + kpqx)]qa 

+k ■ q[vxa{2k - q)p + 7]p^{2k - q)x] - 27]xp{q^ - 2k ■ q)k^ . (5.54) 

Since Axpa{ki, ^2) = ^Apa(^25 ^i) by definition, Axpa{k + q, k) is obtained by changing the 
sign of q in this expression. 

5.4.4 Additional remarks 

This is an appropriate point to consider the cases when 2-type vertices can occur in 
the diagrams contributing to the gravitational vertex. Since the vertices connected to the 
external legs can only be of 1-type, a 2-type vertex cannot occur in any of the diagrams 



of Fig. |5.1| or in the diagram of Fig. |5.2p . However, in Figs. ^.2| A, |5.2| C and |5.2| D, the 
vertex where the photon line connects to the fermion loop can be of 2-type. 

Now, due to the reason given in Sec. |4.5.1| (see after Eq. ([4.92|) ), the unadorned tadpole 
with a 2-type vertex is just the negative of the unadorned tadpole with a 1-type vertex, 
and hence of 0{k) because of the charge neutrality condition. So the diagrams of Figs. 



5^ C and |5.2p , even with a 2-type vertex, are of 0(fi; ), and can be neglected. 
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In Fig. |5.2| A, the occurrence of the 2-type vertex makes the photon propagator a "12" 
propagator. So the photon propagator is proportional to which becomes for 

= 0. However, the value of the graviton momentum is taken to be Q — > 0, and not 
Q = (having Q = in addition to = would mean a potential constant in space and 
time, resulting in zero gravitational field). Consequently, S{Q^) in the "12" propagator 
causes this contribution to vanish. Therefore, 2-type vertices do not play any role in the 
present calculation. 

Next, we make some comments on the expressions deduced in this section. 

Firstly, one can easily check that the expressions given in Eqs. ( |5.39D , (|5.4CI| ), ( |5.43|) , 



( ^.44| ), (|5.45| ), (|5.46|) , and (|5.49|) (coupled with (|5.52|) ), all satisfy the dispersive condition 
of Eq. ( |5.33|) . (Similarly, it can be checked that the terms involving the thermal parts of 



two propagators from T'1^\ T'lf\ ^\p^^^ and T'j^^^ are all absorptive.) 

Secondly, it can be verified that the different contributions to the dispersive part of 
F';^^ also satisfy the charge conjugation relation given by Eq. ( |C.2(]| ). 

Thirdly, the complete one-loop vertex function satisfies the transversality condition, 
which is implied by the gravitational gauge invariance of the theory. This is shown in 
Appendix 0. 

5.5 Calculation of the gravitational mass 

As seen in Eq. ( [5. 251 ), there are three types of O(e^) correction to the gravitational 



mass. One of them is proportional to the inertial mass that was calculated in Sec. |4.5.2| , 
and another one involves the wave function renormalization factor derived in Sec. 4.5.3 . 
In this section we find the contributions from the irreducible one-loop vertex diagrams. 
Since the expressions for those already have an explicit factor of outside the integral, to 
evaluate them we can use the tree-level values for the dispersion relation and the spinors 
associated with the external lepton. 



5.5.1 Terms with the photon distribution from Fig. 5.1 



We first evaluate those terms obtained in Sec. ^]4| that contain the photon distribution 
function. In fact, if the temperature of the ambient medium is low (T <^ rrie) and the 
chemical potential of the background electrons is zero, these are the only terms that 
contribute and they are precisely the ones that were calculated in Ref . [EO] . Since we have 
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performed the calculations in a different way, using 1-particle irreducible diagrams only, 
the following results serve as a good checkpoint between the earlier calculations of Ref. 
| 20|] and ours. 



Contribution (Al) 



From the formula for the gravitational mass given in Eq. ( p.25| ), it follows that we 
need to calculate the vertex only for p = p', in which case 



'^'^^ ^ J (27r)3 4(A;-p)2 

From Eq. ( ^.37] ) it follows that, for any 4- vector y'^, 



which leads to 



(5.55) 



(5.56) 



J (27r)3 (k-py 

n , ^ (k-pY 

+me{kxkp + pxPp) H rjxp 

nil 



k ■ p 



k\Pp 



mi 

UsiP), (5.57) 



where we have used Eq. ( 4.58 ) and omitted the terms odd in k, which integrate to zero. 
Using the notation 

'(^1) = {2v'vP - r^'n [usiP)T';:;'\p,p)usiP)]^^^^^^.^ , (5.58) 

we obtain 



m 



(Al) 



UK) 



1 



mi 



miK 



e^mi dK 

Jo -Y^'^^^ 



(5.59) 



12mi 27r2 

The remaining integral is infrared divergent, but its contribution to the gravitational mass 
is canceled by a similar term that arises from the wavefunction renormalization, as we 
show below. 

Contribution (Bl) 

This term has to be treated carefully because the denominators in the integrand of 
Eq. ( p.43| ) vanish for g = 0. Indeed, by using the formulas 

- I = lim ^ , 
5{x) = lim - ^ 



(5.60) 
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(with fPr denoting the principal value), one can show that |5|] 
^ ( 1 \^/,o Ox \ d 



^2 _ 



(5.61) 



2 dm"^ 

which is ill-defined even for m ^ 0, whereas for m = 0, one can interpret the R.H.S. of 
Eq. ( p.61j ) only as a limiting form with m — > 0. However, a careful evaluation of the term 
given in Eq. (|5.43| ), following the procedure indicated in Eq. ( ^.25| ), shows that the limit 
exists. Denoting 



m 



(Bl) 



= {2v\f - ri^P) lim 



P'^P 



the result is 



m, 



(5.62) 



(5.63) 



The details of the derivation of this result are given in Appendix [E.l| . This result can 
also be obtained by means of a special regularization at g = (following Eq. ( ^.61D ), as 
shown in Appendix [E.4 . 



Contributions (Cl+Dl) 

We can proceed as in the evaluation of m|^j^-j above. Thus, from Eq. ( |5.45| ), 



Us{P)Tf;^'''\p.p)us{P) 



pf = (m.^,0) 



2 [ j£k_ 6{k^)7]^{k) 

6 ^fMuXp j /r)^\ 3 



X 



(27r)3 k-p 

Us{P)Yh''^^s{P) 



using the fact that a^i^xp is symmetric in the indices /x, u. Then using 
it follows that 



(5.64) 



(5.65) 



m 



(^1+^1) ^ {2v'v''-r^'n[ns{P)^'iT^'''\p.P)Us{P) 
= 0. 



(5.66) 



5.5.2 Terms with the electron distribution from Fig. |5.1 



These terms contribute only to the gravitational vertex involving electrons and 
positrons. The integration over fco and the angular variables can be done exactly. The 
remaining integral can be evaluated analytically only for special cases of the distribution 
functions, some of which we consider afterwards. 
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Contribution (A2) 



As can be seen from Eq. ( [5 .401 ), the denominators of the integrand of this term vanish 
as g 0. Consequently, the prescription indicated in Eq. ( 5.25| ) has to be followed 
carefully in this case. As we show in detail in Appendix |E.2| , defining 



m 



P'^P 



(5.67) 



the final result for this term is 



m 



iA2) 



d^K \ 2El -ml [Ek- 2m, df, ^Ek + 2m, df. 



m, J {2t:Y2Ek \ Ek \ Ek — rrie OEk Ek + rn, OEk ^ 
t - Elm, - 5El7 



2E% - Elm, - hElml + 2EKml - 2m^ 



m,El{EK - m,) 
2El + Elm, - 5Elml - 2EKml - 2mt 



fe 



(5.68) 



m,El{EK + m,) 

where Ek is defined in Eq. ( p.6|) . This result is also derived in Appendix [E.4| by means 
of a special regularization at g = 0. 

Contribution (B2) 



From Eq. ( ^.44| ) it is seen that the integrand is not singular in the limit g ^ 0. 
Therefore we can evaluate directly 



V(B2) 



(p,p) 



2/ 5{k' - ml)r^,{k)Y{l/^ + m,)r ^"'^'^[^ _ ^jf ^\ 5.69) 



and the contribution to the gravitational mass is given by 



m. 



pM = (me,0) 



(5.70) 



In the expression for C^y\p, any term having a factor of {p — k)^ or (p — fc)^ does not 
contribute to the integral. This is because, within the spinors, we can write 



+ me)7np - = 7^(1^ + m,){m, - ft) = Yimt - k') , 



(5.71) 



which vanishes because of the 5-function. The argument is similar for (p — k)y. Thus, 



(2t;V - r7"0^.(^)7"(l^ + m,)^''C^,xp{p ~k,p- k)u,{P) 

= Us{P){s{f - 2m,){p -v-k-vf + A{m, - 2k ■ vi;){p - kf]us{P) 



(5.72) 
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ignoring all terms which have a factor of k"^ — m^. Using Eqs. (|4.58| ) and ([4.59 ), we then 
obtain 

,2 



m 



{B2) 



- — / 7^ ^ ^ - ^e)^e{k) ko + 

ml J {2'kY \ ko — nie 



.(5.73) 



Contributions (C2+D2) 

Similarly, for this term we can evaluate directly 

with 



pM = (me,0) 



■^l(C2+D2) 



By straight forward algebra 



{k - 



and using Eqs. ( [4.58|) and (|4.59| ), 



m 



(C2+D2) 



6e' 



d^k 



Si^-mDri^ik) 



1 



kn — nip 



27r2 J Ek 



UEk) UEk) 



Ek — TTie Ek + m, 



5.5.3 Terms from Fig. [5T2 

The contribution to the gravitational mass due to this term is 



m, 



= {2v\P - ri^P) lim 
where, from Eq. ( p.51| ). 



q'J=Q 



qO=0 



Using the expression for A\pa from Eq. ( p.54| ) we obtain 



(2t;V-r/^0^Apa(g) 



d^k 

(2^ 



5{k'^ - 171^)7} f{k)koVa 



(5.74) 



(5.75) 



(5.76) 



(5.77) 



1 



1 



X 



,2ir-g-Q2 2K-Q + Q'^, 
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(5.78) 



(5.79) 



(2^0^ - mj - ig^) 
(5.80) 



where we have omitted the terms that vanish by symmetric integration over K, as well 
as all those terms that are independent of q, because they drop out of Eq. ( |5.79| ), and 
in addition all the terms that are proportional to g^, because in Eq. ( |5.78| ) they yield a 
factor of ^ which vanishes between spinors. Performing the integration over kg, 



(2t;V-r7^0^Apa(g) 



4-«i:s//||(//-//)(2s|-™)-5e^) 



qO=0 



{2tt 



X 



(5.81) 



,2K-Q-Q^ 2K-Q + Q^, 

For the term that contains an explicit factor of in the numerator we use the angular 
integration formula of Eq. ( |E.9| ) , which yields 



{2v^vP-r]^P)Xxpa{q) 



qO=0 



f 



x{2E 



2tt'^ 



m, 



. (5.82) 



,2K-Q-Q^ 2ir-g + gV. 

The evaluation of the rest of the integral is presented in Appendix [E.3| . Substituting the 
results into Eq. ( |5.79|) , the contribution of this diagram to the gravitational mass is found 
to be given by 



m 



5.5.4 Summary 



V 67r2 J 



ff-ff - 



2El~m) d 



2E 



K 



dE 



K 



{ff - ff) 



(5.83) 



Starting from Eq. ( ^.25|) , the total contribution to the gravitational mass of charged 
leptons can be written in the form 



= m£ + + + "^(x) ) 



(5.84) 



where is the contribution from Eq. ( 5.83 ), which is the same for all charged leptons, 
m'^i represents the terms that contain the photon distribution function, and contains 
the terms that depend on the electron distribution function. They are given as follows. 

Substituting into Eq. (|5.25|) the results given in Eqs. ( |5.59|) , ( p.63| ) and ( |5.66|) , and 
using the expression for the wave-function normalization and the correction to the inertial 
mass given in Eqs. ( |4.120| ) and ( [4.105| ), we find 



m 



12m £ 
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(5.85) 



in agreement with the DHR result [2^], quoted in Eq. ( |1.2|) . Notice that the infrared 
divergence contained in the rn'^^^^ term cancels with a similar one that arises from the 
wave function renormalization correction Qi. 

The terms from the diagrams in Fig. |5.1| that involve the fermion distribution function 



contribute only to the gravitational mass of the electron, and therefore 

<2 = <2 = 0. (5.86) 

The individual contributions of this type to the electron gravitational mass appear in Eqs. 
( p.68| ), ( |5.73| ) and ( [5.77| ). Substituting those results into Eq. ( |5.25| ), and using the results 



for the inertial mass and the wave-function normalization factor, given in Eqs. ( |4.113| ) 
and ( |4.123|) respectively, we obtain 



m'o = — / dK —^{\- + -^-- \ feiEx) 

2El~ml f EK-2m, df, + 2m, df, \\ 

2Ek XEK-nie dEx Ek + rrie OEk))' 

The corresponding formulas for the antileptons are obtained by making the substitu- 
tion given in Eq. (^4.115|) , as explained after Eq. ( ^.31|) . Thus, 



M| = mi + m'^^ + ~ ""^x i (5.88) 

where 

<2 = <2 = 0, (5.89) 

while the result for is obtained from Eq. ( |5.87| ) by making the substitution <-> /g. 

5.6 Results for particular cases 

In contrast with the mn and m^^ terms, which depend on the photon momentum dis- 
tribution, me2, Tn^2i ^'e2i ''^'e2 ^'{x) depend on the fermion distribution functions and 
cannot be evaluated exactly in the general case. Therefore, for illustration, we consider 
in detail their calculation for the specific situation in which the background is composed 
of non-relativistic protons and electrons. In this case we can set (for f = e,p) 

ff-{E) ^ (5.90) 
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and 

We consider in detail two cases separately, according to whether the electron gas is classical 
or degenerate. 

5.6.1 Classical electron gas and classical proton gas 

In this situation we can set 

fj{E) = e-'^(^-^/) (5.92) 
for both f = e,p. This implies the relation 

^ = -f^ff^ (5-93) 

as well as the integration formula 

/dA-A-/,^2.V.r(. + i)(^)"',V. (5,94) 

where Uf is the number density, given by 



Let us consider and me2, given in Eqs. ( |4.113| ) and (|4.114| ), respectively. Setting 
/e = and using Eq. ( |5.91 ) to expand the co-efficients of fe in the integrands in powers 



of K, the remaining integrals are evaluated by means of Eq. (|5.94| ) to yield 



6^n 

me-2 = ^ + 0(eVTVm^). (5.96) 



(Note that the subleading terms are small because they are suppressed by powers of T/rrie, 
and T < mg for the electrons to be non-relativistic.) Similarly, from Eq. (|5.83|) we obtain 
for this case 

,2 



If E %^ + 0(eVT/mD 



W 24r^— p rrif 



24meT 



(5.97) 
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where we have used the charge- neutrahty condition which, neglecting terms 0(k), is 
simply J2f Qf^f = 0. Applying the same procedure in Eq. ( |5.87| ), the leading contribution, 
in powers of T/nie, comes from the dfe/dEx term in that equation, and leads to 



K2 



2T2 



(5.98) 



By the substitution indicated in Eq. ( ^.1151 ), the corresponding result for the positron is 



m 



e2 



SrripT 



+ 0{e^ne/m; 



(5.99) 



Therefore, using Eqs. ( 4.104 ), ( |5.84| ) and ( p. 881 ), the inertial and gravitational masses for 
charged leptons i other than the electron are obtained as 



Mg = Mi 
M'- 



m + 



mi - 



± 



(5.100) 



12mi 24meT 

where the upper sign corresponds to the leptons and the lower one to the anti-leptons. 
The corresponding formulas for the electron are 



Me 
M' 



me + 



12me 
12m, 



+ 



2meT 



2T2 ' 



(5.101) 



and for the positron they are 



ML 



"^e + — \- 



12me 



imt 



2e^np 



mp — 



(5.102) 



12me "im^T 

We now estimate how large these corrections could be for the electron. Those due to 



the photon background were estimated by DHR [2C] and were found to be very small. 
This is because the fractional changes in the inertial and gravitational masses in that case 
are 0{e'^T'^ /m^), and T < m^. Therefore, neglecting that contribution, the fractional 
changes are given by 



Me — me 



mp 



M' 



me 



mp 



2mlT 
2meT2 



(5.103) 
(5.104) 
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Although it may seem that the effects are more noticeable as the temperature decreases, 
they are bounded by the conditions that the electron gas is non-degenerate and non- 
interacting, which require that 

^2/3 

T>-^ (5.105) 
vrie 

and 



r 



(5.106) 



since rav ~ |5^. Using the fact that Up = Ue, it follows that the corresponding con- 



ditions for the proton gas do not imply further restrictions, because they are automatically 
satisfied whenever Eqs. ( [5. 1051 ) and ( |5.106| ) hold. 

By writing the right-hand side of Eq. ( ^.103| ) in the alternative forms 

2 



it is seen that 



Similarly, writing 



it follows that 
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1/3^ 



T 
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n: 



2/3 
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2m, 



e nip 



n. 



1/3 



Mp 



nip 



nip 



< 



if n^^ < e^nrip 
Tl2nip, if n^J"^ > e^nip 



'^2/3\2 



n 



1/3 



' e^nie\ 



n 



1/3 



- nip 



nip 



< 



(5.107) 



(5.10J 



(5.109) 



(5.110) 



in either case. Therefore, while the fractional correction to the electron's inertial mass is 
likely to be small in most situations with the conditions that we are presently considering, 
the fractional change in the gravitational mass could be substantial. For example, if we use 
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the temperature and density at the solar core, i.e., T = 1.57 x 10 K, rig = 9.5 x 10 cm 



we obtain 



Mp — nip 



nip 



ML - nip 



nip 



9.8 X 10" 



3.5 X 10" 



(5.111) 



which shows that the correction to the gravitational mass of the electron can at least be 
appreciable in realistic physical situations. 

5.6.2 Degenerate electron gas and classical proton gas 

For a degenerate electron gas. 



„2/3 



nip nip 



where Kp is the Fermi momentum of the electron gas. We assume 



Kp < nip 



(5.112) 



(5.113) 



so that the electrons are non-relativistic. We also assume the electron gas to be non- 
interacting. This requires that the average kinetic energy ~ K\/nip of an electron is 
larger than the average Coulomb interaction energy ~ e^n^J'^ ~ e^Kp. (Unlike the case 
of the classical gas, T contributes negligibly to the average kinetic energy in the present 
case; see Eq. ( p.ll2| ).) This implies that |^ 

Kp>e'^mp. (5.114) 

We now show that, under these conditions, the protons can be treated as a non-interacting. 



non-relativistic, classical gas The only additional condition we have to impose is the 
condition for non-interaction for a classical proton gas: 



(5.115) 



(usmg Up 



Hp 



. Since Eqs. ( |5.112|) and ( |5.113| ) only imply that T < Kp, it is still possible 



to satisfy Eq. ( p.ll5| ). Eqs. ( ^.112|) and (|5.113|) also imply that T < nip. So the non- 
relativistic condition for the classical proton gas is automatically satisfied. Finally, since 
Kp < nip <^ nip, Eq. ( |5.115| ) gives T ^ e'^K^p/nip, so that the non-degeneracy condition 



Kl 



n: 



2/3 



nip nip 



(5.116) 
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is satisfied as well. 

Therefore, Eq. ( p. 921) applies to the proton, while for the electron 



/e = e{KF - K) 



(5.117) 



with 



Kp = (Svr^n, 



(5.118) 



which in turn imply the relation 



dK 



-5{Kf - K) . 



(5.119) 



We repeat the calculation of the quantities me2, Tne2, ^'e2i ^'e2 ^"^^ ^'x ^r this case, 
neglecting the terms that are a factor ~ 0{Kp/ml) smaller than the ones that we retain. 
From Eqs. ( |4.113|) and ( [4.114| ), setting /g = and using Eq. (|5.91| ), we obtain 



me2 
me2 



27r2 



(5.120) 



From Eq. ( ^M ) 



m 



dK 



^^~[^+2K dK 



7e^ Up 
24Tm^ 



(5.121) 



where we have borrowed the result for the proton contribution from Eq. (|5.97|) , while in 
the electron term we have expressed Ek in terms of K and used 

d Ek d 



dEK K dK 

for any function of Ek- Using Eqs. ( p. 1171 ) and ( |5.119| ) this finally yields 



(5.122) 



m. 



(5.123) 



" 12^ • 

Here we have neglected the proton contribution because it is ~ e^K^/[Tmp) < e^Kp 



from Eq. ( p.ll6| ). In a similar fashion, from Eq. ( ^.87|) 



m. 



e2 



2tt^Kp 



(5.124) 
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and by the substitution indicated in Eq. ( [4.115| ), the corresponding result for the positron 
is 

<, = -5g£. (5.126) 

Thus, substituting these resuhs into Eqs. ( |4.104| ), ( p.84| ) and (|5.88| ), we obtain the 
following expressions for the inertial and gravitational masses, retaining only the leading 
terms in powers of Kp/nie- For the charged leptons £ other than the electron. 



Ml = Me = 1711 + 



12me 



with the upper sign corresponding to the leptons and the lower one to the anti-leptons, 
while for the electron 



Me = me + 



12m, 27r2 



and for the positron 



- --12;;^ +12^' ^'-'"'^ 



Mg = me + h 



12me ' STr'^ml ' 



2rp2 2 2 

M' = m,-— (5.128) 

It is interesting to note that Eqs. ( ^.112| ) and ( 5.113|) imply that the photon contribu- 
tions in Eqs. ( |5.126| )-( ^.128| ) are much smaller than the contribution due to the electron 
background in each case. In fact, using Eq. ( ^.114| ), we see that the fractional corrections 
to the gravitational mass can be as large as about 7/127r^ for the electron and l/127r^ for 
the positron and the other leptons. 
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CHAPTER 6 



Some aspects of high-temperature 

QED 



6.1 Introduction 

So far in this work, we have considered low temperatures and densities. Such thermal 
backgrounds ensure the validity of a straightforward expansion in e^, something we have 
always assumed. A negative illustration of this point is provided by Eqs. ( |4.104| ) and 
( [4.105| ) for the inertial mass of the charged lepton in photon background. It is clear from 
these equations that if T exceeds rrii by a large amount, the change in inertial mass can 
no longer be treated perturbation. 

In this chapter, we turn to an entirely different kind of thermal background, namely, a 
background at high temperature. More specifically, we will now consider high-temperature 
QED with exact chiral invariance and at zero chemical potential. In this case, the medium 
effects are not, in general, of O(e^) (see, for example. Sec. |6.5.1| ). Exact chiral invariance 
means that the full Dirac action of Eq. ( |4.71| ) is invariant under ip —>■ e^^^^ip for any real 
constant 9. So, the electron is massless in vacuum, and the self-energy cannot have the 
constant term on the R.H.S. of Eq. ( [4.6|) . (Chiral symmetry also rules out the (J^'^p^Vu 
term in the self-energy, which, anyway, is absent at one-loop, beyond which we will not 
be interested in.) Thus, we now have the general form 

^{k)=a]/^ + bi) (6.1) 

for the self-energy of the electron. We shall use to denote the electron four-momentum 
in this chapter [^]. The scalars a and h are then functions of 

uj = k-v, K=[{k.vf-k''Yl^ (6.2) 
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which satisfy 



(6.3) 



The pole in the full propagator is given by 

f{uj,K) = {uj- K){l-a)-b = Q. (6.4) 

(This is nothing but Eq. (|]3|) at zero mass. Since we now take the chemical potential to 
be zero, C symmetry is unbroken, and the particle and antiparticle dispersion relations 
are the same. So, instead of Ef and Ej, we just use u to denote the energy.) To leading 
order in T and at one-loop, it was shown by Weldon in his seminal work 0] that 



l-^ln(^ 
2K \uj-K 



8K 



---4-11" 



Substitution of Eq. ( |6.5| ) into Eq. 



then gives the electron dispersion relation 



u-K 



8K 



cj\ (uJ + K 



(6.5) 



(6.6) 



An important point is that if we write uj—K = 6/(1— a) ^ b from Eq. ( |6.4| ) by neglecting 
and ab, and then substitute Eq. ( |6.5| ), we do not get the correct answer given in Eq. ( [6.61 ). 
This is an example of what we meant by the invalidity of a straightforward expansion in 
at the beginning of this section. 

In the same work, Weldon demonstrated the gauge independence of this leading order 
dispersion relation, by showing that gauge dependence appears in the self-energy only at 
subleading order in temperature. However, one expects that even this subleading order 
gauge dependence should not show up in the dispersion relation. We are therefore going to 
investigate the gauge independence of the dispersion relation after taking the subleading 
temperature dependence into account [B^. 



6.2 One-loop electron self-energy in a general linear covariant 
gauge 



Re S is defined through the equations 

a{uj,K) = a{i{uj, K) + iaj{uj, K) , 
b{uj,K) = bR{uj,K) + ibi{u,K), 
ReS = aR^ + bR^). 
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(6.7) 
(6.8) 
(6.9) 



To obtain the dispersion relation, we need to consider only Re S and, consequently, only 
the "11" propagators, at one-loop As in Chapter]^, we will denote ReS by just S. 
Because of the additional features of masslessness of the fermion, zero chemical potential, 
and inclusion of gauge dependence in the photon propagator, it will be convenient for 
us to rewrite the "11" propagators for the purpose of the calculations in this chapter. 
Thus, we are going to use the free-particle massless fermion propagator at zero chemical 
potential given by 

1 



Sip) 
fpip) 



jp- -|- %e 

|Po|/T _^ 



and the free photon propagator in a general covariant gauge given by |60| 



^rG(p)+^r(p)' 

^ -27rz5(p2)/^(p) 



—1]' 



(6.10) 
(6.11) 

(6.12) 
(6.13) 

(6.14) 
(6.15) 

Here ^ is the gauge parameter and FG denotes the Feynman gauge (^ = 0). Since the 
calculation will be performed in the rest-frame of the medium, we have used po in the place 
of p ■ f in Eqs. ( |6.11| ) and ( |6.15| ). In this frame, we define the components of the electron 
four-momentum by writing 



-^p^'p" 

o\Po\/T _ 



1 



(p^ + iey 



'^TrifBip) 



d5{p^ 
dp"^ 



= (fc°, K) 



(6.16) 



Note that this choice of frame will not lead to any loss of generality in the proof of gauge 
independence, since the dispersion relation is Lorentz invariant. 
Let us write 



^T=0 



(6.17) 



We now evaluate the real parts of these functions by putting the expressions for the 
fermion and the photon propagators in the expression for the self-energy given by Eq. 
( |4.93| ). Let us first consider the T = part 



^T=0 



,r=o 



(6.18) 
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and write a 



T=0 _ ^T=0 I ^T=0 



pq'^ + a| ^. Dimensional regularization of Spg" in 4 — e' dimensions 



gives 



T=0 
^FG 



(4vr) 



7 



1 + ln(47r) - 2 j\xxln - + O(e') 



where 7 is the Euler-Mascheroni constant. On the other hand, one easily obtains 



1 



(27r)4 (p2 + 



The odd term vanishes on integration. Evaluation of the rest leads to 



r=o 



(47r) 



--7-I+ 



ln(47r) - dx In - x)k^) + O(e') 



(6.19) 



(6.20) 



(6.21; 



Finally adding the counter-term to Qj , cLS fixed by the renormalization condition 



(6.22) 



we arrive at 

where a is the renormalization scale. 

We now turn to the T-dependent part, and write 



(6.23) 



S — SpQ + Eg. 



(6.24) 



To obtain SpQ, one has to put Eqs. ( |6.10D and ( |6.13D in the expression ( [4. 931 ) for the 
self-energy, and consider the relevant terms. It is then convenient to change p to —p — K 
in the /i?-containing term. Finally setting = 0, as allowed by S{p'^), yields 

d^p , „ , , „ 1 



FG 



2e' 



+ lt)fBip) + ^fFip)]2n6ip')Re■ 



{Q.25) 



(277)4^'"' ' ■'"^•'"^^^ ' /-^-v^^J — ' 2p.k + k"^ + ie ' 

Putting ( I6.10D and ( |6.14D in (|4.93|) , one arrives at S^. Its /^-containing part and fs- 
containing part will be indicated explicitly: 



(6.26) 



In S^^, changing p to —p — K and then setting = (allowed by the delta function) 
gives us 



= / ^Ak'^ - 2p.km^6{p')Mp) Re 



1 



{2p.k + P + ie)2 ■ 



(6.27) 
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In S^^^, simplification leads to 



B 



fBiPj- 



(27r 



dp'^ 



j) - {p^lfe + k'^i)) Re ■ 



1 



(6.28) 



{p + kY + ie 

The term, being odd, drops out on integration. To deal with the remaining part, we 
first write down 



dp"^ 

and then use the regularization 

1 



lim 

f^O 2711 



1 



lim 



d 



1 



so that 



(p2 ± ie)2 A^o p'^ - ± ie 



(6.29) 



(6.30) 



(6.31) 



Let us now make use of Eq. ( |6.31|) in Eq. ( |6.28| ). We then commute the integration over 
Po with the limit and the differentiation involving A , and set p"^ = A^ (allowed by the 
delta function). This gives 

d^P d 



Si 



X Re 



(27r)3 A 



Imi ^ j dpofB{p)S{p' - \'){k'i> + A^ 



(6.32) 



fc2 + 2p.k + X^ + ie 

After integrating over po? the operations involving A are carried out. While this last step 
is easily performed for the part proportional to A^ (since 

d 



lim 



:(A^/(A^)) = /(O) 



(6.33) 



with / denoting a function with finite df/dX"^ at A = ||6ll), a more tedious algebra is 
to be worked out for the remaining part, finally giving 



e2 f d^P 
^17 



1 1 



-k'^P-fo Re 



'l + ko/P l-ko/P\ 



27r)3e/3P _ip 



Re 



1 



PP.,^Refl±^+^-^°/^^ 



D. 



e^l3 f d^P 



(27r)3 {e^P - 1)2 p2 



1 



fc2p7o Re 



1 1 



(6.34) 
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Here 

D± = ±2Pko-2P.K + ie. (6.35) 



While the expressions ( |6.25| ) and (|6.27| ) are the same as those in Ref. (except that we 



have carefully incorporated the ie, which will be needed for the calculations in Appendix 
^ , the expression ( |6.34D is a different one. We have found this form of ^ convenient 
for actually carrying out the integration (this, again, is done in Appendix 

Let us now investigate the high-T behaviour of ( |6.25D , ( |6.27D and ( [6.341 ). Following 
Ref. this behaviour can be inferred from the degree of ultraviolet divergence of the 
integral in each expression in the absence of the cut-off of 0{T) provided by or fp 
or e^^/{e^^ — 1)^. One then finds that goes like T^, while and go like 
T. It may be noted that one cannot obtain the correct high-T behaviour of ^ without 
removing the regulator A, and this can be done only after performing the po-integration, 
as in (PI). 

At high T, one can neglect the renormalized compared to S' . Actually TF^^ 
depends on the renormalization scale a, but since the dependence is logarithmic, we can 
still ignore it vis-a-vis the power law dependences on T in the various parts of S'. We 
shall, however, see that it is not necessary to neglect Tp"^^ for proving gauge independence 
ai K cT. We shall also use ( |6.23| ) to incorporate the a-dependence in the equation for 
the effective mass. 

6.3 Equation governing gauge dependence of dispersion relation 
at one loop, and general considerations 



Inverting Eq. ( p.l| ), one obtains 

1 . 



Tr(I^S) +cuTr(^S)], (6.36) 



b = -^^[u;TTm-{^'-K')TTm]- (6.37) 

The dispersion relation is obtained by putting a, b in Eq. (|6.4|). Let us now write the 
function / in Eq. ( |6.4| ) as 

/ = fFG + h, (6.38) 

/fg = (^-ir)(l-aFG)-&FG, (6.39) 

= -{uj-K)a^-bi:. (6.40) 
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Next using Eqs. ( |6.26|) , ( |6.27|) and ( |6.34|) , one readily sees that k"^ factors out from the 



expression for Tr(|^S^). Therefore in view of Eq. ( |6.3|) , one can write from Eq. ( |6.37|) that 

6^ = (c<j — _ft')(l-loop function of u;, K) (6.41) 

(note that 6, being zero at T = 0, is determined by the S' part only). Consequently 

= _(cj - K)[ai: + (1-loop function of K)]. (6.42) 

This is the equation governing the gauge- dependence of the electron dispersion relation 
at one-loop. This equation is of the form 

(/5)i-ioop = /tree X [Hoop ^-dependent function]. (6.43) 

As we shall see in the next section, the fact that /tree factors out on the R.H.S. of Eq. 
( \6.43i) is crucial to the proof of gauge independence. 



Ref. p4| contains a general, nonperturbative derivation of the identities determining 



the gauge dependence of the gluon dispersion relations (see Eqs. (16) and (17) of [p^ ). 
To one loop these gauge dependence identities are shown to reduce to relations (see Eq. 
(20) of IQ) which are analogous to Eq. ( |6.43D above. It is therefore likely that a general 



gauge dependence identity for the electron dispersion relation, reducing to Eq. ( \6.4^ at 
one loop, can be derived and used to arrive at a general proof of gauge independence of 
the dispersion relation (as in the gluon case). In this paper, however, we confine ourselves 
to one-loop calculations. 

6.4 Gauge independence of one-loop dispersion relation at mo- 
menta much larger than eT 

Since we shall consider ^ eT in this section, let us first note that for sufficiently 
small e there always exists a domain eT <^ K <^T, so that the restriction K ^ eT does 
not contradict the high-T approximation. 

The basic premise for the considerations of this section is that when the finite- 
temperature effects are small, the dispersion relation should involve a deviation from 
uj = K primarily by powers of e^. We have used the phrase "primarily" because the 
terms in the deviation may come multiplied with powers of ln(l/e) (originating from 
powers of ln(a; — K)] see, for example, Eqs. ( |6.5| ) and (|6.6|)). Now, to first order in e^, Eq. 
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( |6.6|) becomes uo ~ K = e^T"^ j%K. This, when used in the remaining term on the R.H.S. 
of Eq. (i;!), results in § 

Mo^ Mo^ 

where Mq^ = The terms on the R.H.S. of the above expansion clearly indicate 

that the expansion is valid for K ^ Mq i.e. K ^ eT. This can be understood in the 
following way. As shown in Ref. |^ and as we shall see in Sec. |6.5.3| , Mq is the leading-order 
effective electron mass and thus can be considered a measure of finite-temperature 
effects. So, K ^ eT ensures that the finite-temperature effects are small (even though T 
is large), and only then the expansion in given by Eq. (|6.44|) is valid. (As we will see 
in Sec. |6.5.1| , such an expansion does not exist for K <^ eT.) 

Now we turn to gauge independence. The expected general form of the one-loop 
dispersion relation for K ^ eT when the terms suhleading in T are kept is 

u = K + e^fi{e,K,T) (6.45) 

where the e- dependence of fi involves only powers of ln(l/e). Equation ( |6.44| ) (without 
the last term on its R.H.S., which is actually of order e^ln(l/e)) is a special case of Eq. 
( |6.45| ), where /i turned out to be totally e- independent. Let us now take Eq. ( |6.45| ) to 
be the relation in the Feynman gauge i.e. assume that it satisfies /fg = 0. Then to prove 



that ( |6.45| ) is gauge independent we have to show that it also satisfies /fg + = 0. It 
suffices to show that, at ( |6.45| ), is of order (with or without powers of ln(l/e)) ||63[| . 
This follows readily from Eq. ( |6.42| ), since the portion of the R.H.S. of Eq. ( |6.42| ) within 
the square brackets involves terms of the order of and e^(ln(l/e))" (n is an integer), 
and so also does {u — K) at ( |6.45| ). 

Digressing briefly from high-temperature QED, we mention that the above demonstra- 
tion of gauge independence is actually sufficient also for the case which involves not large 
T but large momenta. Thus for K ^ T, a one-loop relation of the key form (|6.45| ) should 
still hold. Note that the rest of the proof, namely, arriving at Eq. ( IO2D , did not assume 
large T (in particular, we did not neglect the vacuum contribution to a). The proof of 
gauge independence continues to be valid if, in addition, we have a chemical potential /i 
such that K ^ fi a.s well. Then the relevant changes are that fi in (|6.45|) depends on n 
as well, fp is modified and we have fri—p) in place of frip) in ( |6.27|) , but none of these 
affect the proof outlined above. 
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6.5 Gauge independence of effective mass at one-loop 
6.5.1 Subtleties at zero momentum 

The effective electron mass is the value of uj at K = 0. (Here we are following the 
nomenclature used in much of the literature on high-temperature quantum field theory, 
including Ref. 0, though this quantity has been called the inertial mass earlier in this 
work.) The analysis of the previous section does not prove the gauge independence of the 
effective mass because the form ( |6.45| ) does not hold near K = 0. For example, to leading 
order in T, the one-loop dispersion relation for <^ eT is |^ 

K , , 

.^Mo + - + ^ + ... (6.46) 

(This is obtained by substituting in Eq. ( |6.6| ) a trial solution of the form ( |6.53| ) and 
determining the coefficients.) Eq. ( |6.46| ) suggests that the leading order values of a and 
b are not of 0{e^) (with or without powers of ln(l/e)) in the K <^ eT limit; indeed 
the values at i^' = are a = —1/3 and b = —2Mq/3 [Q. This apparently surprising 
behaviour can be understood in a simple way, without performing a detailed calculation, 
as follows. 

At leading order a and b are e^T^ times some function of u and K (T^ being deduced 
from ultraviolet power counting in S'). Therefore, to have the correct dimensions, aK=o ~ 
e^T^/Mo^ and 6k=o ~ e^T^/Mo (remembering Mq = ujk=o at leading order). Putting 
these in the dispersion relation ( |6.4| ) at K = 0, namely. 



Mo - MoaK=o - bK=o = 0, (6.47) 
gives us Mq ~ eT. Using this, it follows that aK=o is of 0(1) and 6x=o is of O(eT). 
6.5.2 Formula for the effective mass 

The traces Tr(|^S) and Tr(^S) are functions of uj and K. The expressions for them 
are obtained from Eqs. ( |6l8D , (|]23D, (|]25D, (Ip^) and ( lOiD . First of all we show that 



these traces are both even functions of K. For the S-^"° part, this is obvious from ( |6.3|) . 
For the S' part, first note that k'^ can occur in the expressions for the traces only through 
k ■ p, k ■ V and k"^ {p^ being the integration variable). Next, let us go to the rest frame 
of the medium (in which K is the magnitude of the three-momentum) and note that K 
can now occur in the expressions for the traces only through K.P and K"^ . Odd power of 
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K can come from K.P = KP cos 9 where 9 is the angle between K and P. This is also 
the only place where 9 occurs in the integrand. So changing 9 to n — 9 together with the 
change K —K, we establish the even nature as a function of K. 
For small K, therefore, we can write 

1 



Tr(I^E) 



4 

I Tr(^S) 



ho + hiK' + h^K" + ■ ■ 

-2 



(6.48) 

go + giK' + g^K^ + ■■■ (6.49) 
where hi and (^j are functions of u). Substitution of the above expressions into Eqs. (|6.36|) 



and (|6.37|) gives 



a 
b 



-^{ho - ^go) - {hi - ujgi) + 0{K^) 

UJ 



{ho - togo) + {ujhi - co^gi + go) + 0{K^). 



(6.50) 
(6.51) 



For the form factors a and h to remain analytic at -fC = 0, the relation 



ho = ugo (6.52) 
must hold. We now put ( p.50| ) and ( |6.51| ), subject to the constraint ( |6.52|) , in the dispersion 



relation ( |6.4| ). Then we use the fact that for any small K, the solution 

M + uiK + U2K'^ H 



UJ 



(6.53) 



where M and uJi are constants, must satisfy the dispersion relation. We also expand hi 
and gi in K by first doing a Taylor expansion of them around = M as functions of tu, 
and then putting ( |6.53| ). Thus 



hi{uj) = hi{M) + h[{M)uJiK + 0{K'^ 



(6.54) 



with a similar equation for g^. The equation that now results from the dispersion relation 
( |6.4D has a series in powers of K alone on the L.H.S. Equating the coefficient of each 
power of K separately to zero will determine the constants M and uji of ( |6.53| ). Thus, 
equating the constant term to zero gives 



M = go{M). 
Multiplying both sides by M, we make use of 

ho{M) = Mgo{M) 
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(6.55) 



(6.56) 



(which follows from ( |6.52D ) to finally arrive at 

= hoiM). 



(6.57) 



One can use either (|6.57| ) or (|6.55|) as the formula for the effective mass. We shall use 
( 16.571 ), which, by using Eq. ( |6.48| ), can be written as []65| : 

1 



= - lim Trim . 

4 K^O,ui^M ^' ' 



(6.58) 



We comment that one should not attempt to compare this formula with the formula for 
the inertial mass of Sec. [4.1| , since the latter was deduced by assuming a straightforward 
expansion in and neglecting various terms accordingly. 

6.5.3 Investigation of gauge independence 



Let us now substitute the expression for the self-energy in Eq. ( |6.58|) . On using Eqs. 
( CT ) and (|08D in Eq. (|638|) , we obtain 



M'^ = - lim Tr[|^S'l 



/2 



(6.59) 
(6.60) 



At T = 0, M' = and so Eq. (|09D is correctly satisfied by M = 0. 

Putting in Eq. ( |6.60[) the expressions for SpG, ^ and S^^^ given by Eqs. (|6.25|) , (|6.27|) 
and ( |6.34|) , and noting that the limit u ^ M translates to /cq — ^ M in the rest-frame of 
the medium, we obtain 



/2 



M'r 



FG 



2e' 



[{po + M)fB{p)+PofF{pmp')Re 



2po + M + ij^ 



M' 



d'^p 



(27r)3 

-Ml 

ie^M f d^P 1 1 



Po/f(p)5(p^) Re 



{2po + M + tj^r 



M' 



2 J {2TifPePP- 
ie^M f d^P 1 



/ 1 1 



Re 



2 J {2Tifei^P-l 
ie^M^(3 f d^P 1 e^^ 



(27r)3 P (e 



f3P 



P \E+^ EJ^ 



E_' 



(6.61) 
(6.62) 

(6.63) 
(6.64) 
(6.65) 

(6.66) 
(6.67) 
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Here 



(6.68) 



We now investigate the high-T behaviour of the three terms on the R.H.S. of Eq. 
( |6.61| ). This can be arrived at from the ultraviolet behaviour, as explained towards the 
end of Sec. E72. However there is a constraint that each term on the R.H.S. of Eq. ( |6.61| ) 
is an even function of M. (This can be seen by changing pq to —po together with the 
change M —M in ( |6.62|) , in ( |6.63| ) and in 

1 



Ml 



(6.69) 



the last equation having been obtained by putting ( |6.28| ) in ( |6.6CI| ).) These considera- 
tions, plus dimensional analysis, tell us that at high T, MpQ goes like T^, M'^ / like 
M2(ln|T/M|)"i and M'^^ like M'^{\n\T /M\Y^ (ni, n2 being positive integers). Since 
after T^, the next allowed term in MpQ^ is M^(ln |T/M|)"^ {n^ a positive integer), the 
general expression for M'^ at high T can be written as 



M' 



/2 



CoT^ + csM^ In 



M 



^M^ ci In 



T' 
M 



"1 / T ^ 



n2N 



(6.70) 



where cq, ci, C2, C3 are constants, and we have dropped the modulus of the arguments of 
the logarithms since T/M is positive. Terms independent of T have been neglected on 
the R.H.S. of Eq. ( |6T0D . 

In principle there could be further constraints ruling out some term(s) in Eq. ( |6.70| ), 
but a detailed calculation, described in Appendix 0, reveals that all these terms are indeed 
present and that rii = 1 = n2 = n^- Thus we actually have the high temperature equation 



T2 T T 
\ ^ In h £ — ^ In — 

8 87r2 M ^87r2 M 



(6.71) 



The detailed calculation may also be viewed as a check on the argument involving the 
degree of ultraviolet divergence mentioned before. 

Using Eqs. ( |6.23| ) and ( |6.71| ) in Eq. ( |6.59| ), we then have 

-y2 



+ 1 



(6.72) 



On the R.H.S. of Eq. (|6.72|) (as also in Eq. ( |6.71| )), we have not given the term ~ e^M^; 
all other terms contributing to the R.H.S. of Eq. ( |6.72| ) (and (|6.71|) ) vanish in the limit of 
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large T/M. It is interesting to note that in Eq. ( p.72| ), the \n{M/ij) term from S^"° and 
the ln(T/M) term from S' have exactly combined to yield just a ln(T/cr) term, because it 
has been observed that similar combination also takes place in the case of the gauge boson 
self-energy in the Yang-Mills theory |^|. A discussion of similar behaviour in the case of 



three-point function, and, in general, A^-point function of gauge boson in the Yang-Mills 
theory, is to be found in Ref. p7|] . 



The first important observation from Eq. ( ti.72| ) is that = to leading order 



in T, which, of course, is a well-known result Now, the gauge dependent part of 
S' goes like T. Therefore in view of Eq. (|6.58|) , one would expect an e^^MT term in 



Eq. (|6.72|) . The absence of such a term shows that the part ofH'^ leading in T does not 
contribute to M. 

To quantify the effect of this absence let us consider the equation 

= e\— + cMT) (6.73) 
8 

where c is a constant. To leading order in T, the second term on the R.H.S. is negligible, 
so that M = 0{eT). This can now be used as an approximation in the second term, to 
give 

= 6^ — [l + 0(e)] (6.74) 
8 

showing that there is a correction of O(e^T) to M. (This conclusion can also be arrived 
at by solving Eq. ( |6.73| ) for M.) So the absence of the e^MT term means that there is no 
0{e^T) correction to the effective mass in any linear covariant gauge. It may be noted 
that the consequences of the presence or absence of various powers of T in Eq. ( |6.72| ) are 
to be taken seriously despite the presence of the scale a, since the a dependence is only 
logarithmic and is not multiplied with any power of T. 

Finally, ( |6.72| ) shows that only the subleading InT dependence of contributes to 
M. Using the leading order result M = 0{eT) as before to approximate the remaining 
terms on the R.H.S. of Eq. (^1% , we infer that the ^dependence in Eq. §^ is ©(e^T^), 



apart from the logarithm. As one certainly expects O(e'^T^) contribution to Eq. (|6.72|) 
from two loops, it is possible that the ^-dependence which we have obtained will be canceled 
by two-loop contribution. But it seems that this issue can be decided only by an actual 
two-loop calculation. It is however clear that an e^^MT term, being of O(e^T^), was less 
likely to be canceled by two-loop contribution. This shows the significance of the absence 
of such a term. 
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CHAPTER 7 



Conclusions 

In this thesis, we studied some aspects of particle physics in the presence of a back- 
ground medium, using the real-time formulation of thermal field theory. We studied the 
O(e^) corrections to the gravitational interactions of a charged lepton in the presence 
of a matter background. We also investigated gauge independence of one-loop electron 
dispersion relation in high-temperature QED. These calculations complement and extend 
previous calculations along similar lines, in various useful ways. The usefulness, as we 
will explain, is both theoretical and phenomenological. 

Wc first discuss the importance of the various aspects of our work from a theoretical 
point of view. The calculation of gravitational couplings in a medium, as presented by us, 
employed several thermal field theoretic techniques that can be useful in other contexts 
also. For example, we made explicit use of only the IPI diagrams in our calculation 
of the gravitational mass. The effect of the one-particle reducible diagrams were taken 
into account by choosing each external spinor to be the solution of the effective Dirac 
equation in a medium, and multiplying the spinor with the square-root of the wavefunction 
renormalization factor. It was then seen that the wavefunction renormalization factors 
arc instrumental in cancelling an infrared divergent contribution in the IPI gravitational 
vertex function. 

Again, we dealt extensively with the limiting case of the finite-temperature gravita- 
tional vertex function in which the graviton carries zero momentum. The ambiguity of 
the finite-temperature Green's functions evaluated at zero momentum is a well-known 
problem in thermal field theory. This property is usually due to the fact that the different 
mathematical limits correspond to different physical situations, so that the resolution of 
the apparent paradox lies in recognizing the appropriate correspondence with the physical 
situation at hand. In our case also, the IPI gravitational vertex at zero graviton momen- 
tum turned out to contain contributions which superficially seemed to be ill-defined, and 
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the problem was resolved by taking the relevant physical situation into account. Thus, 
taking the various limits in the precise order indicated in the operational formula for 
the gravitational mass, given by Eq. ( |5.25| ), we obtained unique and well-defined expres- 
sions for all the contributions. We also reproduced these expressions by using the special 



On the other hand, our work on high-temperature QED gives credence to the idea 
that any gauge dependence present in the self-energy should not show up in the dispersion 
relation. This idea motivated us to focus on the subleading part of the self-energy of the 
electron at one- loop, which (unlike the leading, O(T^) part) depends on the choice of 
gauge. 

We obtained an equation which governs the gauge dependence of the one-loop disper- 
sion relation, stressed its analogy with the corresponding equation in the gluon case, and 
consequently pointed out the possibility of a generalization of our equation to all orders. 
From this equation obtained by us, the gauge independence in the K ^ eT limit followed 
in a straightforward way. We then showed that the effective mass is not affected by the 
leading part (going like T) of the gauge-dependent part of the self-energy and hence does 
not receive O(e^T) correction in any gauge. While the effective mass was found to be 
influenced by the subleading part (going like InT) of the gauge-dependent part of the 
self-energy, it is possible that this will be canceled by two-loop contribution. Exploration 
of this possibility can be a problem for further investigation. 

Our calculations and results are also important from a phenomenological point of view. 
We found that the contribution to the gravitational mass of the electron from a matter 
background with a non-zero chemical potential, such as the Sun, can be quite appreciable. 
An interesting problem, combining the various aspects covered in this thesis, would be 
the study of the gravitational interactions of charged leptons in a background at high 
temperature. The results can then be estimated for the case of a supernova, which is an 
example of such a background. 

Moreover, the matter-induced corrections to the gravitational mass are different for 
the various charged lepton flavors, and are not the same for the corresponding antiparti- 
cles. There are situations in which mass differences, intrinsic or induced, have important 
physical implications, such as the neutron-proton mass difference in the context of the 
nucleosynthesis calculations in the Early Universe. Although our work has focused in the 
case of the charged leptons, similar considerations can be applied to the other fermions 
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as well. Our calculations have provided a necessary ingredient for being able to consider 
them in a systematic manner, and set the stage for their further study on a firm basis. 
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APPENDIX A 



Field quantization using discrete 

momenta 



In this Appendix, we relate the quantization with discrete momenta, used in Sec. 2.2 



to the quantization with continuous momenta, used in Sec. 2.3 onward |68 



Beginning with the complex scalar field, we consider the Fourier expansion given in 
Eq. ( |2.19[ ). The equal-time commutation relations then force 



[a{P), at(P')]_ = 5(3)(p _ p/) ^ [^(p)^ b\P')]_ = 6^^\P -P') 



(A.l) 



all other commutators being zero. Also, Eq. ( |2.19D leads to 

H = j (PPEp(a\P)a[P) + b\P)b{P)) , 



(A.2) 



If we instead use quantization in a volume V, the plane-wave states can only possess 
the discrete momenta 



P 



27r 

\n/3 



n . 



(A.3) 



where rii, n2, are integers. Integrals are then replaced by sums 

d^P 1 



and delta functions by Kronecker symbols 



V 



E 



(A.4) 



(2-Kf5^^\P-P')^V5pp, 



(A.5) 



Note that the left-hand sides of Eqs. ([A.4|) and ([A.5|) taken together give unity, and so 
also do the right-hand sides. 
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Let us now redefine the ladder operators through 



V 



{2n 



i2n) 



3 "P 



(A.6) 



Eqs. O, (^) and (Q) give 



(A.7) 



while Eqs. (Q), ([A^) and O give 



(Alternatively, one can use Eqs. ( |2.19| ), ( |A.4| ) and ( |A.6| ) to arrive at 



(A.9) 



the discretized Fourier expansion for the field. Eqs. ( |A.7| ) and ( |A.8D then follow from this 
expansion.) Eqs. ( |A.7| ) and ( |A.8D are the key expressions used in Sec. |2l^ . 

Turning next to the fermion field, we consider the Fourier expansion of Eq. ( p.37| ) and 
the spinor normalization given in Eq. ( ^.3^ ). The equal-time canonical anticommutation 
relations then force 

[cs{P),cUp')U = Sss'6^'\P-P') , [4(P),4(^')]+ = Sss'6^'\P-P') , (A.IO) 

all other anticommutators being zero. (In more general terms, the choice A^i = 1 = A''2 
forces A/'g = 1 in Eqs. ( pl6|) - ( [BlSl) .) Also, Eq. (^1371) leads to 



H = d'PEpY.{c\{P)c,{P) + d\{P)d,{P) 

•' s 

Q = Jd^PY.{cl{P)cs{P)-dl{P)ds{P) 



(A.ll) 



Use of Eqs. ( [A.4| ) and ( |A.5| ), and redefinitions similar to those in Eq. ( |A.6| ) then give the 
fermionic counterparts of Eqs. ( [A.7| ) and ( |A.8| ), which were used in Sec. pT2 . 
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APPENDIX B 



Full propagator and wavefunction 
renormalization factor in vacuum 

We start with the definition of the full fermion propagator in vacuum 

is)^p{x) ^ {mMx)Mm^) (B.i) 

= ^(xo)(0|7/;,(x)V^^(0)|0)-^(-xo)(0|V^^(0)V^,(x)|0) (B.2) 

where ip{x) is the Heisenberg picture field, and |0) the vacuum of the full interacting 
theory. The wavefunction renormalization factor Z is defined through 

(0|V^40)|/(P,s)) = v^M,«(P), (B.3) 

where the 1-particle state |/(P, s)) is again an eigenstate of the full interacting hamilto- 
nian. The tree-level spinor Us{P) is actually the analogue of Us{P) (see Eq. ( [4.26| )), since 
the sum of the bare mass and the self-energy correction in vacuum is set equal to the 
physical mass of the fermion. 
Eq. ( |B.3D implies that 

(0|^/^„(x)|/(P,s)) = ^Usa{P)e-'^-\ (B.4) 

with po = Ep = \J P"^ + '^z- (We used an identity similar to the one in Eq. (|5.6| ).) It then 
follows that 

(/(P,s)|V'a(x)|0) = y/ZuUPy"-'' ■ (B.5) 

Similar relations involving the 1-antiparticle state are |/(P, s)) are 

(/(P,s)|V'„(a;)|0) = v^t;,,(P)e^f-, (B.6) 
(0|7/>«(x)|/(P,s)) = VZvUPy-'""' . (B.7) 
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In Eqs. ([B.6|) and ( |B.7|) also, po = Ep. It may be noted that in vacuum, Z is independent 
of momentum and also the same for the fermion and the antifermion, properties which 
may not hold in a medium. 

We shall take the normalization of 1-particle states to be 



{f{P\s')\f{P,s)) = {2r:f5^'\P~P')5,,,. 



(B.8) 



with a similar normalization for s)) as well. The completeness relation is then given 
by 



r SP ^ ^ r <fP - ^ - ^ 

i = |o)(o| + E j ^l/(^.^))(/(^'^)l + E j + 



(B.9) 



where the ellipsis stands for multiparticle state contributions, which we shall not be 
interested in. 

Let us now insert Eq. (|B.9|) in each of the two terms of Eq. (p.2| ). A look at Eqs. 
( P.4|) -( pl7| ) then reveals that the 1-particle intermediate states contribute to the first term 
in Eq. (|B.2|) and the 1-antiparticle intermediate states to the second term (there is no 
contribution from |0)(0| since the fermion field cannot have a vacuum expectation value). 
Thus, 



r d^P ^ ^ - 

= ^ME / 7^(0|^.(x)|/(P,s))(/(P,.)|V^,(0)|0), (B.IO) 

1-particle g •' [^T^j 

r d^P - - ^ - ^ 

-^(-^o)E j T^,mM\Kp,s)){f{p,s)\Mxm . (B.ii) 



1— antiparticlc g J (2vr) 

In Eq. (IRToD , use Eqs. and (|R5D , and also 



Changing to po — Ep then leads to 



2-71 Po + ie 



d'p ^_,^.,ZEsUs{P)Us{P) 



[B.12) 



(B.13) 



l-particlc {2tt)'^ Pq - Ep + it 

In Eq. (|B.11|) , use Eqs. ( p.6| ) and (p.7| ), and also put 6{—xq) from Eq. ( p.l2|) . Changing 
Po to — Po — Ep and P to —P then leads to 



1— antiparticle 



d'^p 
(2^ 



-ip-x 



ZEsVs{-P)Vs{-P) 

Po + Ep- ie 



(B.14) 
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We are now going to check Eqs. ( |B.13|) and ( |B.14| ) at the tree- level, for which Z = \. 
We first note that Eqs. ([B^ID and (IKSD imply that 



«1(PK(P) = 1. 



[B.15) 



(This can be seen in the following way. Let us introduce the factors A^^i, A^2, ^3, ^4 
through 



Cs{P)u,{P)e-'P-^ + 



\f{p,s)) = jN,ci{pm 



(B.16) 
(B.17) 

(B.18) 
(B.19) 



omitting the antiparticle parts for brevity. Then, [ip{x,t),ilj{y,t)]_^_ = S^{x — y) gives 
N^N-i = N2, Eq. (^) (with Z = 1) gives N^^'^N^ = (27r)3iV2, and Eq. (gj) gives N^N^ = 
(27r)3. It then follows that A^i = 1.) 

Similarly, vl{p)vs{P) = 1. For these normalizations, the spin-sum relations are given 
by Eq. (|^). Since 



SfiPo, P) = —2 ^2 



(B.20) 



pl-El + ie 2Ep {po - Ep + ie) 

near the pole po = Ep ~ "^e, Eq- ( [B.13| ) (with Z = 1) checks out. Similarly one can check 
Eq. 



In a medium , Eqs. ( B.13| ) and ( [B.14 ) are generalized to Eqs. ( |4.38 ) and ( |4.52| ) respec- 
tively, while maintaining the relations given by Eq. ( [B.8| ) (see remark after Eq. ( p.l8[ )) 



and Eq. (|B.15|) (see Eq. ( [4.27] )) which do not depend on the wavefunction renormalization 
factor. 
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APPENDIX C 



Charge conjugation and the IPI 
gravitational vertex 



The aim of this Appendix is to derive the charge conjugation property of the IPI 
gravitational vertex, given by Eq. ( |C.20| ), or Eq. ( |U.21| ), below. 

We begin by considering the case of a C invariant Lagrangian and a C symmetric 
background. Since /i = 0, Eqs. (|2.2|) - (|2.4|) reduce to 



Tr e-l^H 

Using = 1, this can be rewritten as 



Use of CHC = H (see Eq. (g)) and = 1 then gives 

(O) = (COC) . (C.3) 

It will be instructive to first apply this equation to the definition of the full fermion 
propagator S'f{p), given by 

^{2nr5^'\p-p')S'J^M = J d'y dh e^^-ye-^^'^ {T[My)i'mi^)]) , (C.4) 

where the fermion fields are the Heisenberg picture fields. (This is equivalent to the 
definition 

{T[My)^mm = ^ J (0e-^^-(^-^)^}_(A;) (C.5) 
used elsewhere in this work.) Using Eq. ( |(J.ci| ) and = 1, we can write 

{T[My)i^m{^)]) = {T[CMy)c Ci>Uz)C]) . (c.6) 
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Now, CipC is given by Eq ( [4.60| ). Also, CifjC = {CipCY'-yo, since (as explained below Eq. 



(g)) C is hermitian. So, using Eqs. {^^^, and ([O^ ), we obtain 

C^jC = -V^^C-i . (C.7) 

One can now use Eqs. ([4.6CI|) and ( |C.7|) to determine the R.H.S. of Eq. (|C.6|) . Consequently, 
the R.H.S. of Eq. (CJ) equals 



The double integral is just i{27c)^5^'^\p — p')S'j-^i^{~p), by the definition (|C.4|) . So the 
expression given in Eq. ( p.8| ) equals 

z(27r)^5W(p - p') [CS']{-p)C~'U . (C.9) 

Since this equals the L.H.S. of Eq. ( |C.4| ), we arrive at 

S'f{p) = CS']i-p)C-'. (C.IO) 

As a check, we note that the use of Eq. ( [4.1| ) reduces Eq. ( |C.10| ) to Eq. ( |4.74| ), which was 
deduced by demanding the action to be C invariant. 

Now we shall apply Eq. (|C.3D to the definition of the IPI gravitational vertex function 
T^P{p,p'), given by 

d^xd^yd^ze-'''-''e'P-ye-''p'-'{T[f^P{x)^n{y)^rn{z)]) . (C.ll) 



Using Eq. (^), = 1, and CT^pC = T^p, we can write 

{T[f^p{x)My)^m{z)]) = {T[f^'{x)CMy)cci>Uz)C]) . (c.12) 

One can now use Eqs. ( ^4.60| ) and ( \C7\ ) to determine the R.H.S. of Eq. ( |C12D . Conse- 



quently, the R.H.S. of Eq. (|C.11D equals 



Cnb{C-')am / d'x d' z d'y e'^'^-^ 6-'^' c'^'y {T[f^P {x)Mz)My)]) . (C.13) 



The triple integral is just {2'k)^6^^\p - p' - g)5}„„,(-p')r^f„,(-p', by the 

definition So the expression given in Eq. (|C.13| ) equals 

(2vr)^5W(p - p' - q)[CS']{-p)T'P^{~p', -p)S']{-p')C-'U • (C.14) 
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Since this equals the L.H.S. of Eq. (|C.ll ), we arrive at 



S'f{p)T,,{p,p')S'j{p') = CS']{-p)Tl{-p\ -p)S']{-p')C-' . (C.15) 

Using Eq. ( p.lCI| ), this leads to 

T,,{p,p') = CTl{-p\-p)C-K (C.16) 

Now, consider the case of C symmetric Lagrangian and C asymmetric background. 
Then 

instead of Eq. Since CQC = — Q (see Eq. ([4.81| )), we now have 

(O) = [{COC)]„ (C.18) 
in the place of Eq. ( |C.3| ) . It is easy to see that this results in 

S'fip) = [CS']{-p)C-']„, (C.19) 

r,,(p,p') = [ctU~p\-p)c-']„. (C.20) 



Since fi enters through the thermal propagator of the fermion given by Eqs. ( |2.50D and 
( p.51| ), and /i —fi in the propagator is equivalent to v ^ —v, one can also write Eq. 
( 10201) in the form 

rxpip,p',v) = CTli-p', -p, -v)C-' , (C.21) 
where we have explicitly indicated the dependence of the vertex function on the vector 
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APPENDIX D 



Transversality of the gravitational 

vertex 



It is useful to verify that the complete vertex function to O(e^), obtained in Sec. |5.4| , 
satisfies the transversality conditions 

q^UiP') T,,{p,p') U{P) = = q^UiP') T,,{p,p') U{P) (D.l) 

to this order. The wavefunction renormalization factors associated with the external lines, 
being scalars, do not play any role in these relations (in other words, they just cancel out). 
Also, since the vertex is symmetric in the Lorentz indices A,p, either of these relations 
guarantees the other. In order to simplify the notation, in this Appendix we omit the 
subscript s in the spinors. 

In order to verify this relation, the important point is that we must include all the 
terms upto O(e^). Since the one-loop terms in the induced vertex are already O(e^), for 
them we can adopt the tree-level definition of the spinors, i.e. 

^u{P) = mu{P) , m(P')/ = mu{P') , (D.2) 

as well as the tree-level on-shell conditions 

In this appendix, as well as in Appendix |^, we use the tree level mass m without any 
subscript, implying m^, rUe or mj which should be understood from the context. Also 
note that the photon distribution function as well as the associated 5-function are even 
in k, and therefore those terms which are odd in k in the rest of the integrand do not 
contribute. 
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We first show that the vertex contribution from Fig. 5^ A is transverse by itself. From 
Eq. (|31, 



q^Axpa{k, k - q) = -{q^ - 2k ■ q){Akpka + k ■ qrjpa) 



(D.4) 



where we omit the terms that are proportional to q^ because, in Eq. (|5.49|) , they will yield 
^ which vanishes between the spinors. Changing the sign of q in the last equation yields 



q^Axpa{k + q,k) = {q^ + 2k ■ q){Akpka - k ■ qr]pa) 



(D.5) 



and as a result q^u{P')T^^\p,p')u{P) turns out to be proportional to I]/<5/(^/ — n^), 
which is zero to this order. 

As for the other diagrams, straightforward algebra gives the following results: 



qMn^T\p,p>ip) 



d^k 

4 J {2^ 



Kk')v,{k) 



X u{P') 



Amkp - (p + 5p')p^ + 2k ■ p^ip , a 
k^' V'^'P) 



4 

X u{P'] 



(27r)3 

4(^ - 2m)kp + ^{p + p')p - 2k ■ p'^ip 



— k ■ p 



q^u{P')T'TM)u{P) 



X u{P') 



{p ^ p) 
d^k 

mkr 



uiP), 
S{k')v,{k) 



mkp - ftpp 



k ■ p' 



k ■ p 



u{P) 



q'u{P')TiP{p,p'HP) 



— e 



d^k 



5{k'^ -m^)r]f{k) 



X u{P') 



(27r)3 

{f - 2m)kp + mpp {jt - 2m)kp + mp'^ 



m? — k ■ p 



m? — k ■ p' 
1 \ 



u(P), 



(27r) 

X u{P')hqp + k-q-ip]u{P) 



k ■ p' k ■ p 



q^u{P')vf;^^'\p,p')u{P) 



^2 d^k 



{2nY 



5{k^ -m^)7]f{k) 



X 



1 1 

+ 



m? — k • p' m? — k ■ p J 
X u{P') \{lt - 3m)qp + k ■ q-fp] u{P) 



(D.6) 
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Therefore, adding all the one-loop contributions to the vertex, we obtain 



(2vr)= 



X u{P- 



k ■ p' 

2 " 'd^k 



Xp 



7 (3pp - p'p) - [p^p' 
6{k^ - m^)7]f{k) 



u{P) 



(D.7) 



X uiP' 



4 J (2Txf 

- 2m){3p' - pp) - 2k ■ p-ip 



im? — k ■ p 



- (p ^ P') 



u{P). 
(D.8) 



We need to add to these the tree-level contribution to the gravitational vertex that 
appears in Eq. ( p. 951 ). In this case, we must include the O(e^) corrections to the equation 



for the spinors, which arise from the self-energy diagrams of Sec. [4.5| . Thus, for this part, 
using Eq. (|4.26| ) and its hermitian conjugate 



U{P')U' -m-T.{p')) =0, 



we obtain 



U{P')iU{P) = U{P')(^\p) - ^\p'))U{P) , 



(D.9) 



(D.IO) 



which in turn yields 

q'^U{P')V^p{p,p')U{P) = \u{P') 



{3p' - p)p^'{p) + p^-ip - (p ^ p') 



U{P). (D.ll) 



This can be cast in a different form by multiplying Eq. ( [4.26|) from the left by U{P')'~fp{'^ + 
m) and Eq. ( p.9|) from the right by + m)'jpU{P) and taking the difference of the 
resulting equations. This gives 

{p' - P") U{P'hpU{P) = U{P') hp{i> + m)S'(j9) - + m)^p] U{P) , (D.12) 



and substituting this result into Eq. ( p.ll|) , we obtain 



q^UiP') Vxpip^p') U{P) = -U{P') iSp' - p)pj:'ip) - (Sp - p')pS'(p') 



U{P). (D.13) 



Since S' is explicitly of O(e^) while we are interested in results to O(e^) only, we can use 
the tree-level spinors on the right-hand side. Using Eq. ( p.3|) in Eq. ( |4.95| ), we can write 
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the self-energy contribution involving the photon distribution function as 



disregarding terms odd in k. Similarly, from Eq. ( [4.96D , the part containing the Fermi 
distribution function can be written as 

Substituting these forms into Eq. ( p.l3| ) and using the identities 

u{P'hpi:^ + m)^u{P) = 2k-pu{P')-ipu{P) 

u{P')lt{f + m)-fpu{P) = 2k ■ p' u{P')-fpu{P) , (D.16) 

we see that Eq. ( |D.ltj| ) cancels the contribution from the loop diagrams given in Eqs. 
( p. 71) and (|D.^ ) to this order. This proves the transversality of the effective vertex. 
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APPENDIX E 



Apparently ill-defined contributions 
to the gravitational vertex 



E.l The Bl contribution 



We start from the formula given in Eq. ( p.43| ), from which it follows that 



V(B1) 
Ap 



where we have put 



- 17 Ijm 



k-p{2K-Q-Q^ 



Yjj) + m)YC^^xpik + g, k) 

{k-p + K-Q){2K-Q + Q^) 



(E.l) 



(E.2) 



In order to take the limit Q ^ 0, our strategy is to expand the coefficients of the factors 
1/{2K ■ Q ± Q"^) in powers of Q. Of the resulting terms in the coefficients, those which 
are quadratic in Q do not contribute in the Q — ^ limit and therefore we need to keep 
only the terms that are at most linear in Q. 

Using the property C^y\p{k + q,k) = Cy^\p{k, k + g), we can write 



Cfj^u\p{k, k — q) — C^uXpik, k) — C'^^^p{k, q) 
CfivXpik + q,k) = Cp_yXp{k,k) + C^^)^p(k,q) 



(E.3) 



where 



C'^uXpi.^^ Q) = VxpiVpiyk ■ q - q^,k^) - r]p^{kxqp + qxkp) + k^{r]xpqp + VppQx) 

+q^{Vxukp + r]pykx) - k ■ q{r]Xf,r]p^ + iixuVpp) ■ (E.4) 
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To first order in Q, we can also put 

1 _ 1 K-Q 

k-p + K -Q' k-p {k-pY' 

This enables us to decompose T'^^^^ {p, p) in the following four terms: 



p/(Bla) 



X 



1 



_2K-Q-Q^ 2K-Q + Q^_ 



— lim 

2 Q^O 



(2^ 



5(A;2)r/,(A;) 



C,,.Ap(fc,fc)7"^7^ 
k-p{2K-Q + Q^) 



(E.5) 



-./(Blc) 
Ap 



r 



/(Bid) 



Ap 



{P) 



— lim 



(2vr) 



Cp.Ap(fc,fc)7'-(|^-|^ + m)7^ 
(A;-p)2 



— lim 

2 



X 



K-Q 



2K-Q + Q^ 



(E.6) 



2K-Q-Q^ 2K-Q + Q^_ 
We carry out these integrals one by one. 

Eliminating the manifestly k-odd terms from the integrand and performing the kQ- 
integration, we obtain 



Tif'^\p) = e'hra (■ 



Akxk 



1 



(E.7) 



.0^ {2n)^2K k-p [2K-Q-Q^ 2K ■ Q + Q\ 

using k^ = 0. The expression within the square brackets is finite for Q —>■ 0. Therefore, 
in the spinors we can set p = p', and using Eq. ( |4.58| ) we then obtain 



m 



(Bla) 



4e2 , f d^K 
lim / - — -- 

m q-^qJ [2txY 



f,{K)K 



_2K-Q-Q^ 2K-Q + Q\ 



(E.8) 



We can perform the integration over the angular variables in K, the integral being un- 
derstood, as usual, in terms of the principal value part. That gives 



dn- 



2K-Q-Q^ 



dVL 



TX 



SO that 



(Bla) 



TT^m 



6m 



dK f\{K)K 
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(E.9) 



(E.IO) 



As for the next contribution, it is straightforward to verify that 



(E.ll) 



using q ■ V = qo = 0. So 



'-^ + 2k-vi)) 
k ■ p 



(E.12) 



Now using Eqs. (|4.58| ) and (|4.59|) , carrying out the integral over ko, and finally putting 

P = 0, we get 



m 



(Bib) 



Similarly, 



and 



(2t;V-r^^0r'ir^(p) 



m J {27t)32K'^'^ ' 12m' 

f,{K)Y{^ + m)YC,,,,{k,k) 



2 J {2tiY2K 

2e2 r cPK 



{k ■ py 



(E.13) 



(E.14) 



X 



m J (27r)32K ^^^^\k ■ pf 
-{k ■ pf - 2m\k ■ + 4:{k ■ p){k ■ v){p ■ v)] (E.15) 



so that 



2e' 



2rp2 



(E.16) 



(^^■^^ m J {2Tif2K 12m 

For T'I^^'^^ we first perform the integral over k^. Remembering that in the remaining 
integral we can change K to —K and using the fact that C'^j^xp{~^i l) = ~C'^i/Ap(^) i)j 
obtain 



^ Ap 



lim 



2 {2'nY2K''''"' \k-p^ 



2K-Q-Q^ 2K-Q + Q^_ 



K.x,ik,q)-Cl^,^{k,q)] 



lim 



— — + — — 

2K-Q~Q^ 2K-Q + Q^ 

'''' ' 



[C[,u\p{k,q) + Cl^xp{k,q)] 



2 Q^oJ (27r)32ir'^' 'k-p 



X • 



7 



K-Q 



C',.xpik,q)+Cl^,^{k,q)\+0{Q) 
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(E.17) 



We now use 

(2t;V - v'nrh'[C',.x,{k, q) + Cl^x,{K q)] = 8{k ■ v){K ■ Q)f . (E.18) 

Then, using Eq. (|4.58| ) and putting P = 0, we get 

4e2 r d^K 



(E.19) 



^^^^^ m J {2tt)^2K-' Qm 
Adding the results given in Eqs. (lElOl ), (|El3|) , and ( ^lOD , we get the total 



contribution from the Bl term presented in Eq. (|5.63| ). 



E.2 The A2 contribution 



For this contribution, we start from Eq. ( 5.40 ). Using Eq. ( [E.2|) , we can write it as 

-2 f d^k 
m / ■ 

Q 

Axpik + q, k) 



2 q^qJ U-k)'^ 



X 



(27r)3 
AAp(fc, k~q) 



X2K-Q- Q2)(^2 ^k-p) {2K-Q + Q'^){m? - k ■ p') 
Following the strategy stated below Eq. ( |E.2| ), let us now write 



(E.20) 



Axp{k,k~q) = Axp{k,k) + A\Xk,q), 
Axp{k + q,k) = Axp{k,k)+Al^{k,q), 

and, in the denominator, expand — k ■ p' in powers of Q: 

I = 1 + ^'-g +0(Q^). 

— k ■ p' 711? — k ■ p {m? — k ■ pY 

Then we can decompose T'l^'^'^\p^p) into the following terms, omitting higher powers of 
Q which anyway will not contribute: 



(E.21) 



(E.22) 



V(A2a) 



e2 , r d^k 
— lim / - — — 

2 Q^J (27r)3 



X 



di.k'^ - m^)rif{k) 
Axp{k,k) ( 1 



1 



w? - k-p\2K -Q-Q"^ 2K-Q + Q\ 



V(A2c) 
Xp 



e 

— lim 

2 



(27r)3 
d^k 



[m? — k ■ pY 



S{P - m^)rif{k) 



X ■ 



0^ (27r)3 

1 / ^'xpik,q) 



Kp{k,q) ' 

n? -k-p\2K -Q-Q^ 2K ■ Q + Q"^ , 
We discuss these contributions one by one. 



;e.23) 
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The A2a contribution 



Using Eq. ( |5.56| ) and the 5-function appearing in tlie integrand, we can write 

kxk,{^ -2m)/ 



p/(A2a) 
^ Ap 



-2e^ lim / ^^4^ 6{k^ - m^)rif{k) 



1 



1 



X- 



-k-p \2K-Q-Q^ 2K-Q + Qy' 



(E.24) 



As argued before Eq. ( |E.8|) , we can put Q = in the spinors, and use Eq. ( [4.58|) . Per- 
forming the /cQ-iiitegration, we obtain 



(2i;V-r7^0rv'"^(p) = —lim 



2e2 



m Q-,oJ (27r)3 



F{K)\ 



1 



1 



2K-Q-Q^ 2K-Q + Q^, 



(E.25) 



where the expression on the left is understood to equal the one on the right only between 
the spinors, and 



FiK) 



2El 



m 



2E 



K 



1 + 



— k-p 



/e + 1 + 



-\- k ■ p 



f-e 



(E.26) 



with ko = Ek- Since the integrand contains K ■ P, and we must set P = only after 
taking the limit Q ^ 0, the angular integrations cannot be performed using Eq. ( [E.9|) . 
So we shift the integration variable to K ± in the terms having 2K ■ Q T in the 
denominator. This gives 



{2v\P - r]^P)r[f''\p) = —lim 



2e' 



d^K Q ■ VkF 



m Q-,oJ (27r)3 2K-Q 



(E.27) 



Clearly the magnitude of Q now cancels out. The derivative with respect to K can be 
taken easily, using 



V rEk 



K 



E 



K 



K 



m? ± k ■ p 



[m? ± k ■ pY \ ^ Ef 



(E.28) 



The term proportional to P from the last derivative does not contribute because it mul- 
tiplies a factor whose integrand is odd in i^' at P = 0. Putting P = in the other terms. 
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we obtain the contribution to the gravitational mass: 



2e^ 



m 



(A2a) 



m J {2'nf2EK 



X 



- m2 [Ek - 2m dfe Ek + 2m df, 



+ 



2Ek \ Ek — m dEx Ek + rn OEk 
m „ m 



2El - m^ 



r/e 



2Ek \{EK-my {EK + mf 
2E'^ + m^ f Ek - 2m EK + 2m 



2EI 



Ek — ni 



■fe + 



Ek + rn 



/e-)|. (E.29) 



The A2b contribution 



The integral in the (A2b) term is independent of Q. So, in a straight forward way, we 
obtain 

^2 - 



m 



(A2b) 



m^ J (2'nf2EK 



{2Ef^ - m' 



EK-2m EK + 2m 

;Ie[EK) - j-^ — r~Zv^.Js[EK) 



{EK-mf 



{EK + my 



(E.30) 



The A2c contribution 



For the (A2c) contribution, first we use the expression for K\p from Eq. ( |5.37|) to find 



Alp(fc,g) = r/Ap(g^-2A;-g)(|^--2m) + (A;Agp + MA)(|^-4m) + A;A|^7p^ + A;p|^7A^, 
A^p(fc, q) = T]\p{q^ + 2k- q){^- 2m) - {kxqp + kpqxjift - 4m) - kx^lpf^: - kpi^ixf: , 

(E.31) 

dropping irrelevant 0(g2)-terms and using fc^ = m^ . In the rjxp terms, the integrand 
becomes independent of q. Thus, these terms give a regular contribution. Let us denote 
it by (A2r): 



2e2 



m 



(A2r) 



<PK 



m J {2'Kf2EK 



Ek — m 



Ek + m 



(E.32) 



The terms which appear next will be called (A2s). For these, we use the fact that 



{2v\P - v^')ikxqp + kpqx) = -2k ■ q = 2K ■ Q , 



using g ■ t> = go = 0. The Q ^ Q limit can then be taken easily, and we obtain 



(E.33) 



m 



(A2s) 



mJ {2tt)^2Ek 



Ek — m 



Ek + m 



[EM) 
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Finally, we come to the terms with three 7- matrices, which we denote by (A2t). For 
these, first we note that 



2m^ 



(E.35) 



and a similar expression with the other term. Since the ^ term vanishes between the 
spinors, we can write 



m 



{A2t) 



2e^ lim lim 



p^oQ^oJ (2tt 
ko 1 



6{k' - m')r]f{k) 
u{P')(m + m)u{P) 



we obtain 



(E.36) 



m? - k - p 2K ■ Q 

omitting the terms in the denominator since they will not contribute for Q ^ 0. Using 
the identity 



(E.37) 



+ 



2K -Q^ 



(E.38) 



between the spinors, since q -v = and ^ terms vanish. Putting this back into Eq. ([E.36|) 
and using Eqs. ( [4.58|) and (|4.59| ), we obtain 



2e' 



m 



{A2i) 



m J {2tt)^2Ek 



Ek — rn' 



Ek + m' 



(E.39) 



The sum of Eqs. { WM , O), (|K3^) , (|E;3|) and gives the total contribution of 

the A2 term, given in Eq. ( |5.68D in the text. 



E.3 The X contribution 



The part of the integral from Eq. 

d^'K 



(2vr) 



F{Ek) 



that we consider here is given by 
1 1 \ 



^2K-Q-Q^ 2K-Q + Q^, 



(E.40) 



where 



F{E)^{ff{E)-fj{E)){2E 



2 2\ 

m ] 



(E.41) 
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Shifting the variables, the integral can be written as 



Q) 



(E.42) 



(27r)3 2K ■ Q 

We have to expand the numerator to 0{Q'^) in order to obtain the integral to 0{Q^). 
Writing di to denote a partial derivative with respect to K^, 

1 /I 



The derivatives we need to use are: 



g^g^) d,d,F ± - (-g^g^g') d,d,diF . (e.43) 



EdE 



d.d^diF = {6'^K' + 6''K^ +6''K')(^^-^^ F + K'K^K'I^^ F. (E.44) 



Using 



3 



(E.45) 



within the integrand, we have 



j(/)(g) 



(2^ 



1/19 



1 g 



1 d 



1 9 



2\EdE 



3! 8 



Therefore, the quantity that we must substitute in Eq. (|5.79|) is 



/(/)(Q)_/(/)(g^O) 



3! 8 y (27r)3 I [eOE^ 



F + 



EdE 



1" v^aE 



.(E.46) 



F kE.47) 



We now use the identity 



1 d 



dK K"" {-—] F = -{n-1) dKK''-^- — 



EdE 



1 d 



EdE 



(E.48) 



which holds for > 2, so that the surface term vanishes. It is obtained by using Eq. 
( ^.122| ) and performing a partial integration. Using it repeatedly, we can rewrite Eq. 
( lOTl) as 

Q2 



/(/)(g)_/(/)(g^O) 



487r2 



\EdE 



(E.49) 



Putting this back into Eq. (|5.82|) , we obtain the total X contribution given in Eq. ( |5.83|) . 
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E.4 Alternative evaluation of the Bl and A2 contributions 



If we straightaway put g = in Eqs. ( [5.43[ ) and ( [5 .401 ), they are seen to involve singular 
expressions of the form ^r{l/x)6{x) with x = k"^ and x = k"^ — m? respectively. Using 
Eqs. ( p.6(J| ), we can obtain the following regularized form for such expressions 

(— \ 5{x) = lim 



xj tt(x^ + e'^y 

e(x - A^l 



lim lim 

€^OA^o n[{x - + e^] 

lim lim - - 



2 



.^OA^o2 7rdA2 (a;-A2)2 + e2 

We shall find that the use of Eq. (|E.5CI| ) leads to expressions for m'^^Bi) '^{A2) which 
are identical to the ones deduced earlier in this Appendix. Let us also compare Eq. ([E.50| ) 
with Eq. ( p. 611) . It is clear that Eq. ( |E.5CI| ) for x = k"^ is identical with Eq. ( |5.61|) for 
m — ^ 0. For the case of m 7^ 0, Eq. ( p.61j ) is quite convenient to use if the rest of the 
expression in which it occurs does not involve (so that d/dm? can be performed at the 
end, as in Ref. But in our problem, the rest of the expression involves not only m^, 

but m and m? as well (see Eqs. ( |E.51| ) and ( |E.52| ) below). Consequently, we have found 
it convenient to use Eq. (F30|) with x = k'^ — m . 

It may be noted that it is the form ( [E.50|) which was essentially employed by DHR 
to regularize !3^r{l/k'^)6{k'^) in T'j^p^\p,p) [^, who, however, made a detour into the 
imaginary-time formulation for this purpose. Here we first give the /«(^2) calculation, and 
then, for the sake of completeness, summarize the ^(^1) calculation. 

Putting g = in Eq. (ISliOl ), we get 

r'if\p,p) = 2e7^Av(A:,fc)^r(^^)^r(^^)5(P-m^)r,,(A;) 

d^K , 1 



/a K d f 



x6{k^ -m^ - X^)r]f{k), (E.51) 

using Eq. (|E.50| ) and also setting k"^ = + A^ as dictated by the delta function. 

Making use of the delta function and Eqs. (|4.58|) and ( |4.59D , we obtain, from Eq. 

(2t;V - r]^P)u{P)Axp{k, k)u{P) = -8k^ + 20mkl - 8m^ + 2ko{2m^ + A^ - 2mA;o)(E.52) 

119 



at f ^ = (1, 0) and p^^ = (m, 0). 

Use of Eq. (|E.52|) , and of = (m, 0), in Eq. (|E.51| ) then leads to 



(A2) 



e^(-8/i + 20mJ2 - Sm^/g + 2/4) 



(E.53) 



where 



h 
h 
h 
h 



[d%. 



1.3 


2m? + A2 - 


2mko 






2m2 + A2 - 


2mfco 


1 




2m2 + A2 - 









(E.54) 
(E.55) 
(E.56) 
(E.57) 



We have used the notation 

j [(i^/c]Afunction(A;o, A) 



hm^ I dko5ikl-El,)vfik) 



(27r)3 A^O rfA2 
X function (/co, A) 



(E.58) 



with 



(E.59) 



Performing the /cq integration first, and then removing the regulator A in Eqs. ( [E.54D - 
( IE.57D (it is useful to make note of 



A^o dX 

for this purpose), we arrive at 

1 r d^K 



!i™;3T^/e,e-(^i^A) = OF fe,e{.EK) 



h 



Srr? J (277)3 
1 



^EkOEk 



(EM) 



Ek — ra 



+mE 



K 



dfe 



+ 



Ek + m 



+m\ 



y Ek — m dEx Ek + m SEk 
d'K 

Of. 



(E.61) 



(27r 



m 



1 E, + m Ek 

tjK \ ~ ^ + ^ 



f-e{E 



K 



1 dh 

Ek — m OEk Ek + m OEk 



(E.62) 
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1 



m 



(27r)3 
1 



1 

Ek 



3 (UEk) + MEk)) 
1 df, \ 



E]^ \ Ek — rn dE^ Ek + rn dE^ , 



(E.63) 



1 r d'^K 1 d 

Substitution in Eq. ( |E.53D finally gives ^(^2)' which is found to be the same as the result 
given by Eq. ( p. 681 ). 

For T'j^p^\p,p), given by Eq. ( |5.43|) , similar steps are to be carried out. Thus the 
analogues of Eqs. (|K5l|) , (|K5^ ) and (WM) are 



lim 



d 



X 



^ -die - X^)vj{k), 



-2p-k 



(E.65) 



{2vy,Vp - 7]xp)u{P)Y{i> + m)j'^u{P)C^,xp{k, k) = 2X^{ko - 2m) + Amkl - 8fc3,(E.66) 



m'(Bi) = e\2Ji + 4m J2 - 8J3). (E.67) 

Here 

The notation [c/^/cJa is now defined by Eq. ( [E.58| ) with rjf —>■ rj^, and the notation Ekx by 
Eq. ( [E39|) with m = 0. 

The evaluation of the above three integrals is straightforward. Each one reduces to 
the standard integral 

dKKf,{K) = —-, (E.71) 
6 



while J'? also involves 



00 

^^^ \ePK_^y =-9pl dKKf,{K) (E.72) 
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We then have Ji = -TV24m, J2 = -T'^/ASm^ and J3 = T'^/ASm, so that the value of 
m|^j^^ is the same as the one stated in Eq. (|5.63| ). 

Finally, we comment that the regularization presented above is not useful for evaluat- 
ing the X-contribution, as we need the 0{Q'^) contribution from the fermion loop in that 
case, and not the the value at g = 0. 
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APPENDIX F 



Thermal contribution to the 
equation for the effective mass at 
high temperature 



In this appendix we arrive at Eq. ( IETTD , which expresses the high T contribution to 
the equation for the effective mass, by exphcit calculation of the expressions given in Eqs. 
(PI), ( 16:631) , (IHD, ( lOeD and 



We begin with M^q^, given by Eq. ( |6.62|) . After the po integration is carried out 



the integrand is a function of P (the magnitude of the three-momentum) alone. Then 
integrating over the angles trivially, we are left with a one- dimensional integral. On doing 
some rearrangement, and making use of the identity 

fsiP) - friP) = 2fB{2P) (F.l) 



(here and elsewhere in this appendix, /f(-P) and /b(-P) denote ( |6.11[ ) and ( |6.15| ) with \p, 
replaced by P), we obtain 



ol 



M' ' 



Now the first part of the R.H.S. of Eq. ( [l''.2| ) can be evaluated exactly by using 



/ dPPMP) = (F.3) 

JO D 

dPPMP) = (F.4) 



The second part of the R.H.S. of Eq. (|F.2|) involves the integral 
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Here S^r denotes the principal value. The functional dependence of / on TjM was 
determined on a computer for large T /M, yielding 

1 T 

— In—. (F.7) 

Note that this logarithmic behaviour is in accord with the ultraviolet behaviour of (|F.5|) 
without the distribution function. Thus at high T 

M'^g' = eH— + —Ari-) (F.8) 
V 8 87r2 ^ ' 

neglecting terms independent of T. 

Next we turn to M'^ J^, given by Eq. (|6.63| ). Integrating over po and the angles, we 
are led to 



M' ' 



rdPPMP)Re I 

dPPfF(P)Re : — . (F.9) 

327r2 Jo ^ ' (p2 _ _ V ; 



Moving on to M'^^^'^, integration over the angles in Eqs. ( |6.65D , ( |6.66D and ( |6.67| ) lead to 



M' 



r ^^^^^) p2 (F.io) 

M'am = - \ / dPP^^ ^Re -. (F.12) 

Using the identity we then arrive at 



2 /\^2 





f°° 1 



/■oo 

/ ciPP/B(P)Re 

JO 



47r2 Jo ' ' ' (P2 _ M2 -2ie)2' 

(F.13) 

While the first part again involves ( [!''. 5| ), the second part, being ill-defined, is regularized 
in a way similar to ( |6.3CI| ): 

/; dP PHP) Re ^p2_M^_,,,y - 1- ^ /; PMP) Ke p. _ ,1 _ • 

(F.14) 
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The integral on the R.H.S. of Eq. ( |t'\14 ) is similar to / in Eq. ( |1'\5| ) and equals — | ln(T/A) 
for high T. So the second part of ( |F.13| ) equals —^e'^M'^/lQir'^ which can be neglected as 
it is independent of T. Thus we are left with the first part of ( |F.13| ), so that 



In M^^g(jjj)^, we encounter the integral 

1 p^/^ 1 

where Mi = M/2. Converting J into a function of T /Mi alone (as we did in the case of / 
in Eq. ( [!''. 6| )), the functional dependence was found out for large T /Mi using a computer, 
which gave J ^ — i. This constancy again agrees with the expectation from the argument 
involving the ultraviolet behaviour. Alternatively one can note that 



J = Tt^L dPP (F.17) 



Then using our knowledge of / (see Eqs. ( [b'.5| ) and |i^'.7| )), we obtain J ^ — ^ at large T. 
Thus equals ^c^M^/IGtt^ and can be neglected. 

Therefore in the high-T limit, M'^ receives contribution from ([F.8|) and ( |F.15| ), thereby 
leading us to Eq. ( |6.71| ). 
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